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ON THE MOTION OF A SPINNING SHELL* 


BY 


K. L. NIELSEN f (Louisiana State University) AnD J L. SYNGE (The Ohio State University) 


1. Introduction. Next after the problem of the motion of a particle in a resisting 
medium, the problem of the motion of a spinning shell is the major problem of ex- 
terior ballistics. Many crude treatments have been given, but the problem was first 
discussed exhaustively by Fowler, Gallop, Lock and Richmond.':? Reference may also 
be made to treatments by Cranz* and Moulton.‘ 

An exact treatment of the motion of a spinning shell as a hydrodynamical problem 
is obviously out of the question. The problem must be treated aerodynamically. This 
means that the forces exerted on the shell by the air must be regarded as dependent 
only on the instantaneous motion of the shell. The connection between the aero- 
dynamic force system and the motion cannot be deduced logically. It must appear in 
the mathematical theory as a hypothesis, preferably supported by experimental ob- 


servations. 


But although mathematical theory cannot supply the aerodynamic forces, it does 
give us some information about them. Two basic ideas are important here. 
First, the shell has an axis of symmetry. This fact has been used in all existing 


theories. 

The second idea is a little more subtle. It concerns the connection between the 
position of the mass center (or center of gravity) of the shell and the aerodynamic 
force system. In one manner of speaking, there is no such connection. For two shells, 
moving with identical motions but with different mass-distributions, the aerodynamic 
forces are the same. But we cannot introduce the aerodynamic force system into the 
mathematical argument without expressing that force system mathematically as a 
force and a couple (or something equivalent). To do this, we must use a base-point, 
* Received January 22, 1946. This paper was written in 1942, when one of the authors (J L. S.) was 
at the University of Toronto. It was issued asa restricted report in January 1943 by the Ballistic Research 
Laboratory, Aberdeen Proving Ground, with permission of the National Research Council of Canada. 
Later work by other authors, issued in restricted reports, has improved on some of the theory, but it has 
been thought advisable to publish the paper in its original form. 

ft On leave with U. S. Naval Ordnance Plant, Indianapolis, Ind. 

1 R.H. Fowler, E. G. Gallop, C. N. H. Lock, H. W. Richmond, The aerodynamics of a spinning 
shell, Phil. Trans. Roy. Soc. London (A) 221, 295-387 (1920). 

2 R. H. Fowler, C. N. H. Lock, The aerodynamics of a spinning shell, Part 11, Phil. Trans. Roy. Soc. 
London (A) 222, 227-249 (1921). 

*C. Cranz, Lehrbuch der Ballistik, J. Springer, Berlin, 1925, p. 358. 

*F. R. Moulton, New methods in exterior ballistics, University of Chicago Press, Chicago, 1926, 
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and reduce the force system to a force at that base-point, together with a couple. 
It is well known that, for a given force system, the force is independent of the base- 
point, but the couple is not. 

Also, to describe the motion of the shell mathematically, we must use a base- 
point. The motion is described by the velocity of that base-point and an angular 
velocity. The angular velocity is independent of the choice of base-point, but the 
velocity is not. 

Now it is natural to use the mass center as base-point. If there are two shells, 
S; and S:, with mass centers O; and Oz, we may use OQ, as base-point for S; and O:2 as 
base-point for S2. Suppose that the two shells are of identical geometrical form (but 
O; and O2 are not geometrically corresponding points) and that their motions at the 
instant are the same. (This means that geometrically corresponding points have equal 
velocities; the velocities of O; and O2 are not the same.) Then the force systems on the 
two shells are the same. But the moments about Q; and Q: are not the same. 

If we set out to formulate aerodynamic laws, using the mass center as base-point, 
we must exercise great care. We must ensure invariance with respect to shift of mass 
center. We must make sure, in the case described above, that when we apply our law, 
first to S,; and then to S:, we get equivalent force systems. 

Unfortunately, Fowler et al.' paid no attention to this fact in formulating their 
aerodynamic laws (pp. 302-305), although they draw attention to the necessity for 
invariance (p. 305), and in fact make use of it. By considering a special case, it is 
easy to see the fallacy in their basic laws. 

Consider the two shells described above. Let the velocity of O; be directed along 
the axis of the shell, and let the shell have an angular velocity represented by a vector 
perpendicular to the axis (plane motion). The yaw is zero, and the effect of the air 
is a drag along the axis. But now consider Sz. On account of the angular velocity, the 
velocity of O2 is not along the axis; there is a yaw, and hence a cross wind force in 
addition to a drag. It is easy to see that the force systems on the two shells are not 
equivalent, as they ought to be since the motions are the same. 

Thus the theory of Fowler et al. contains a logical contradiction. It is very diffi- 
cult to discuss critically a theory containing a logical contradiction, for from incon- 
sistent hypotheses we may arrive almost anywhere (at 1=0, for example.) It may 
well be, however, that the logical contradiction does not invalidate the physical con- 
clusions of their paper. In the example given above, the yaw of S:; may well be very 
small indeed in cases of practical interest, and the logical inconsistency may be no 
more serious than that involved in writing 7 =3.14. Used in one way, this statement 
leads to 1=0; used in another way, it leads to important practical results. 

Nevertheless it is sound policy, in building up a theory in applied mathematics 
to make it logically consistent as far as possible. In the present paper we shall take 
care to state the aerodynamic laws in such a way as to avoid logical inconsistency. 

Apart from the thorough treatment of the theory of the aerodynamic force system 
in sections 3 and 4, the following features of the present paper may be summarized 
here. 
The exact equations of motion of the shell (independent of any aerodynamic hy- 
pothesis) are given a very compact form in (2.6). In section 5 it is shown how the aero- 
dynamic functions may be found from high frequency photographs of a shell. Such 
observations should provide the ultimate test of the validity of the aerodynamic 
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method. In view of the success of the cruder jump card method of Fowler et al., 
it seems probable that the aerodynamic hypothesis is valid, and, if so, the pro- 
posed method of observation should give us all information required concerning 
the aerodynamic functions. 

There are three conditions for the stability of a spinning shell (section 7), but 
they are too complicated to interpret in the general case. If Magnus effects are 
absent (section 8), they become much simpler, and in fact there is then just one 
stability condition (8.19). In this condition the effect of the position of the mass 
center is shown explicitly. The condition is stronger than the usual condition 
(8.13b) based on the stability factor; a shell which is considered stable on the 
basis of the usual condition may in fact be unstable. We are very much indebted 
to Professor E. J. McShane for his critical comments on this paper in its original 
form. He has informed us that the existence of second stability condition, 
stronger than the usual one, has already been pointed out by R. H. Kent (Re- 
port No. 85, Ballistic Research Laboratory). This condition is implicit in the 
paper by Fowler et al. (1.332, equation 3.6234, and 4.12); this is discussed in 
section 10, where their method is brought into line with the more general method 
of the present paper. 

Some well known facts are confirmed by theory in section 9. For a stable 
shell, after the oscillations have been damped out, the axis of the shell always 
points above the trajectory and to the right if the spin is right-handed. The 
phenomenon of trailing is explained; the axis of the shell turns downward at a 
rate approximately equal to the rate of turning of the tangent to the trajectory. 

Drift also is discussed in section 9. A general condition (9.17) is obtained for 
standard drift, i.e., drift to the right for right-handed spin. When we specialize 


to subsonic velocity and flat trajec- q 
tory, this condition simplifies to (9.20). 

When the numerical values of Fowler Ak 

et al. are inserted, this inequality is 

liberally satisfied, so that the present e 


theory is in agreement with the ob- 
served facts. 

2. Exact equations of motion. We 
shall now develop the equations of mo- 
tion of a shell in convenient form. No as- 
sumption is made here regarding the fe) 
aerodynamic forces, and the only as- 
sumption regarding the shell is that it 
has an axis of dynamic symmetry (i.e., 0 
the momental ellipsoid at the mass cen- 
ter is a spheroid). Thus our equations 
would apply, for example, to a homo- J 
geneous projectile of square section or y 
to a bomb with three or more fins, ; 
placed symmetrically. 

We shall use the following notation, Fic. 1 
the motion being referred to a Newtonian reference system: 








na 
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O=mass center of shell, 
m =mass of shell, 
A, C=transverse and axial moments of inertia at O, 
q= velocity of O, 
w@ =angular velocity of shell, 
h=angular momentum of shell about O, 
F = vector sum of aerodynamic forces acting on shell, 
G=moment of aerodynamic forces about O, 
F’ = weight of shell. 
Then the equations of motion are 


mq =F + F’, h =G. (2.4) 


We introduce a right-handed unit orthogonal triad, i, j, k, fixed neither in 
space nor in the shell (Fig. 1). We take k along the axis of the shell, and i, j 
perpendicular to k, but the final choice of i, j is deferred for the present. LetQ 
be the angular velocity of the triad. 

We may now resolve the vectors as follows: 


= uit w+ wk, 





q ] 
a= ait wot ol, | 
Q= Ait Ajt Ak, | 
bh = sod + Aug + Cok, | (2.2) 
> id Fad Fe | 
C=. G+ Gis ce | 
Fo = Ff/i+ Fij+ Fik. | 


Clearly 2)}=@;, Qe = we. 
In scalar form the equations of motion (2.1) then read 


m(u — 23 + wwe) = Fi + Fi, ) 
m(i — wo, + uQs) = F2+ Fi, } (2.3) 
m(w — uwe + v0;) = Fs +F3, } 


A (w, a w2Q3) + C302 _ Gi, ) 


Il 


, . y 4 
A (we + w1Q3) = Cw 30 = Go, (2.4) 
Cas = G3. ) 
It is now convenient to introduce complex variables. We write 
ut+tiw =i, 
Wi + 1W9 = 7, 
F, + iF; =F, } (2.5) 


G,; + iG, =G, 
FY + iFg = F’. } 
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We multiply the second equation of (2.3) by i and add it to the first, and deal 
similarly with the equations (2.4). Thus we reduce the equations of motion to 
the form: 
E+ itQ3 — iwn = (F + F’)/m, 
4 + inQs — Com =G/A, (C’ =C/A), 
(Fs + Fs)/m, 
3 = G;/C. 


(2.6) 


W— Uwe + Vw 


These equations are exact; no approximations have been made. 

3. The general aerodynamic hypothesis. What is here set down is probably a little 
more general and explicit than previous statements about aerodynamic force 
systems. There is no implication that the hypothesis is physically accurate in 
all cases. All we can hope is that deductions from these assumptions lead in suit- 
able cases to results in fair agreement with observation. But it seems best to 
make the hypothesis mathematically clear. 

First we consider a fluid, at rest or in motion. We are not particularly con- 
erned with the properties of the fluid. The important thing is that it defines 

(i) a scalar field of density p; 

(ii) a scalar field of local sound velocity c; 

(iii) a vector field of velocity W. 

This last field defines two other vector fields, vorticity (V=1/2 rot W) and ac- 
celeration (a=dW/dt). 

Usually in ballistics we deal with the static case in which W=0 and p, c are 
functions of height only. A more accurate model is that in which W is horizontal, 
but in different directions at different heights to allow for changes in the direc- 
tion of the wind with variation of height. 

Now suppose we wish to investigate the motion of a solid through this fluid. 
To treat the problem adequately we should of course consider the disturbance 
produced in the fluid by the solid. But we do not do this. We use the fluid merely 
to compute from its undisturbed motion the aerodynamic forces acting on the 
solid. 

Let O* be the centroid of the solid, i.e., the position its mass center would 
occupy were the solid of uniform density. Let the motion of the solid be de- 
scribed by the velocity q* of O* and the angular velocity o*. 

The basic hypothesis is then as follows: 

Aerodynamic hypothesis: The aerodynamic force system exerted on the solid 
by the fluid consists of 

(i) an aerostatic force; 

(ii) an aerokinetic force system. 

The aerostatic force acts at O* and equals 


pV o(a — P) (3.1) 


where p is the density of the fluid at O*, V> is the volume of the solid, and P 
is the body force per unit mass acting on the fluid at O*. (Note that ifa=0and P 
is gravity, this is simply the Archimedean buoyancy.) The aerokinetic force sys- 
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tem is represented by a force F* at O* and a couple G*; these are functions of p 
and ¢ at O* and of the vectors 


¢- FF, =. (3.2) 


If q*=W and w*=V, then F*=0 and G*=0. 

Henceforth we shall assume W=0, and so F*, G* depend only on ap, c, q*, w*, 
while the aerostatic force is —pVoP. If we were discussing the aerodynamics of a 
dirigible, the aerostatic force would be very important. For a shell it is quite 
trivial and we shall omit it altogether. 

Thus for our purposes the aerodynamic force system consists of the force F* 
at O* and the couple G*; they are functions of p, c, q*, and w*. 

It will be observed that our base-point O* has been chosen in a definite way 
with respect to the geometry of the solid, and not with respect to its mass-dis- 
tribution. This frees our laws from the objection raised in the Introduction to 
the laws of Fowler et al. 

It is to be noted that it is by no means essential to select the centroid as 
base point. But it is least confusing to choose, once and for all, a point simply 
related to the geometry of the solid, and the centroid seems the most natural 
point to take. 

4. The aerodynamic force system for a shell with an axis of symmetry. We now 
consider a shell with an axis of aerodynamic symmetry. By this we mean that 
its exterior is a surface of revolution. We might proceed for the present without 
introducing the mass-distribution of the shell, but it seems simpler to proceed 
at once to the case of complete symmetry. We shall therefore suppose that the 
shell has a common axis of aerodynamic and dynamic symmetry. All that is 
stated in section 2 is then valid and we shall use the same notation. 

The mass center of the shell is at O and its centroid at O*. Let us write 


— 
OO* = rk, (4.1) 
and 
q* = velocity of O*, } 
‘ } 
w* = angular velocity of shell, | (4.2) 
F* = vector sum of aerodynamic forces, sis 
G* = momentofaerodynamicforcesaboutO*. } 
Then . 
q* = q + o X rk, o* = o, 
(4.3) 


F*=F, G*=G+FX rk} 


In the notation of (2.5) with asterisks attached to the symbols referring to O*, 


we have in consequence 


gP = § — Gry, wt = w, | 
| 
n* = Un w3" = W3, | 
; (4.4) 
F* = F, I 3" ao id ae 


G* = G — irF, G;* = G3. 
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Now F* and G* depend on q* and a@*. It follows from the aerodynamic sym- 
metry that if the pair of vectors q*, w* is given a rigid body rotation about the 
axis of symmetry, then the pair of vectors F*, G* is also rotated rigidly about 
the axis through the same angle. Hence the following ten scalar quantities are 
unaltered by such a rotation: 

F3, G3, 


u*F * + o°F * o*F * = u*F*, u*G;* + v*G.*, o*G;* —u*G*, ) 4 5) 
‘4 ° 


wf F* + wr, we F * — oF e, w*Gy* + w8Ge*, wil * — wl". | 


But, to within such a rotation, the vectors q*, w* are determined by the quanti- 


ties 
W, Ws, UF? + oF? wy? + wok? uFw* + vtwo*, u*wo* — v*w;*, (4.6) 
between which there exists the identity 
(u*? + y**)(wi*? + wo¥®) — (u*wi* + v*wo**)? = (u*wo* — v*w,*)*. (4.7) 


Therefore the quantities (4.5) are functions of the quantities (4.6); in fact, for 
a shell of given size and shape, (4.5) are functions only of (4.6) and the air scalars 
gp. ¢ ato". 
We now write 
u*F* + o*Fo* = 5, v*F,* — u*F.* = So. (4.8) 
Multiplying the second equation by 7 and subtracting it from the first, we get 
E*F* =3$,;-— iSe, (4.9) 


the bar denoting the complex conjugate. Dealing similarly with the other quanti- 


ties in (4.5), we see that 
cope EG* ) 
aa Aiea (4.10) 
Fr, 7*G*, | 


are complex functions of the real quantities in (4.6). 
We cannot proceed further without an additional hypothesis. We shall as- 


sume that 


F}, F.*, Gi, Ge 
are linear functions of 
u*, v, w;"*, wo* 


This is certainly a reasonable assumption when the latter quantities are small. 
We can then write 

F* = ayu* + agu* + Bywo:* + Bows", ) (4.11) 

. f ' 

G* = yu* + you" + bywr* + Sowe*, | 


where the eight complex coefficients are functions of w, w;, p and c. When we 
form the quantities (4.10) and use the fact that these must be functions of the 
quantities (4.6), we find ag=ia;, Be =18;, etc., and so 

F* = g*P* + n*O*, 

G* = ¢* p/* + n*Q’* 

s ’ 


(4.12) 


where P*, Q*, P’*, Q’* are complex functions of w, ws, p, c. 
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The components F3, G3; are functions of the quantities (4.6). We shall assume 
that they are functions only of w, ws3, p, c. This also is a plausible assumption when 
u*, v*, wi*, we* are small. 

To sum up: There are ten real aerodynamic functions of w, ws, p, Cc, contained 
in the set 

-. O*, Pe: oO”, F;*, G;*. (4.13) 

Let us see what these assumptions amount to in the case of a shell in a wind-, 

tunnel. We think of the shell as moving and the air at rest. We put 

o* = Q, a," = wo” = w;* = 0, 
and (4.12) gives 
Fi + iF = u*P*, G;* + iG* = u*P"™. 


In this simple case we must have, by symmetry since w;* =0, 


Fs = G;* = G3* = 0, 
and so we have 
FY = u*P*, iG3* = u*P'*. (4.14) 
It is easy to see that these equations imply that (for small yaw), the cross wind force 
and the moment are proportional to the yaw. This is the usual assumption. 
We now pass from the centroid O* to the mass center O by the transformation 
(4.4). We get for the force system F, G on the shell 


Fz Fy + iF > = EP + nQ, F; aed F;*, } 


7 ; : 5 : (4.15) 
G= Gy + iG2 = Ep’ oe nQ’, G3 = G;*, ) 
where P, Q, P’, Q’ are complex functions of w, ws, p, c, given by 
P = P*, QO = Q* — irP*, 
(4.16) 


P’ = P* + irP*, Q’ = Q* — irP’* + ir(Q* — irP*). 


This gives the transformation of the aerodynamic functions when we pass from the 
centroid O* to the mass center O. Actually this is the transformation for passage from 
any base-point to any other, provided of course that both lie on the axis. 

To show the real and imaginary parts of the aerodynamic functions, we shall 
write (with similar equations in asterisked form) 


P = P, ss iPs, Q - Q—”; + 102, | 


Pi = Pi{* — rP#*, 
Pi = Ps* + -rP;{, 
Or = Qi* + rPz* + r(— OF + rP3*), | 
Of = Of* —rP{*4+70*+4+9P*). | 


c 


; : (4.17) 
P’ _ PY + iPS, QO’ — OY + iO? . J 
The transformation (4.16) then gives 
P, = Pr, | 
P, = P#, 
’ | 
0, = O* + rP*, 
Qo = Qs* — rPi*, | 
if 4.18) 
| (4.18) 
| 
| 
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The method used above for the resolution of the aerodynamic force system is not 
the usual one. Three important vectors are involved: k the axis of the shell, q the ve- 
locity of the mass center, w the angular velocity. In resolving vectors, it is necessary 
to pick out one of these three as a fundamental vector and build a basic triad on it. 
The traditional plan is to pick out q as fundamental and take k as a secondary vector, 
so that q and k together give one of the planes of the basic triad. Resolution of F 
along q and perpendicular to q in this plane gives the usual drag and lift. However 
convenient this may be for wind-tunnel work in which q is fixed while k is altered, 
it certainly appears less convenient than the method of the present paper for a simple 
mathematical formulation of the problem of the spinning shell. There is a further 
objection to the usual plan; the direction of q depends on the mass center. 

The conventional terminology does not suit the present resolution. The following 
is suggested. The asterisk indicates that the centroid is used as base-point. The same 
notation without asterisks refers to the mass center. 


u*i + v*j = cross velocity, 


wk = axial velocity, 





| 
(4.19) 
wi + wej = cross spin, | 
w3k = axial spin. } 
P;*| | = cross force due to cross velocity (—), 
P# | | = Magnus force due to cross velocity (+), | 
Q,*| n*| = Magnus force due to cross spin (+), ? (4.20) 
Q#| n*| = cross force due to cross spin (+), | 
F; = axial force (—). 
P{*| &| = Magnus torque due to cross velocity (—), 
P{*| t*| = cross torque due to cross velocity (—), | 
Q{*| n* | = cross torque due to cross spin (—), (4.21) 
Q4* | n*| = Magnus torque due to cross spin (+), 
G3; = Magnus axial torque (—). 


It is a consequence of symmetry that where the word “Magnus” is included above, 
the quantity in question changes sign with 3; where the word “Magnus” does not 
occur, the quantity in question does not change sign with w;. For uniformity, we have 
called the axial (viscous) torque “Magnus”; there is justification for this in the fact 
that it is the viscous torque that sets up the circulation which is responsible for the 
other Magnus effects. The signs in parentheses indicate probable signs of the various 
quantities when w; is positive, assuming a center of pressure in front of the centroid. 
Since | 


| &*| = g* sin (q*, k), | n*| = w sin (o, k), (4.22) 


it is clear that the usual sine law of variation is implicit in (4.20), (4.21). But since 
we suppose the angles in question to be small, the sine, tangent and circular measure 
are not distinguishable. 

It is convenient to introduce positive dimensionless aerodynamic functions, as is 
done by Fowler et al. So we write, paying attention to dimensions and signs, 
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P* = — pa*wfi*, PP? = parwsfo*, | 
O* = patwsg,*, Q.* = pa wes", | 4.23 
P{* = — pa los fi*, Pi* ‘ie alae pawfs*, ( . 23) 
Oi* = — pa'wg;*, QO7* = pa®wsee *. 


Here p is the air-density and a the radius of the cross section of the shell. The func- 
tions (f*, g*) depend certainly on w/c, and possibly also on aw;/c and the Reynolds 
number. The above equations may be regarded as definitions of the eight aerodynamic 
functions (f*, g*), which are analogous to the fz, fw, etc. of Fowler et al. To the above 
equations we may add 

F; = — pa*w’ fs, G3 = — patwwsg3, (4.24) 


where f; and g; are dimensionless; fs; is the usual drag except for the slight difference 
that we resolve along the axis of the shell and treat w as basic instead of g*. 
As the notation is necessarily somewhat complicated, let us summarize as follows: 
Askerisked quantities refer to the centroid, unasterisked to the mass center. 
The aerodynamic force system is denoted by 


F* = F;* + iF *, G* => G;* + iG, F3, G3. 
There are ten real aerodynamic functions contained in the set 
aaa C*, Fr, ?, F;, G3, 


and these may be expressed in terms of the ten positive dimensionless aerodynamic 

functions 

- ee 
v4 


+ of pk x rr 
Few Je"; £1 82", 


t* t* 
J2°5 Sas 82 


, Ts, §3.- 


The same notation may be used with reference to the mass center, but since the aero- 
dynamic force system has nothing to do with the mass center as such, the asterisked 
quantities are the more fundamental. If we wish to pass from O* to O, we must trans- 
form by (4.18) and (4.23). Thus f*¥ =fi, fe* =fo, fs* =fs, gs*=gs, but the other functions 
change. 

One more notation will be introduced for convenience in (6.4). 

It is clear from (4.20), (4.21), (4.23) that if the dimensionless aerodynamic func- 
tions (f*, g*) are constants, we have the following proportionalities, 6 denoting the 
small yaw: 

cross force due to cross velocity « w’é, 
cross torque due to cross velocity x w6, | 
axial force <« w’, | 


axial torque < ww3. | 

The first three of these are in agreement with experiment for subsonic velocities—the 
effects vary as the square of the velocity. The last (axial torque) requires comment. 

The form of G; in (4.24) agrees with Fowler et al., but one may ask why (apart 
from the theory of dimensions) the factor w should be present. The following is a pos- 
sible explanation. The rotation of the shell generates a rotating wake. If this wake has, 
throughout, the same spin as the shell, it has angular momentum $7pa‘w; per unit 
length. In unit time a length w of wake is generated, and so, by the conservation of 


1946] ON THE MOTION OF A SPINNING SHELL 211 


angular momentum, the rate of loss of angular momentum of the shell is 
= G3 = 3 rpa'wws. 
This argument not only confirms the form G; of (4.24); it gives 


gs = $m. (4.26) 


A crude argument of this sort must be accepted only provisionally in the absence of 
experimental check. 
5. Determination of the aerodynamic functions by observation. Fowler et al. 








stressed the importance of avoiding the simple empirical assumptions previously em- 
Z, 
i> | 
0 Z) 
I 
O, ‘ Y, 
° 











Fic, 2 


ployed. As in the case of the drag function, it is necessary to determine the aerody- 
namic functions experimentally. What follows is a refinement and generalization of 
the jump card method of Fowler et al. Unless there are technical difficulties, or unless 
the basic aerodynamic hypothesis is wrong, the following method should yield all the 
aerodynamic functions quite simply, except perhaps g;, and no doubt a method could 
be devised for it also. 

Let a shell be fired horizontally and observations made of it not long after it leaves 
the muzzle. These observations consist of high-frequency photographs, one set of 
photographs being taken vertically and the other set horizontally from the side. These 
photographs show successive positions of the shell at short intervals of time. 

We now turn to the exact equations of motion (2.6). There is some indeterminacy 
in these because we have not yet chosen the vector i definitely. Let us choose it in 
the vertical plane through the axis of the shell (k), pointing downward (Figure 2). 
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Ne 
Ne 


Then 
fF’ = mg cos 8, (S. 2) 
and the first two equations of (2.6) may be written 
F = m(é + iQ; — iwn) — mg cos 8, | 
G = A(q + inQs — iC’ws). j 


These equations are exact. We may put cos @=1, since the axis of the shell is approxi- 
mately horizontal. Then 23;=0 by (6.2). 

Now m, A, C’ are known for the shell; w may be found from the observations or 
otherwise (muzzle velocity), and ws; deduced from the rifling. To find £&, 7 as functions 
of t, it is merely necessary to measure on the photographic plates the linear displace- 
ments of the mass center and the angular displacements of the axis of the shell, cor- 
responding to the short intervals between successive photographs. Smooth graphs 
might be made showing u, v, wi, We as functions of ¢ or the complex quantities £, 7 
might be plotted on an Argand diagram with the values of ¢ marked in. In any case 
it should not be difficult to obtain £ and 7 also as functions of ¢ from these graphs. 

When these functions are inserted in the right-hand sides of (5.2), we have F 


and G as functions of t. By (4.15) we have 
tP + 70 =F, iP’ + nQ’ =G. (5.3) 


If we use two values of ¢, each of these equations yields two complex equations, and 
from them P, Q, P’, Q’ can be found. Here we have a good test of the aerodynamic 
hypothesis, for the values of P, Q, P’, Q’ should be independent of the particular in- 
stants chosen. 

It may be advisable, as a refinement, to allow for the decrease in w between the 
two instants in question. This can easily be done from our knowledge of the drag 
function. 

By repeating the experiment on the same shell, but using different muzzle veloci- 
ties and riflings, we obtain P, Q, P’, Q’ as functions of w and a. 

The next step is to transform from the mass center to the centroid. This is done 
by (4.16), and we obtain P*, Q*, P’*, Q’* as functions of wand w;. Finally, the dimen- 
sionless aerodynamic functions (f*, g*) are found from (4.23). 

It should be stressed that these last functions are characteristic of the form of the 
shell and completely independent of the mass distribution. Indeed, to a certain ex- 
tent they will be independent of the size of the shell, but this must be accepted with 
caution. 

6. Plan of solution and partial linearization of the equations. We now introduce 
fixed axes Opxoyo20, Oozo being directed vertically upward. Let @ be the inclination of k 
to the horizontal (Figure 2), and ¢ the inclination of the horizontal projection of k 
to Ooxo. We have already made the vector i definite in section 5. We have then 


F’ = mg cos 0i — mg sin 6k, ) 


. is a } (6.1) 
Q = — ¢cos Ai — 6) + ¢ sin 6k. | 


Hence ’ ’ 
n = — (¢cos 6+ 2), Q3; = ¢ sin @. (6. 2) 
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We substitute from (4.15) in (2.6), and the equations of motion become 
£+ iQ; — iwn = EX + nY + g cos, | 
7 + inQs — iC’w3n = EX’ + nV’, | 


, (6.3) 
WwW — Uwe + 20; = F;/m — g sin 0, | 
és = G;/C, 
where 
X = P/m, Y = Q/m, X’ = P’/A, Y’ = Q’/A. (6.4) 


If we substitute from (6.2) for n, Q; and regard X, Y, X’, Y’, F;, G3; as known func- 
tions of w, ws, p, c, we have six real equations for the dependent variables u, v, w, 
6, d, ws. But unless we assume p, ¢ to be constants, we must bring in further equations. 
Let us assume them to be functions of height (zo) only. By resolution of velocity we 
have 
Xo + ivo = (u sin 6+ iv + w cos Oe, | 
ir } (6.5) 
o= — ucosé+ wsin 8. 


t- 


When the last of these equations is associated with (6.3), we have seven real equa- 
tions for seven unknowns, namely, those stated above and zo. When they have been 
solved, the trajectory of the mass center is given by (6.5). 

We now make the following two assumptions: (i) the vertical plane through the 
axis of the shell turns slowly; (ii) the angle of yaw is small. The first assumption 
implies that ¢ and hence Q, is small; the second implies that £/w is small. On account 
of the smallness of 23; we reject the second terms in the first two equations of (6.3), 
and on account of the smallness of £/w we reject the second and third terms in the 
third equation. 

Our partially linearized equations now read 


— — iwn = EX + nY + g cos8, |) 
—— sn = x* 4. Y’, 
me 4 6.9 
w= F;/m — g sin @, | 
3 = G;/C, ) 
where : 
n = — (@ cos 0 + #6). (6.7) 


7. The stability of a spinning shell. In discussing rapid oscillations of the shell, 
we treat w and ws; as constants in the first two equations of (6.6). Consequently 
X, Y, X’, Y’ are constants. In rapid oscillations differentiation with respect to ¢ 
greatly increases the importance of a term. Hence we shall treat cos @ as a constant in 
the first equation of (6.6); the term corresponding to a small change in 6 will be neg- 
ligible in comparison with the terms in 7. 

We have then linear equations with constant coefficients, which have solutions 


of the form 


a 
s 


A yeu! + A 2e%! + A3, ) 
Byes! + Boe + Bs, | 


(7.1) 


II 


U] 
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where aj, a2 are the roots of the equation 


a? — (iC’ws + X + Ya + i(C'wsX — wX') + X¥'-— XV =0, = (7.2) 
and 
As = — — cos 0(iC’ws + Y’), | 
E 
g (7.3) 
B; = : cos 6-X’, 
E 


E = i(C'w3X — wX’) + XY’ — X’Y. 


The condition for stability is that both roots of (7.2) should have non-positive real 


parts. 
If we write 
Ki= X,+ Fi, 
Ky = Clos + X2o+ ¥?, | (7.4) 
Kin «Cat + atl + EE — BE = X67 + KEY. | 
Ke= CoarX,— wX{i + X¥:¥i + X.¥{ — X{¥2—Xi¥», | 
then (7.2) becomes 
a? — (Ki + iKs)a + (K3 + ik, = 0. (7.5) 
The condition for stability may be written 
Ki+¢cosx <0, (7.6) 
where ¢, x are defined by 
p* = (Ki — K;— 4K) +4(KiK2—2K), {$20 
¢ sin 2y = 2(KiK2 — 2K,), (7.7) 
colke R,-Ki- ks, -iwSxS ie. 


It is immediately evident that there is instability if K,>0. If K; SO, then the con- 


dition (7.6) is equivalent to 


K; =f cos x (7.8) 
or 
2K; = ¢ (1 + cos 2x). (7.9) 
On substituting for ¢? cos 2x from (7.7), this becomes 
ve " x, + 4K;2 C. (7.19) 


Thus there is instability if K, $0, K? + K? +4K; <0. If Kis0, K? +K? +4K3;20, the 


condition (7.10) is equivalent to 


(Kit Ke+4Ks) 2¢, (7.11) 

















1946] ON THE MOTION OF A SPINNING SHELL 215 


and, on substitution from (7.7), this becomes 
KiKs + KiK2K,— Ky 0. (7.12) 


To sum up, the motion of the shell is stable if, and only if, the following three condi- 
tions are all satisfied: 


K, < 0, (7.13a) 
Ki + K: + 4Ks 2 0, (7.13b) 
KiK; + KiK:K, — K,2 0. (7. 13c) 


The K’s are given by (7.4). 

These conditions are more general than any given previously. 

If there is strong stability (i.e., if the real parts of a, a, are negative and large), 
then the first terms in (7.1) die away quickly. In fact, the rapid oscillations are 


damped out, and we are left with 


g oy = | 
E = — 3 cos 6- (ac ‘ws + } ), 
| (7.14) 
§ ry 
n= —cos6-X’, 
With these we associate the last two equations of (6.6), viz. 
w = F3/m — g sin 8@, 
(7.15) 


W3 = G;/C, 


and also n= —(@ cos 6+76). 

In (7.14), (7.15) and the last of (6.5) we have seven real equations for the seven 
quantities u, v, w, 0, b, ws, 20. E is a function of w and w; as in (7.3); it also involves 20, 
since the properties of the air depend on 29 and aerodynamic functions X, Y, X’, Y’ 
depend on the properties of the air. The above equations determine the motion of the 
stable shell. 

We note that the equations (7.14), (7.15) are simply (6.6) with the terms &, 9 de- 
leted. To test whether this treatment is valid, we should solve (7.14), (7.15) for &, 7, 
calculate £, 9 by differentiating these solutions, and compare these calculated values 
with the other terms in (6.6). They should, of course, turn out to be small. 

8. Stability in the absence of Magnus effects. If we accept the linear law (4.11), 
the aerodynamic force system (4.13) is the most general possible. As we shall see in 
section 10, the force system of Fowler et al. is a special case. The system (4.13) con- 
tains ten real functions, and it appears impossible to make any deductions of physical 
interest without introducing some simplifications. We shall retain a force system a 
little more general than that of Fowler et al.; our system satisfies the fundamental 
condition of invariance with respect to shift of mass center, whereas theirs does not. 

Let us refer to (4.20), (4.21), and assume that all Magnus effects vanish, except G3; 


this means that 


Pt = Oy" = Pi* = Qi* = 0. (8.1) 
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This leaves us with four real aerodynamic functions, in addition to F; and G; 
P*<0, Q*>0, Pi*<0, Qi* <0. (8.2) 
There can be no doubt that these inequalities are physically valid. 
We now transform to the mass center O by (4.18). We find 
P» => OQ; => P{ = O2 = 0. (8.3) 
Thus the Magnus effects do not reappear with change of base-point; in fact, the 
vanishing of Magnus effects 1s an invariant condition. For base-point O there are again 
just four real aerodynamic functions in addition to F; and G;: 
P; _ ri } 


Oz = O* — rP;*, 


> (8.4) 
Py = Pj*+rPf, | 
Qi = Qi* + rPi* + r(— 0% + rP*). 
Then by (6.4), (7.4) and (8.3), 
K, = X,+ Yi = Pi/m+ QO1/A, ) 
Ks = Cas, 
wP3 1 | 
K; = wX{ + XiVi + X27¥.=- r + ; (PiQi + PQs), + (8.5) 
d mM. 
- - C'wsP 
K, = C’w3X, = —— - 
m 
The stability conditions (7.13) read 
Xi+ Yi $0, (8. 6a) 
(C’w3)? + 4wXf + (Xi + Vi)? + 4(XiVi + XY.) 2 0 (8. 6b) 
XiV{ (C’w3)? + (Xi, + Y{)*(wX7 +XiV{ + X27 V2) 2 O. (8. 6c) 


These are the stability conditions in the absence of Magnus effects. Now by (4.23), 
(6.4), (8.4), we have (since A* =mr?+A) 








: pa*w 
X, = —-—fi*, 
m 
pa*w r 
Y¥,= —(g*+— jt), 
m a | 
oe paw / , . | 
Xz; =-- - fi* + — fi 
y, a | 
A , (8.7) 
# patw - r P a aa" 
Vi = ———| gi* + — (gs + f2*) + — fi |, 
A a a? 
, r? pa'w / . 1%) A* * 
Xit+ Yi = ——| gi* + — (g.* + fz") + ~ fi ; 
A a ma* 
“f ys rs , p*aew* / / 
XiVi + X2V2= (fi*gi* _ go*fe*). 


mA } 
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If we substitute these expressions in (8.6) we get stability conditions in terms of the 
functions (f*, g*). However, these conditions are somewhat complicated, and we shall 
make approximations. 

The f’s of Fowler et al. hardly exceed 10 in value. Our (f*, g*) functions are defined 
in a slightly different way, but it certainly seems legitimate to assert that the dimen- 
sionless quantities 

pa* 
e=—f (8.8) 
m 
are much less than unity, f standing for any one of the (f*, g*) functions. Then it is 
clear that 
(X, a Y;/ XiV{ + X? Y> 
are both small relative to wX7. Consequently our stability conditions (8.6) may be 


simplified to 


X.+ Fi s@ (8.9a) 
(C’w3)? + 4wX/ = 0, (8. 9b) 
XiY/ (C'ws)? + (Xi + Y{)*wXs 2 0. (8. 9c) 


It will be noticed that Y2 has disappeared from the stability conditions in the last 
approximation. This aerodynamic function corresponds to cross force due to cross 
spin relative to the mass center [cf. (6.4) and (4.20) |. Thus it might be asserted that, 
for the discussion of stability in the absence of Magnus effects, cross force due to 
cross spin may be neglected. But this statement is not entirely correct, because this 
cross force contributes to the moment Y/,and Yj remains in the stability conditions. 

Let us examine the first stability condition (8.9a). On substitution from (8.7) it 
reads 

r F * 
— (g.” + ff") + gi? + — fr’ =O. (8. 10a) 
I ma’ 


( 


If r is positive (so that the mass center lies behind the centroid), this inequality is 
certainly satisfied; it is also satisfied for some negative range of r. But an interesting 
question arises: Can we make the shell unstable by pushing its mass center forward 
towards the nose? This is hardly to be expected on physical grounds, and it may well 
be that (8.10a) is satisfied for all permissible values of r, i.e., all values which place 
the mass center inside the shell. 

It is tedious (and perhaps of little physical interest) to discuss the other stability 
conditions for sufficiently large negative values of r. We shall therefore assume either 
that 7 is positive, or, if it is negative, it is such that (8.10a) is satisfied and also 


Xe <8, ry. <6. (8.11) 
Let us write 
(C’w3)? 


9 ae (8.12) 
— 4wX¢ 
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This is essentially the same as the usual stability factor.> Then the second stability 
condition (8.9b) takes the familiar form 


IV 


i” (8.13b) 


‘ 
while the third condition (8.9c) may be written 


(X¥,+ YY) ; 
s2— . (8. 13c) 
4X,V{ 
Since the fraction on the right is never less than unity, this condition replaces (8.13b). 
Let us substitute in (8.13c) from (8.7) and sum up as follows: 
STABILITY CONDITION. The following assumptions are made: 
(i) Magnus effects are negligible (except that G; may exist). 
(ii) The quantities € of (8.8) are very small. 
(iii) The mass center is behind the centroid, or, if in front, its negative coordinate 
r is such that (8.10a) is satisfied and also 


r 
f:* +—f*> 09, 


} (8. 14) 
gi* + — (go* + fe*) + fee fi*> 0. 
a a? 
Then the motion of the shell is stable if, and only if, 
> ma? [ei* + (r/a)(gs +") + (A*/ma?) fe) (8.15) 
— 44 fi*[ee* + (r/a)(gs* + fs") + (72/02) f*] 
where s is as in (8.12), or equivalently 
— : ea , 1 = A* — mr’. (8.16) 


~ Apar. Aw?|fs* + (r/a) ft 
To show the dependence on 7 more explicitly, we introduce the dimensionless quantity 


9 
i W3 


7. = : , (8.17) 
4pa*mw- 
so that 
a i r 
p= s— fi* + fir}. (8.18) 
ma? a 


Then the sole condition for stability reads 


Ge* + (r/a) fr ‘) Lgi* + (r/a)(go* + fe") + (A*/ma*) ft |? 


_ (8.19) 
4f*[ei* + (r/a)(gx* + fi*) + (r°/a%) fe] 


’T, J. Haves, Elements of ordnance, J. Wiley and Sons, New York, 1938, p. 418. 
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Since A* is the transverse moment of inertia at the centroid, the position of the mass 
center is involved in this formula only in the symbol r. 

We see therefore that the usually accepted criterion for stability (8.13b) is not 
the truce one; it must be replaced by one of the inequalities (8.13c), (8.15) or(8. 19), 
which are of course equivalent to one another. As we remarked in the Introduction, 
the existence of a second condition for stability has been noticed by R. H. Kent. We 
shall refer to stability again in section 10. 

9. The trajectory of a stable shell in the absence of Magnus effects. Let us as- 
sume, as in the preceding section, that Magnus effects are absent, except that G; 
may exist. Then, using (8.3) and (6.4) with (7.14), we get for the trajectory of a stable 
shell, after the disturbance has been damped out, 


g . 
—— 7 cos A( iC’ ws + Y/ }s 


' 


wtr 
| 


(9.1) 


n=1 ~ cos 6-X¢, n= — ¢cos0 — #6. | 
E ] 

Here F is as in (7.3); let us make the approximation indicated above (8.9), so that 
E = wX3 + iC’w3X}. (9.2) 


Splitting (9.1) into real and imaginary parts we get 


g ade | 
u=— — Cos 6[X4(C’ws)? + wX/} {] | 
E |? | 
(9.3) 
§ 7 rr 
—- +o 0-C’wswXe , 


4 |< 


(where we have dropped a term X,¥7 in comparison with wX¢) and 


(9.4) 


—— cos 0-w(X)?. 
| Fi? 


We shall assume, as in section 8, that X,, X?, Y/ are all negative. Further, since the 
shell is stable, we have as in (8.9c) 
X,V/ (C'os)? + (Xi + Vi)?*wXf 2 0. (9.5) 
But 
(X,.4+¥ipore, Xx<6 
and therefore 
X,V{ (C’w3)? + Y{2*wX!s 2 O. (9.6) 


It follows at once from (9.3) that u is positive. This means that the nose of the shell 


points above the trajectory. 
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From (9.4) we see that ¢<0 if w;>0. Thus for positive (right-handed) spin the 
vertical plane through, the axis of the shell turns to the right.6 For negative spin it 
turns to the left. 

These two facts are well known to be true in practice. 

There remain two outstanding physical facts to explain. These are (i) the trailing 
of the shell along the trajectory, (ii) the drift. 

We see from (9.4) that 6 is negative, i.e., the inclination of the axis of the shell 
to the horizontal decreases steadily. But does it decrease at that rate required for 
trailing? We must be careful to avoid a circular argument. We have assumed that 
trailing takes place—otherwise the yaw is not small, and all our arguments are based 
on the smallness of the yaw. We must now verify that 6, as given by (9.4), is approxi- 
mately equal to the rate of turning of the tangent to the trajectory of the mass center. 
The theory of the plane particle-trajectory gives, on resolution along the normal, 

g COS Ao 


66 = p] (9. 7) 


w 
where 4 is the inclination of the tangent to the horizontal. To establish the required 
result, we must compare this with (9.4), and show that 

| E 
(wx)? 


approximately. Now by (9.2), (8.12), (8.7), this fraction is 


af Clos \! 4sX1° 
1 + X\ ) = 1 eet 


ry 
WAos wx 2 


9 








= 1, (9.8) 














3 *2 
wit Le sail (9.9) 
m ma? f*’ + (r/a)fi* 
The last expression here is of the order of se, where € is as in (8.8). Hence, unless the 
stability factor s is very great, this expression is very small, and the condition of 
trailing is approximately fulfilled. 

It is interesting that if s is very great the verification breaks down, for this is just 
what we would expect. If, by some mechanism, an enormous spin were imparted to a 
shell, the gyroscopic stability would be so great that the direction of the axis would 
remain fixed and the shell would not trail. 

To discuss the drift, we write down (6.5) again: 


Xo + ivo = (usin 6 + iv + w cos A)e**. (9.10) 
This is the horizontal velocity of the mass center in complex form. Consider the com- 
plex quantity 
£o + t¥o 
— (9.11) 
Xo + 190 





a+ iB= 


It is obvious that the vector %9+7 turns to the left if 8B is positive, and to the right 
if 8 is negative. It is our business to investigate the sign of B. 


6 Hayes, op. cit., 420. 





j 
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We differentiate (9.10) logarithmically and simplify the result by the fact that 


u/w and v/w are small. This gives 


d v e 
meets : Q 9.12 
B As cos ;) t¢ ( ) 


With the approximation (9.8), we have from (9.3), (9.4) 





v C'w ‘ Cw 
aoe ae es ote ’ q=-— lin Xi, (9.13) 
w cos 6 wX7 wXg 
and so 
B d . 
— = — log |Z| + Xi, (9.14) 
z dt 
where 
C’w; 
Ee oi (9.15) 
w'X9 


As a terminology, let us say that a shell has standard drift when it goes to the 
right (8 <0) for right-handed spin (ws;>0), and vice versa. Now Z has the same sign 
as w;. Hence we get a standard drift if 





d 
— = — log |Z|+ X1 < 0. (9.16) 
Z dt 
Substituting from (8.7), we see that this condition for standard drift reads 
d Ww e 
(9.17) 





ji* > ~ tog | 
; pa’w dt pw*(fe* + (r/a)fi*) 


Let us look into the meaning of this inequality, assuming that the dimensionless 
aerodynamic functions are constants. This corresponds to a subsonic velocity [cf. 
(4.25) |. Further, let the trajectory be flat, so that p is constant and @ so small that 


it may be neglected. 


Then by (6.6) and (4.24) 


. pa* ; pa‘ 
SS re og = Cc W383. (9.18) 


m 
Let s be the are length of the trajectory (do not confuse with the stability factor). 
Then w=ds/dt, w=wdw/ds, and so we have 
1 dw a? 1 dws pa‘ 
— : (9.19) 


iendieecad a ie ee ere 
w ds m W3 ds é 


The right-hand side of (9.17) becomes 


m a : W3 34 ma? 
ote se Se ae a er en 
pa? ds w* ( 
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No 
nN 
to 


and so the condition for standard drift reads 


ma” : 
f;* + C £3 > 3/3. (9.20) 


We may write f; in place of f*—they are equal. We note that ma?/C will lie between 
1 and 2. 

We observe from (9.20) that cross wind force and axial torque tend to give stand- 
ard drift, but axial force acts the other way. Let us use the numerical values of Fowler 
et al. in (9.20). We have ([1], pp. 306, 309) 


fi = fut fre = fy = 3.34, 
fs _ Tr = 0.34. 


(9.21) 
We see that (9.20) is liberally satisfied, even if g;s=0. Thus the present theory appears 
adequate to explain drift without bringing in Magnus effects. 

10. The aerodynamic force system of Fowler, Gallop, Lock, and Richmond.' In 
the preparation of this section we are very much indebted to Professor E. J. McShane, 
who read our paper in its original form and pointed out in detail the connections be- 
tween our work and that of Fowler et al. 


— 
a 


Velocity of 0 











Fic, 3 


The axis of the shell is indicated in Figure 3; Ois the mass center and 6 the yaw. 
The aerodynamic force system of Fowler et al. is represented by seven vectors—three 
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forces (plain arrows) and four couples (arrows with crossbars). Their terminology is 


as follows: 


= drag, 
= cross wind force, | 
| 
} 


= swerving force, 





moment tending to increase yaw, (10.1) 
= yawing moment due to yawing, 


f 
. } 
= axial moment, 
} 


ll 


= swerving moment. 


We shall use the notation of the present paper for velocity, angular velocity and the 
radius of cross section of the shell (@), and consider only the case of small yaw 
(6=|£)|/w). Then the dimensionless aerodynamic functions of Fowler et al. are de- 


fined by 


R = pa’w*fr, } 
L = pa*w*bf, = pa’w| —| f1, | 
K = pa*ww35fx = pa*ws| &| fx, | 
M = pa*w*6fy = pa*w| é\ fur, (10.2) 
H= pat | n | tu, | 
IT = pa‘wws3fr, 


J = pa‘ww3bf7 = patw3\ &| fo. } 


Let i’, j’, k be an orthogonal triad of unit vectors, with k along the axis of the shell. 
The vector i’ lies as shown in the plane containing k and the velocity of O. Then, to 
the first order in 6, 


= Roi’ — Rk, | 3 


= — Li, 


(10. 3) 


~ £SAe ww 
\ 
Kaa 


I = — Jk, 
J = IT, 
where w/, w¢ are the components of w along i’, j’. 


Let i, j be any orthogonal unit vectors, perpendicular to k, so that the triad i, j, k 
is that considered in the present paper. It does not matter at present whether i lies 
in the vertical plane through k. We have 

ui + vj . . -“8ee 


ea eee 
s 1s 


wf i! + ws jf! = wri + wa. (10.4) 


The total aerodynamic force is 
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nN 


u u U 





’ v v _ u 
+i(- Ro ~ Lye} £~- 
| é| | é| | é| 
— ER, (10.5) 
and so, since wb = | | j 
R L iK ) 
P= y+ ifs = e(- —7—-+7,), | 
w |é | ¢| | (10.6) 


The total aerodynamic couple is 


v w4 u 
G=M+H+1+J=i( M -H 7". - 7") 


f | | n| |g 
u We v 
+i(-w pe - as ) 
| é| | | | é| 
— k/, (10.7) 
and so 
, ” g 7 E 
G= Gy, + 1G2 = iM 7 “lies H ee ee J — mo 
| | || le] | (10.8) 
G3; = — f, J 


Certain quantities are defined as follows: 


L K Mw 
k= ; Xr =— an - B= a ; | 
m)\ &} mw; | & E 
re. > (10.9) 
R H I Jw | 
yv=xKk+—) h=— r=—, ¥=—7° | 
mw A|\n Cus Cus | } 
Then (10.6), (10.8) give 
F = £(idmw; — my), F; = — mw(v — x), ) 
Bn Cas; (10. 10) 
G= :(- i—— =") — Ah, G3; = — Ca,I". 


Comparing these with (4.15), we see that the force system of Fowler et al. is a par- 


ticular case of our general system, with 


P = — mv + imdws, ) 
yCw np } 10. 11) 
a ee ae ( 
w w | 
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The general system has eight real parts in these terms; the system of Fowler etal. 


has only five: 


P, = — my, Ps = MAW, 
yCws » | 
P{ ee ia’ Racca P? =_——; (10. 12) 
Ww w 
OY => , —"A h. 


It is clear from the transformation (4.16) that Q=0 is not invariant with respect to 
shift of mass center. Thus Eqs. (10.11) describing the aerodynamic force system can- 
not be valid in general. It may happen of course that they are true for one particular 
mass center, but they cannot remain true when we shift the mass center. 

Fowler et al. find little evidence for the existence of the Magnus effects J, K, or 
equivalently y, \. If we put them equal to zero, the survivors in (10.12) are 


P, = — mv, PJ =-—») Qi = —Ah. (10. 13) 


These should be compared with (8.4), which are the general survivors in the absence 
of Magnus effects. We note that Q2 is absent from (10.13), which means that the mass 
center is chosen so that Qs*—rP;}* is zero, or at least negligible. 

3, (6.4) we obtain from (10.13) 


XxX, = — DV, Xs = | ’ Yi ==-_— h, (10.14) 
and so the stability condition (8.13c) reads 
(vy + h)? 
:< . (10.15) 
duh 


We have referred in the Introduction to a second stability condition implicit in the_ 
work of Fowler et al.; it is 
(x + h)? 


-—- (10. 16) 
4kh 


$2 


The difference between (10.15) and (10.16) does not appear to be very great in prac- 
tice. It is a question of replacing v by x, and by (10.9), (10.2) 


Rae! R \¢} fr 1 
maniseh i Seas Sn ae hee heey (10. 17) 


roughly. 

There are very simple relationships between the dimensionless aerodynamic func- 
tions in the two theories. We take the mass center O as base-point, and use (4.23) 
without asterisks, together with (10.12), (10.9), (10.2); we find 

fi =fvntfi = fn, fo = fx, ) 
- F ‘i s t (10.18) 
ji =fa, fz = fu, gi = fu.) 
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The functions gi, ge, g? are zero in the theory of Fowler et al. 

As the paper of Fowler et al. is one of the basic papers of modern ballistics, it will 
be useful to summarize our criticisms as follows: 

(i) Their aerodynamic force system is not the most general system consistent 
with 

(a) the aerodynamic hypothesis, 
(b) linear dependence on the cross components in the case of small yaw, 
(c) the symmetry of the shell. 

(ii) Their system does not satisfy the fundamental requirement of invariance with 
respect to shift of mass center. 

(iii) If only shells with mass centers near their centroids are considered, it may be 
that the above theoretical objections are of small practical importance. 

We believe that our exact dynamical equations (6.3) provide a clearer approach 
to the problem of the spinning shell than do the dynamical equations of Fowler et al. 
But it is frankly admitted that our simple treatment of the equations of motion in 
section 7 does not appear to be as satisfactory mathematically as their method. We 
have made the plausible but rather crude assumption that it is permissible to regard 
cos 6, w, wz; as constant during the oscillation. It would be interesting to apply their 
more refined methods to our differential equations, but this we must defer for the 


present. 
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TRANSFORMATION GROUPS OF THE 
THERMODYNAMIC VARIABLES* 


BY 


WALLACE D. HAYES 
Lockheed Aircraft Corporation 


Abstract. A certain class of transformations on the thermodynamic variables E, H, F, G, S, T, P, and 
V which leave the fundamental equations invariant is investigated and found to form a group of order 
thirty-two. The quotient group with respect to a normal subgroup of order four gives the octic group ob- 
tained by other investigators, the normal subgroup containing trivial but non-excludable transformations. 
In contradistinction to previous investigators, it is not necessary to use absolute values or a rule of signs. 
Examples are given of the application of the transformations. 


Certain transformations on the fundamental thermodynaniic variables will change 
members of a large class of thermodynamic equations valid for reversible processes 
into other valid equations of similar form. These transformations have been investi- 
gated by Koenig! and Buckley? and found to form the group of order eight called the 
octic group. Koenig restricted his transformations to pure substitutions, or permuta- 
tions, took care of a difficulty in sign by introducing absolute values and a rule of 
signs, and discussed a geometric method of exhibiting the transformation group. 
Buckley showed that Koenig’s group could be derived in part by Lie’s theory of con- 
tact transformations, and listed a number of families of thermodynamic equations 
to which Koenig’s transformations apply. Although of course mathematically correct, 
the application of Lie’s theory is not essential in this case. 

In order to eliminate the inconvenient and somewhat disturbing use of absolute 
values and a rule of signs, the following exposition of the theory of these transforma- 
tions is presented. The transformations considered are not limited to pure permuta- 
tions but allow changes in sign, and the octic group is finally obtained out of a larger 
transformation group as a quotient group without the necessity of using absolute 
values or a rule of signs. The transformations are represented by matrices whose ele- 
ments in any single row or column are all null except for one element which equals 
1 or —1. 

The thermodynamic quantities involved are*: the internal energy E, the en- 
thalpy 77, Helmholtz’ function F, Gibbs’ function G, the entropy S, the absolute 
temperature 7 (intensive), the absolute pressure P (intensive), and the volume V. 
All these quantities except J and P are extensive quantities. The quantities H, F, 
and G are defined relative to E by 

* Received Aug. 26, 1945. 

1 F, O. Koenig, Families of thermodynamic equations, I, Journal of Chemical Physics, 3, 29 (1935). 

2 F. Buckley, Transformations of the fundamental equations of thermodynamics, Journal of Research 
NBS, 33, 213 (1944). 

® This notation is perhaps the most common in scientific literature, with the symbol U often used in 
place of Z. The functions F and G, invented by Helmholtz and Gibbs respectively, are generally known as 
the free energy and the thermodynamic potential, while the American Standards Association has used 
the term free enthalpy for the function G. The standard usage of American chemists is that of Lewis and 
Randall, where the functions F and G are denoted by A and F, respectively, and are termed the work func- 
tion and the free energy. 
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H = E+ PV, (1a) 
F= E— ST, (1b) 
G= E+ PV — ST. (1c) 


These definitions, together with the fundamental thermodynamic equation for dE 
give the equations 


dE = TdS — PdV, (2a) 
dH = TdS + VdP, (2b) 
dF = — SdT — Pdv, (2c) 
dG = — SdT + VdP. (2d) 


The transformations considered are all transformations which leave Eqs. (1) and 
(2) invariant, such transformations preserving the validity of any equations derived 
from Eqs. (1) and (2). The class of equations to which the transformations apply is 
therefore the class of equations thus derived. If the symbol x is used to denote un- 
determined matrix elements, the transformations will be of the form 








eo = a [* x x id r° 0 0 * "7 eo" 
H 22 € = 0 0 0 0 il 
F’ ££ *% 2 0 0 0 0 F 
G’ Lx 2 « #5 0 0 OC ‘@) G 
= (3) 
a ' ). & -@ sf [2s s £ 7 § 
z 6-0 8 @ 4 2-8 * T 
Pp’ 0 6 0 O xX xX 4 P 
Rie [ 0 0 0. @F bs = & <3 | iV J 


























From the invariance of Eqs. (2), the following seven observations on the trans- 
formations are made: 

I. The off-diagonal 4 by 4 submatrices are necessarily null. This fact allows the 
transformations to be put in the separated form 


Mm be x x a1 277 
pi x x x x H 
= 5 (4a) 
F’ x x x x F 
ra Lx x X Ri iG ; 
ba il x ‘] l* ‘]] he i 
+ jas x x 2 r 
= . ( 4b ) 




















ii | x ‘] | x ‘] P 
ae tLe & x 2 | bee | 
II. If one of the variables (ZJ/ FG) is changed in sign, all of them must be thus 


changed. Such transformations may at this point be excluded as trivial. This exclusion 
limits the transformation (4a) to pure permutations. 
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III. The invariance of the equation 
E-H-F+G=0 (5) 


derivable from (1) may be used to limit the transformations (4a) to eight in number, 
all of which are of the type considered. 

IV. Since S is always associated with T in Eqs. (1) and (2), as is P with V, the 
form of the transformations (4b) must be as shown with two diagonally opposed 2 
by 2 submatrices null. 

V. The 2 by 2 submatrices of Eqs. (4b) are necessarily of one of five forms, which 


are abbreviated thus: 


2 ek: 








= ¢, (6a) 
| Oo 
pe 07 
= = ¢, (6b) 
| © =) 
2 
= §, (6c) 
Pe ‘in 
r @ =e (6d) 
== £, 
| 1 O 
r 0 00 ; 
= 0. (6e) 
eae 


VI. Transformations of the type in which, for example, both P and V are changed 
in sign are admittedly trivial but cannot be excluded because they are necessary for 
closure of the group of transformations (4b). 

VII. A given transformation (4b) defines at most one transformation (4a). The 
converse is not true, however, and the correspondence is found to be four to one. 
Thus the number of transformations (4b) is thirty-two. 

The eight transformations (4a), the corresponding thirty-two transformations 
(4b) expressed using the abbreviations in Eqs. (6), together with eight symbols repre- 
senting group elements, are listed in Table I. 

The transformations (4a) form a group of order eight which is designated as M. 
The transformations (4b) form a group of order thirty-two which is designated as G. 
The group G is four to one homomorphic to M, the normal subgroup 


ie ae ee es eae fm 


of G corresponding to the indentity element of M. From this correspondence is estab- 
lished the isomorphism 


M ~G/N (8a) 


or the congruence 


M =Gmod N. (8b) 
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This group is the octic group, whose multiplication table in terms of the group ele- 
ments shown in Table I is given in Table II. This multiplication table is consistent 


with matrix multiplication of the matrices representing either transformations (4a) 


TABLE I: The transformations of the thermodynamic variables. 


| 


— ———— — = —=s pn nn ——— — = 7 ———— 
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or (4b). The non-identical transformations of the subgroup N are those of the trivial 
type mentioned in observation VI, and their elimination in the process giving Eqs. (8) 


is tantamount to disregarding a change in sign of both S and T or of both P and V. 
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To illustrate the general application of these transformations the equations ob- 


tained by transforming two given thermodynamic equations are shown in Table III. 
The two given equations are those shown in the table opposite the symbol mm, which 


TABLE IT: 

my, 
m my 
Ms ms 
m m3 
ms ms, 
m Ms 
me me 
m7 m7 
m Ms 


TABLE III: 


rransformation 


( as ) 
m - 
oT 
ae) 
mM 
os 
(— ) 
m7 
OV /s 
(= ) 
Ms ae 


ms m3; 
Me ms 
mM ms, 
m4 my 
m3; m 

me Mm; 
Ms ms 
ms ms 
m7 Me 


ms 


ms 
m3 
me 


ms 
m7 
m6 
ms 


Multiplication table for the octic group. 





ms 
Me 


m7 
my 
Me 
ms, 
m3 


me my ms 
ms my ms 
ms; ms m7 
m7 mes ms 
ms ms me 
m2 m3 ms 
Mm, ms, m3 
ms me my; 
ms, my Me 





Examples of equations obtained from the transformations. 





First seein 


iy 
\ 7 oT 


r), 





1 aG 
I ead s(Gs), 
1 /aG ° 
7 * V (=). 
1 /dE \ 
s an f (). 
“1 (0H 
P % i G ), 
1 (OF \ 
ees! 
+1 /a#HY 
= é (F), 


Second mene 








_ Vv - ) 
‘<i Rj S\@P Js 








§ 
S17G\y 
ss V Ar ), 
a me 
x P a P = . 
“VO o«1faH\ 
- T i T aP a 
— ve 
ge es sr yy 
ft _ 
F3(2), 
— en 
ic ea P (or ), 


represents the identity transformation. As an example of the carrying out of one of 
these transformations, the m7; transformation of the second equation of Table III is 


here given in detail. 


0 


IT. For the ( 


to be 


used, 


8 @ E 
0 0 1 H 
ooollFl, 
0 1 0d Le 


The transformed value of E is shown by 


H 


(9a) 


STPV) transformation any one of the four matrices given may be 
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i) 
w 
to 


Hence the equation 


aP P 1 /dE 
(22) - £442) 
OT Jy 7 (a OV T 
(<. ) T 1 (=) 
wis VV \@lc 


Since m;=msms; from Table II, Eq. (10b) can also be obtained by 
to the equation obtained by the m; transformation of Eq. (10a). 


is transformed into 


(9b) 


(10a) 


(10b) 


applying ms; 








CONSTRUCTION OF A COMPLETE SET OF SOLUTIONS 
OF A LINEAR PARTIAL DIFFERENTIAL EQUATION IN 
TWO VARIABLES, BY USE OF PUNCH 
CARD MACHINES* 


BY 
STEFAN BERGMAN** 


Brown Unwersity 


1. Introduction. Many problems in Engineering and Physics lead to the deter- 
mination of a function U satisfying the partial differential equation 


vse 1 (= 4 “) 1 ou 4 1 8 au oa y , 2 
AU) = - —— } + — a(x, y) — — 2, ¥) — x,vU =0 : 
; > ol x, 4 rm ; (x, ay Y ; (1.1) 


Ox 0 y? 


“ 


in a domain D, which function assumes certain prescribed values on the boundary of 
D. Here a, 8, y are polynomials in x, y (or can be approximated by such polynomials). 

Recently a procedure, “the method of particular solutions,” has been developed 
for the solution of problems of this kind [2].¢ The idea of this method is to determine 
at first a “complete” set of particular solutions of (1.1), i.e., a set of functions each of 
which satisfies (1.1); this set [denoted by p,(x, y), (v=1, 2,---+)], is chosen in such 
manner as to possess the property that every solution U can be approximated in any sim- 
ply connected domain by a conveniently chosen finite combination 

N 


> ap,(x, y) (1.2) 

v=l 
of the above particular solutions. In the case of the Laplace equation (0?U/dx*) 
+(0°U/dy*?) =0, such a set can be obtained by taking the real and imaginary parts 
of the powers (x+iy)’"!, (v=1, 2, -- + );i.e., prp=1, po=x, ps=y¥, Pa=x?—y*, ps =2xy, 
etc. In the following, a procedure will be described for finding analogous solutions 
for any equation of the form (1.1). The second step of the method consists in indicat- 
ing a rule for determining the a, so that (1.2) assumes on the boundary of D values 
which approximate the prescribed values of U. In order to apply the method of par- 
ticular solutions to an actual problem and obtain numerical results of interest to an 
engineer or physicist, it is frequently necessary to carry out lengthy computations. 
These computations can as a rule be performed most efficiently by the use of special 
computing devices, such as punch card machines. Before the computations can be 
carried out on such machines, however, it is necessary to organize the computations 
so that they can be given to the operators of the machines. This organization is often 
a problem in itself, as in the case of the example given below. 

In the present paper there is described a working procedure for computing the 

above set of particular solutions, and for carrying out the associated numerical com- 
putations by the use of punch card machines. 


* Received Oct. 11, 1945. 
** Now at Harvard University. 
+ Numbers in brackets refer to references at the end of this paper. 
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2. Notation. In the following, it is convenient to use the complex variables z and 
= instead of x and y. We have s=x-+iy, 3=x—iy, where x and y are cartesian coordi- 


nates in the plane. We then have 


au 1 fal dU rey a 1 /oU dU 
Oz 2 Ox Ov OZ 2 Ox ov 


l 1 f/oru ri ad Of 
U.z= AU = ( + ). 
4 4 Ox Oy? 


In terms of z and 3 Eq. (1.1) then assumes the form 
U,stAU,+ BU; +CU = 0, t2.3) 


where 2A =a+tif, 2B =a—if8, C=y7. 

REMARK. In the above case, B =A. However, in the following we derive formulas 
without making this assumption and we note that in general (2.1) is equivalent to a 
system of two equations, one corresponding to the real part and one to the imaginary 
part. 

3. Some previous results on particular solutions. The method to be described in 
the following is based on the following theorem, proved in [2] p. 542. 


THEOREM. For every partial differential equation (2.1) there exists a function 


E(z, 2, t) = exp| - f ade || + >. 2rgrQ'”)(z, 2)| (3.1) 
0 1 


Pp 


such that the functions* 


l 
Povr41(3, 3) = Re l* f E(z, 2, (1 — #)"" at] | 
= F  @=OL A), G.2) 
| 
| 


el 
Poi2(, 2) = Im l= f E(z, 2, (1 — #)° vat] 


1 


form the required complete set of particular solutions,** 1.e., a set possessing the properties 
mentioned in §1. The O°” aret given by the recurrence formula 


* As has been proved in [1], the series (3.1) converges for |x| <<, |y| <, |t] $1 and therefore 
the functions p, are entire functions. We note that they are independent of the domain for which the 
boundary value problem is considered. 

In the N.A.C.A: Technical Notes Nos. 972, 973, 1018 and 1096 the functions Q(z, 2) have been com- 
puted for the compressibility equations. ; 

** If E(z, Z,¢) is real, then p2(z, 2) vanishes identically and we have to change the numeration of the 
p(s, 2) accordingly. The proof that (7.1) where f is an arbitrary function of a complex variable is a solu- 
tion of (1.1) is given in [1] §1. Substituting Q instead of ¢@») /s? into (1.16) of [1] p. 1174 and integrat- 
ing with respect to Z we obtain (3.3). (Note that there are some misprints in [1] which are indicated in 
the Trans. of the Amer. Math. Soc. vol. 57, p. 311, Footnote 15.) 

+ Unless otherwise indicated the arguments of the functions Q, F, A, etc., are (z, 2). 
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1 . | 

) - 2f Fdz, 

0 | 

) - ' 
p < »—1) (p—1) ; (p—1) (p—1 + 5.3) 

QO =- g —Q,  (, 0) +f (D0; + FQ e al | 

2p esi 1 a | 

(¢ = 2, 3,4,---). 
0" (s 0) = 0 (p = 1, 2, 3, ), 
where 
D = -{ 1.d= + B, F=-AB-—A,+( (3.4) 
0 


If the coefficients A, B, C of (2.1) are polynomials 


\ \ Vv N N N 
A => > amet", B= > > banat", C= DoD Can8E", (3.5) 
m==0 n=-( m=0 n=0 m=-0 n—0 
then for the coefficients d,,, and fm,, of the developments 


ueseM 


M Vv 
D ts Zz p dn x32", F - 7 > | ie LD (3.6) 


m=-0 n=0 m=0 n=0 


we obtain the relations 


(me + 1)Gm+1,n—1 


do.n = bon; : = + Cncte (3.7) 
n 

Juv = | (m + 1) @m41,n + > >» amines] + Cm,n- (3.8) 
i=0 j=0 


4. Determination of the Q‘”). Since the A, B, C (see (2.1)) are assumed to be 
polynomials, D and F are also polynomials in 2, 2, as is indicated in (3.6). (Both F 
and D are assumed to be of the same degree; this is always permissible, since some 


f's or d's can be assumed to equal zero.) We now write* 


(p—1) (p—1) i_j (p—1) 
O° =Vdas 22, go =9, (4.1) 
’ j 


aa 


and proceed to the representation of g@), in terms of the gj¥~”. Formal computations 


yield 


(p—1) #7 


of = DL alee = DDG t pala e's, 
i j i j 


po . (p- . , . 
* The relations q?, ” =0 follow from the equation Q‘°-))(s, 0) =0. See (3.3). Note that if one of the 
subscripts m, 2 of d»,» becomes negative or larger than J, it is necessary to substitute 0 for the corre- 


sponding dm.n. The same holds for fnn- 
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rg” = | Lee (fro + ford + fusd +> || DY 8 (giak + 48 +--- | 
=) z( pp et ae a ere’ (4.2) 


i=] j=] 


ll 


t 


po.’ = | Y 2(dr0 + dyad + dye + - | DY 2 (gia + 296.28 + 3qi.58 +--- | 


= > z( > XG + Ddn-in aiina )e°8 


s=0 j7=—0 


m 


n+l 
3 x( > oh ar) a 


i=? j=l 


>» (p--1 m = 
Fo’ + DQ: a Zz, 2 ‘E 14m 1,09i,1 + [2d » i,09i,2 = (dm 1,1 + Su i,0) gird | 


1=0 


+ 13800008 + leat fa-iadia t+ Geaat faeidealt +-*:3 
= ~ . e (p—1) <e (p—1) mon 
= z. } | a he (Jd m-i n—j+1 + J m—i,n—5) Qi, j + (n + 1)>> Gu 09: = | SS a (4,3) 
m n i=0 j=l i=0 
Therefore 
y f (4.4 
m i eS m,k—ly . ) 
q BR: 1 : 
(p) 2 = (p—1) (p—1) 
ae er ee ; dm—i,0Gi 1 ++ (m +> 1)Qm+1,1 ’ (p = Z. 3. eo 
2p — 1L ‘Ho 
(p) 2 1 m n—-l (p—1) 
°° OI = > (Rigs + ; ae ~1) i,k | 
2p —1L Mm iso kai (4.5) 
| 


(p-1 


= ) (p—1) 
5 :» (dim 1,09: n -+- (m + Danan) | 
i=0 
(m = 0,1,2,---,n = 2,3,4,---,p =2,3,---).]J 


5. Evaluation of the g”,. If we assume that the d,,,, and fm,, are determined, the 
evaluation of the qi), represents a considerable computational task. Merely writ- 
ing down the various pairs which have to be multiplied together and adding the re- 
sulting products is quite laborious. Therefore, it is convenient to organize the 
determination of the g®, in such manner that one can see automatically what 
quantities have to be multiplied and added. After such a procedure has been de- 
veloped it can be conveniently employed for performing the computation on punch 
card machines. 

We shall describe such a procedure, assuming for the sake of simplicity that M 
in (4.1) is a given number, namely M=3. The changes that have to be made when M 
assumes other values are in many instances clear. If however this is not the case, 
then we shall write these formulae out explicitly for the general case, i.e., retaining M. 


The expressions g?), (p >1) to be evaluated are given in (4.5) in the form of sums of 
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products of two sets of numbers; the first set are the coefficients of qi, , and the 
second set are gi". This suggests the construction of a set of cards and stencils. 
The coefficients of gi’? are entered on the cards. The stencils are cards with holes cut 
in them, and the quantities gi," are entered on the stencils in such a way that when 
a stencil is placed over a card properly, quantities on the cards appear in the holes 
in the stencil beside the gf") by which they are to be multiplied. 

A card is constructed for each value of n; these cards are called the [S, 2] cards. 
A stencil is constructed for each value of p; these stencils are called the [Q, p] 
stencils. In’order to compute g®, and hence construct the [Q, p] stencil, we require 
the [Q, —1] stencil and the [S, | cards corresponding to* n=1, 2, ---, p(M-+1). 

Each [S, »] card contains a matrix {s} of M@+2 rows and M+1 columns; 


ry indicates the row and k the column. Also when M =3, 


1 
5; = [(n 4. = 5)d, i b—r + | are | 


n 
n » (n) (nm) (n) 
11> fa.s/m i= fo.g/m Si.3 Sis/n Sia= fo.s/n 
21> [(n—3)ds.st+fs.0]/n 52.2= [(n—3)des+fe2]/m $23= [(n—3)dist+fie]/n Sa= [(n—3)do.s+fo.2|/n 


31= [(n —2)ds.o+fs.1 } n s32= [(m —2)d2,o+fe i] n Si3= [(n—2)diethfia |/n ssa= [(n—2)do.2+fo.1 | n 
“i= [(n—1)d3,1 +fs.o] Nn a= [(n —1)deit+fe0!/n “= [(n— 1 \dirthiol/n se [(m—1)do,:+fo,0!/n 


n) (n) (n) (n) 
$, = d;.9 si dso 55,3= dio s4= doo 


Fic. 1. A [S, n] card** in the case M=3. 


For n= M-+1 (i.e., in our case for » $4), the elements of the first M+2—n lines 
are to be set equal to zero, i.e., in our case sf =0 for kS5—n. 

A [Q, p] stencil consists of (6+1)(M+1)—p columns and (p+1)(M-+1) rows. 
The gi), in the last M columns and M-+1 last lines are equal to zero. In the case 
M =3, typical stencils are as shown in Figs. 2 and 3 (toe interiors of the rectangles 
should be cut out). 

The number of remaining columns in any stencil will be three greater and the 
number of rows four greater than in the stencil for the preceding value of p. The 
[Q, 1] stencil has 7 columns, 7 “deleted” columns, and 8 rows; hence the [Q, 2] 
stencil will have 10 columns, 10 “deleted” columns, and 12 rows; etc. (The [Q, 1] 
stencil has 7 columns (+7 “deleted” columns) and 8 rows and the S cards have 
4 columns and 5 rows.) See Fig. 2. 

Each member q®), on the [Q, p] stencils is specified by three numbers, one super- 
script and two subscripts. The superscript remains constant for each Q stencil, i.e., 
all numbers computed for the [Q, 2] stencil will have the form q®),, for the [Q, 3] 
stencil g®,, etc. The subscripts give the position of the number on the paper, indicat- 
ing column and row respectively; g{} means that the number is on the [Q, 1] stencil 
in the fourth column and third row; g®} that the number is on the [Q, 2] stencil in 


the fifth column and second row; etc. 


* We note that in increasing p to p+1 the number of [S, n| cards increases to (p+1)(M+1); 
however p(M-+1) of these cards are exactly the same cards used at the pth stage, so that it is necessary 
to prepare only M+1 new [S, »] cards. 

** The numbers on the S stencil are the coefficients of . gs in (4.5). 
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Fic. 2. Stencil [Q, 1]. 


Using (4.4), we determine ght’,, (m=0, 1, 2, 3, p=1, 2, 3, 4); gv, (m=4, 5, 6 or 
p=5, 6, 7, 8) we set equal to zero. Thus the stencil [Q, 1] has the appearance indi- 
cated in Fig. 2. The stencil [Q, 2] consists of 10 columns and 12 rows (besides the 
deleted columns, i.e. the holes in the stencils). See Fig. 3. We now proceed to the com- 
putation of the [Q, 2] stencil. To compute g®, we proceed as follows: the second sub- 
script (m, n) indicates the g®), on which we must fix our attention. Since there are 
four more rows and three more columns on the [Q, 2] stencil than on the [Q, 1] 
stencil, there will be no corresponding e. for these last columns and rows—we fill 
in these last 3 columns and last four rows on the [Q, 2] stencil with zeros. To con- 
struct the remainder of the [Q, 2] stencil, we place the [Q, 1] stencil on top of the 
(S, n] card so that the number s{?,.4,, on the [S, m] card occupies the space to the 
left of g®,.* There are then 20 numbers on the S card adjacent (on the left) to 20 
numbers on the [Q, 1] stencil, which indicates the products (20 in all) which have to 
be obtained; once the products are computed they are to be summed. 

In this manner we may tentatively fill in the whole [Q, 2] stencil. Now we com- 
pute the product of (m+1)q™. and add this to each of the “tentative ge.” If this 
number is multiplied by —2/3 [in the case of an arbitrary p by —2/(2p—1)], we ob- 
tain the “final g®),” of the stencil [Q, 2]. 

Having completed the [Q, 2] stencil we may repeat the above operations but with 
[Q, 2] and qe, replacing [Q, 1] and . respectively, to compute the [Q, 3] stencil. 
Similarly, we can compute as many Q stencils as desired. 


* See Fig. 4, where the [Q, 1] stencil covers the [s, 4] card, and the arrangement for the computa- 
° (2). « ° 
tion of q,, is indicated. 
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a a) a ee) 
sf sf sf sf 
5 va e 5 Soe 
re ee 


Fic. 4. [S, 5] card. 


6. Example. In this section, we shall illustrate our general descritpion by a spe- 
cific example. We shall indicate the operations to be performed on punch card ma- 
chines only. The arrangements and the methods of how the computation is to be per- 
formed on these machines can be found in books on punch card methods e.g. in[3]. 
More specifically we refer for operations concerning complex numbers to [4], and for 
computations used in the present paper to [5]. 











5) 1 
(5) |.(2)/ (5) q)) 3!) , (2) a! (2) ko ae 0) 
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wo!) [aq (2)} 65) Jq(2)) 9 £5) [C2] 5 (5) |g (2) jy (2) jg (2) (1) 
31 | 0,3 32 | 1,3] 33 | 2,3] 34 | 3,3 4,3 5,3 6,3 
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. = ’ . . e . (2 (6 
Fic. 5. When stencil [Q, 1] is placed on the [S, 5] card in order to compute q,.,, the s,., appear 


Phe » ©: ‘ ‘ (2) yy 
as indicated above. + WA CE neds gives the “tentative q, 5. 


As an illustration of the above method, we now give dmin, Om.n, Cm.» the specific 


values, 
2mn + 1t(m? — n?*) 
aoo = 1, Gun = ee Pee 
(m* + n?) 
1 
bo.o = 1, bo» = , 


(m? + n?) 
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' (m + n) + i(m — n) 
Coa = ‘ Caen ® - _ _ —— ° 

(m? + n*) 


In Table 1 the values of d,., and fm,, which have been determined according to 
(3.7) and (3.8) are tabulated for m=0, 1, 2, 3, 4, »=0, 1, 2, 3, 4. 


TABLE 1. The values of d,,,, and fn,,. 











m n Re(dy»,n) Im (dn ,n) Re(fn,n) | Im(fm,n) 
0 0 1.0000 0.0000 0.0000 — 1.0000 
1 1.0000 0.0000 — 1.0000 0.0000 
2 0.2500 0.0000 —0.5500 2.1000 
3 0.1111 0.0000 —0.3778 2.7167 
1 0 1.0000 0.0000 0.0000 — 2.0000 
1 0.0000 — 2.0000 —2.1000 —1.2000 
2 —0.8000 —0.6000 —3.4000 1.6500 
3 —0.6667 0.0000 —3.1962 3.5581 
2 0 1.0000 0.0000 0.2500 —4.5000 
1 0.0000 —3.2000 —3.2000 —4.0550 
2 —0.9000 —1.2000 —5.4942 —1.1538 
3 —0.9231 —0.3846 —6.0154 1.6216 
3 0 1.0000 0.0000 0.2222 —5.9157 
1 0.0000 — 4.0000 —3.2324 —6.3183 
2 —0.9412 —1.7647 —6.2154 —4.0216 
3 — 1.0667 —0.8000 —7.6534 —1.1200 


To determine g?), by use of punch cards we proceed as follows: for each number 
on an S card, we make one punch-card, and separate into groups. Thus, in order to 
compute the [Q, 2] stencil which has 8 rows we would need 8 groups of S cards. See 
Table 2. In the group [S, 1] there would be four cards, in [S, 2] eight cards, in [S, 3] 
twelve cards, in [S, 4] sixteen cards, in [S, 5] and in all succeeding groups twenty 


cards. 
We arrange the numbers in the first line of the [Q, 1] stencil (see Table 3) as 


follows: 


re 0 0 0 
qin gor 0 0 

G21 gin qo. 0 

qs gan qi qo 
981 G31 G21 gin 
G51 qs G31 Gan 
q . ; qs q a : 4 . : 


For each of these we punch one card and separate into groups as indicated. We do 
the same for each row of the [Q, 1] stencil, so that we will have seven such groups in 


total. (There are then four q$), q{'}, - - - cards and a total of six zero cards.) 
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TABLE 2. The values of i. for n=1, 2, 3, 4, 5.* 
n 1 J Re(s,) Im(s,") n | rcs i Re(s*’) (Im i.) 
—— - = ES eles = a = Z a a a a oe . —— = 
1 5 4 1.0000 0.0000 4 5 3 | 1.0000 | 0.0000 
5 3 1.0000 | 0.0000 4 | 3 | 0.0000 —2.0000 
5 2 1.0000 | 0.0000 3 3 | —0.9250 —0.6000 
5 i 1.0000 | 0.0000 213 | =1eu7 0.4125 
: i ne - - | 5 2 1.0000 0.0000 
2; 5 | 4 1.0000 0.0000 | 4 | 2 | 0.0625 —3.3750 
4 | 4 0.5000 —0.5000 3 | 2 | 1.2500 —1.6125 
5 | 3 1.0000 0.0000 2 | 2 | —1.6043 | 0.3846 
4 3 0.0000 —2.0000 5 | 1 | 1.0000 | 0.0000 
5 2 1.0000 0.0000 = &. | 0.0555 | —4.4792 
4 | 2 0.1250 ~—3.7500 3 | 1 | -1.2787 | —2.4619 
5 1 1.0000 0.0000 | 2 1 | ~—1.8205 —1.2054 
4 t 0.1111 —4.9583 = a —}— _ 
—_—}— an ios 5 s12 1.0000 | 0.0000 
3 5 4 1.0000 0.0000 4 4 0.8000 —0.2000 
4 4 0.6667 —0 .3333 3 4 | —0.0500 0.0000 
3 | 4 —0.2500 0.0000 2/4 —0.0655 0.4200 
5 | 3 1.0000 0.0000 1 4 —0.0755 0.5433 
4 3 0.0000 — 2.0000 5 3 1.0000 0.0000 
sia —0.9667 | —0.6000 4 | 3 0.0000 | —2.0000 
5 2 1.0000 | 0.0000 3: i 3 —0.9000 | —0.6000 
si3 0.0833 | —3.5000 2 | 3 | —0.9467 0.3300 
3 2 —1.3667 | -—1.7500 1 3 —0 .6392 0.7116 
5 1 1.0000 | 0.0000 5 2 1.0000 0.0000 
4 1 0.0741 —4 .6389 4 - 4 0.5000 — 3.3000 
3 1 —1.3912 | —2.6943 3 2 —1.1800 —1.5300 
|  —— ————___|——____— 2 2 —1.4681 | —0.3846 
4 5 4 1.0000 | 0.0000 1 2 —1.2031 0.3243 
4 4 | 0.7500 —0.2500 5 1 1.0000 0.0000 
3 4 —0.1250 | 0.0000 S-.45% 0.0444 —4.3833 
- is —0.1097 | 0.5250 3 1 —1.2112 —2.3225 
2 1 —1.6697 —1.1243 
1 1 —1.5307 —0.2240 
TABLE 3. The values of g°’,. 
m k Re (qo) Im (qo, ) m 4 Re(q.., ) | Im (q..% ) 
0 1 0.0000 2.0000 2 1 —0.5000 9.0000 
2 1.0000 0.0000 2 3.2000 | 4.0500 
3 0.3666 —1.4000 3 | 3.6628 0.7692 
4 0.1889 | —1.3583 4 | 3.0077 | —0.8108 
| 
1 | 1 | 0.0000 | 4.0000 3 | 1 ~0.4444 | 11.8333 
} 2 2.1000 | 1.2000 2 |. 3.2m 6.3183 
| 3 2.2667 | —1.1000 3 4.1436 2.6810 
4 1.5981 —1.7791 4 | 3.8267 | 0.5600 





* The s;;’s for n=6, 7, 8 have to be computed in a similar manner. 
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Starting with the group of cards for line one, we take the “dot” product for each 
aay he with the [S, 1] cards.* Let the numbers in the first row of the [S, 1] card 
be denoted by s{?, s{?, s&, s&); we then have to compute the following products: 


sO +40 -s8+0 s2l+0 -sf 
gi). st) + of!) 5 0 mo +0 -si 
Gan s qt: - + gon Sp2 +O +55, 
ales + ols + alli + als 


1) 1) (l a (1). (1) a). 
q a ) +4 ). si 


G41 ° 55,4 3,1 ° 55,3 92,1 ° 55,2 1,1 


The actual procedure is to form the products gf)- sf}, gi)-s&l, @}-s,- ++; 
0-sf3, gf)-sf&, g{)-s& and enter these together with appropriate signs etc. on the 
corresponding cards. These cards are then assorted in sub-groups as indicated and the 
products summed to yield the desired “dot-product.” 

To compute the numbers which are to appear in the second line of the [Q, 2] 


stencil, the numbers in the second line of the [Q, 1] stencil are arranged as follows: 


gh) 0 0 0 
q's qo.2 0 0 
2, q + se 0 
= qs'2 g!'2 gs. 
rots qs'3 qi'3 G53 


For each of these numbers one card is punched and separated into groups of eight 
cards as indicated below, using the cards previously punched for the first line of the 


[Q, 1] stencil: 


qo2 gf 0 0 0 0 0 0 
qi'2 gi} gi qi!) 0 0 0 0 
qh'2 Qn qi'2 Qi Go.2 Con 0 0 
G32 G3 Go'2 qo qiiz Gil Co Got 


The dot product of each of these groups in the [S, 2] cards are then taken. Thus: 


Go2 2 Sha + gas aa + 0 ‘ssa + 0 sa + 0 ssa + 0 2 +0 -sht +00 -st 


ds) + gf 62h + 5 + lo + a a + ols + old s8l + ofl a 
The actual products are again taken as indicated before, and then reassorted and 
summed. 

This process is repeated for each line of the [Q, 1] stencil until each group has 
twenty cards in it, after which the first line is discarded when a new line is added so 
that there are never more than twenty cards. In this way the “tentative [Q, 2] sten- 





(p-)) (n) 


* Note 2 that, in general, both g,, and nas are complex, and in crder to compute qe wee on punch 


card machines we have to use rules for the evaluation of complex numbers. 
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cil” is computed and it is necessary only to add (m+1)q®,,,, and multiply the result 
by —2/3 to obtain the actual values of the elements of the stencil [Q, 2]. This pro- 
cedure may then be repeated to obtain the values of the elements of the stencils 
[Q, 3], (ie. the quantities g®,) [Q, 4] ete. 

In Tables 2 and 3 the values of the elements of the cards [S, »] and those of the 
stencil [Q, 1] are given. 

7. Conclusion. Application of the differential analyzer to the construction of par- 
ticular solutions. As soon as a sufficient number of [Q, p] stencils (i.e. g®,) have 
been prepared we can, by (3.1) and (3.2) and standard methods of evaluation of 
polynomials on punch card machines, determine the values of p2,_; and pz, for a suff- 
ciently dense lattice of points. 

The second step, the determination of the coefficients a, in the expression (1.2) in 
order to obtain a solution of (1.1) which assumes the prescribed boundary. values, 
will be discussed in a subsequent paper. The basic idea of the procedure to be em- 
ployed has already been indicated in [2]; we shall however discuss this in more detail 
from the point of view of numerical analysis. 

Sometimes we need solutions of (1.1) in connection with questions other than the 
boundary value problem, and it is convenient to apply the method of particular solu- 
tions in the following slightly modified form. 

As was proved in [1], §1, 


U(z, 2) -{ 


l 
E(z, 2, t)f[(2/2)(1 — &) |dt/(1 — &)'/*, (7.1) 
1 


where f is an arbitrary analytic function of a complex variable, which is regular at the 
origin, is a solution of (1.1). 

There are instances in which a large number of solutions of the same equation 
are required, and the corresponding functions f are known. (See (7.1).) (This situation 
occurs, for example, if an “atlas” of solutions has been prepared.) In these cases it is 
then very convenient to employ graphical methods. In the following we shall indicate 
two graphical methods for the evaluation of (7.1). Both can be performed con- 
viently by use of a differential analyzer. 

I. One prepares once and for all for a given equation (1.1) diagrams in which the 
curves 


Y = E,(z,, 2,, #) = Re [E(z,, z,t)|, —1sS¢31, (%,, %,) fixed, —_ 
ee 4 
VY = E,(z,, 2, t) = Im [E(z,, 2, 2)], —1S5¢21, (2, 2) fixed, 
for a number of points (x, y) =(x,, y,), y=1, 2, 3, - - - are drawn. Further one has to 
prepare tables for the values }32,(1 —f), for tz,=—1, —1+a, —1+2a,---, 1, where 


ais a sufficiently small positive constant. z,=x,-+iy, denote the coordinates of points 
mentioned above. If now the function f(z) =u(z, 2)+iv(z, Z) is given, say in the form 
of two diagrams for curves u(z, 2) =const. and v(z, 2) =const., we draw (using the 
tables mentioned above) the curves 


Y = u[z(1-),211-)], Y= 201-),20-8)], -1s¢s1. (7.3) 


Using these diagrams and those mentioned above and employing a differential ana- 
lyzer (or simply an integrator) we compute the real part of (7.1), 
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t=1 
f | Er(z,, 2, d)u[3z(1 — &), 42,1 — 2)] 
t 


1 


— E2(z,, 2,, tv[4z.(1 — &), $31 — &2)] }de, (7.4) 


and analogously its imaginary part. 

II. Sometimes it is not sufficient to determine the values of (7.1) at a set of points 
(x,, ¥») which are prescribed in advance. Then one can apply the following procedure 
which was suggested to the author by Mr. Hans Kraft. 

One prepares (once and for all) diagrams 


E\(2, 2, t.) = const., Ee(z, 2, ¢.) = const., ¢, const. (7.5) 


for a set of values ‘(j= —1, —1+a, —1+2a,---, 1. 

Using these diagrams and the tables (described in method I, for every required 
value of z we can easily determine the curve (7.3) and evaluate the real and imaginary 
part of (7.1). 

Remark. The procedure I can be performed by the use of punch card machines. 
In this case instead of diagrams (7.2) it is necessary to prepare master cards. 

The author should like to thank Professor George E. Hay for his exceedingly help- 
ful advice and friendly criticism. 
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THE LIFT OF A DELTA WING AT SUPERSONIC SPEEDS* 


BY 


H. J. STEWART 
California Institute of Technology 


1. Introduction. The use of the two dimensional linearized theory of supersonic 
flows in the solution of airfoil problems as introduced by Ackeret! has been extremely 


successful in solving these problems and the results have generally been completely 
satisfactory for engineering purposes. The generalization of these results to the three 
dimensional finite span problems has, however, progressed rather slowly due to math- 
ematical complications. The flow near the tip of a rectangular wing was given (in- 
correctly) by Schlichting.? The drag of a “delta” wing (a wing having an isosceles 
triangle for its planform with the symmetric vertex pointing into the oncoming flow 
as in Fig. 1) has been determined by Puckett.* These two flow patterns and many 
other technically interesting finite span flow problems are particular cases of conical 
flows. A conical flow is one for which the fluid properties (pressure, velocity, etc.) 

are constant along each radial line emanat- 


MACH CONE AT THE ing from the given origin. The concept of 
TRAILING EDGE si ee ay a conical flow was given by Busemann‘ who 
\ developed certain general techniques for 


treating these flows and who applied the 

method to several problems including 

Schlichting’s problem. 

j The methods of analysis used by Buse- 
mann have, however, proved to be rather 
obscure, and it has been found difficult to 
follow these methods in the solution of 
additional conical flow problems, in par- 
ticular the currently very interesting prob- 
lem of the lift of a delta wing. A new 
Fic. 1. Delta wing in a supersonic flow. method of treating these conical flow air- 

foil problems which uses the well known 
theory of conformal transformation has been devised. It is the purpose of the pres- 
ent paper to discuss this method and to apply this method to the problem of the 
lift of a delta wing. In this application it is only necessary to consider the case for 
which the leading edges of the delta wing are within the Mach cone from the vertex. 

The other case for which the leading edges are outside the Mach cone has already 

been solved by Puckett. 

In the present method no essential mathematical difference is found in the solution 
of the two cases. 





* Received May 21, 1946. 
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3 A. J. Puckett, Aero. Sci. (To be published shortly). 

4 A. Busemann, Schriften der Deutschen Akademie fiir Luftfahrtforschung, 7B, 105 (1943). Also 
Luftfahrtforschung, 12, 210 (1935). 
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2. General theory of conical flows. It is well known that the linearized theory of 
steady supersonic flows is based on the Prandtl-Glauert equation, 


oP FF PF 
(1 — M?)- —= 
Ox® dy? = az” 


0, (1) 








where the undisturbed flow of Mach number MM is taken to be parallel to the x axis. 
Here, P may denote a velocity or acceleration potential, or one of the velocity compo- 
nents “, v, w in rectangular Cartesian coordinates x, y, 2, or a property of the state 
of this fluid such as pressure or enthalpy. It can be seen that the coordinate trans- 


formation 


x? + «eee R cos w 
pa) ——— —(¥ +8) = SS 








M? —1 p/MI— 1 
y? + 227? 2 
u= E — (M* — 1) —— = [1 — (M? — 1) tan? ow)", 4) 
x” 
6 = tan! (y/z), 


where R= (x?+y?+2?)"/? and w=tan~![(y?+2?)!/2/x |, transforms the Prandtl-Glauert 


equation, Eq. (1), into 


o°P oP re) oP 1 o@P 
. 2 — (2 = gf —— — ae ee os @, 3 
"ar . "ar . mak “>| ” P (3) 


The surfaces on which @ is constant are the meridional planes through the x axis; 
the surfaces on which yp is constant are cones about the x axis; and the surfaces on 
which r is constant are hyperboloids. It may also be noted that r=0 and w= & on 
the Mach cone through the origin. Both yw and r are real within the Mach cone and 
complex outside it. The harmonic solutions of Eq. (3) may be written in the form 


p=TA 4 BF sae: ae (4) 
= ee Q;(u)S Usin (mays ’ 
by the well known theory of the Laplace equation. Here, Py and Q? denote Legendre 
functions of the first and second kind, respectively. By introducing the normal spheri- 
cal coordinates as given in Eq. (2), Eq. (4) is seen to give the harmonic solutions of 
the Prandtl-Glauert equation in spherical coordinates. 

Busemann’s conical flows are included in the general solution of Eq. (4) as a special 
case. For example, if P is a velocity potential, then »=1. On the other hand, if P 
is one of the Cartesian velocity components (u, v, w), a property of the state of the 
fluid such as the pressure or enthalpy, or the acceleration potential, then n=0. It 
is the latter case which is of particular interest here, for P is then independent of r, 


and Eq. (3) becomes 
7) oP o*P r 
(n* — 1) — | @* — 1) — | + — = 0. (S) 
Ou Ou 06" 


It is apparent that this may be reduced to the Laplace equation in two dimen- 





sions;* in fact, if 





* This result was first communicated to the author by Mr. W. D. Hayes. 
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u—l 
5 = - ’ 
4 wt 
Eq. (5) becomes 
F) oP 0°P 
Ss (s — + = (), (6) 
Os Os ae? 


This is the normal form of the Laplace equation in two dimensional polar coordinates. 
It is seen that s is a function only of yu and is thus constant on any one of the cones for 
which w is constant. The relations between s and w are as follows: 


\/M? — 1 tanw 


s= oo a ON a eel (7) 
1+ 7/1 —-— (M? — 1) tan’w 
see as 2s 
\/ M? — 1 tanw = ——-- (8) 
1+ s? 


It may further be noted that s=1 on the Mach cone through the origin. 

Since the reduction to Eq. (6) is possible, any of the quantities which P may 
represent can be written as the real (or imaginary) part of an analytic function of 
the complex variable ¢ where 

¢ = se”. (9) 


Furthermore, all the methods of treatment of such functions, in particular the method 
of conformal transformation, may be used in the analysis of these quantities. If P is 


the harmonic conjugate of P and 





P+ iF = P(e), (10) 
the Cauchy-Riemann equations for these conjugate functions may be written 
oP aP ; oP 
AY a = say = (ue? — 1) . 9 
Os 06 Ou 
- nf (11) 
oP oP. aP 
— — = s— = (y? — 1) — 
06 Os Op 


In the direct airfoil problem, the airfoil geometry is given, and if the 2 axis is 
taken normal to the airfoil plane, the boundary conditions for determining the flow 
are thus given in terms of the disturbance velocity component w. It is desired in this 
case to compute the pressure distribution which may be easily expressed in terms of 
the axial disturbance velocity component w. In the inverse airfoil problem, a pressure 
distribution is defined, and it is desired that the airfoil shape be computed. In either 
case the boundary condition is given in terms of one velocity component and another 
velocity component gives the desired result. For a conical flow there are simple rela- 
tions between the complex functions representing the various Cartesian velocity com- 
ponents. The use of these relations is the essence of the present method of treatment 
of conical flows. These relations between the complex functions corresponding to the 
Cartesian disturbance velocity components u, v, w which will be written 
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u+ in = U(e), v+ id = V(¢), w+ iw = Ws), (12) 
are essentially the vorticity relations. 


The fundamental linearized relations governing the steady flow of a fluid at super- 
sonic speeds are the vorticity relations 


Ov Ow ou Ow Ou av H 
—=—; (13) —=-—-; (14) —-s—) (15) 
Os Oy Oz Ox Oy Ox 
and the linearized equation of continuity 
Ov Ow Ou 
-- +— = (M* — 1)—- (16) 
dy as ax 


If these are solved simultaneously, it is easily seen that each of the velocity compo- 
nents obeys the Prandtl-Glauert equation, Eq. (1). For any conical flow each of the 
velocity components must be a function only of the coordinates y and 6. By means of 
this fact, Eqs. (13) to (16), respectively, may be written as follows: 


ov : Ow Ow 2 Ov ‘i 

u(w? — 1) (cos 4 — sin @ —]} = ( cos @-—— + sin @ — }, (17) 
Ou Ou a0 a0 

(u2 — 1)9/2 dw Ou : ou ; 

= —— = p(w? — 1) cos 6 — — sind —; (18) 
VM? —1 op Ou 00 
(u2 — 1)3/? dv ; Ou Ou 

as —— = p(w? — 1) sin 6 — + cosd—» (19) 
\/M?—1 op Ou 00 


Ou ; Ov Ow 
— \/M? — 1 (yp? — 1)8/?? —- = p(e? — 1) (sino —— + 98 — 


Ou Ou Ou 
ov : Ow 
+ [cos @ — sin 0 ——}. (20) 
06 00 


If Eqs. (18) and (19) are combined, it is seen that 


Ou (u? — 1)3/? Ov . Ow 

= — cos 6 — sin 6 , (21) 
00 VM? —1 Ou Ou 
Ou (u* — 1)}/* . Ov Ow 

= — sin 0 +. cos 6 —— }. (22) 
Ou \/ M? — Ou Ou 


Furthermore, Eqs. (20) and (22) show that 
w—1/. ov aw dv Ow ; 
— sin 6 + cos @—-} = — {cus 8 — sin 0 ——}. (23) 
m Ou Ou 0 00 


If the derivatives with respect to w are eliminated from Eqs. (17) and (23) by 
means of the Cauchy-Riemann equations [cf. Eq. (11)] for v, w, i and #, these equa- 


tions may be written as follows: 
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1 d ov Ow 00 Ow 
— {sin 6 — + cos 6——]} — [ cos 0 — sin 6 — ) = 0, (24) 
7 06 06 06 00 
1 : Ov Ow f Ov Ow 
- sin @ og + 08 0) 4 (cos @ — sin 6 -) = () (25) 
be 06 06 00 00 
If Eq. (25) is multiplied by i and added to Eq. (24), it is seen that 
1 ‘ ov ow av _ OW 
sin 6 —— + cos 0 =) +i (cos @ —— — sin 6 —} = 0. (26) 
06 06 00 00 


jm 
Since V and W are functions of the complex variable ¢, this may further be written 


dV dW iu sin 6 — cosé 





a oe (27) 
dt dg sin #@+ ip cos 0 
and, by the definition of s and Eq. (9), 
dV 1-77 dW 
—— = j——_—_ —_ (28) 


= 1 . 
dg 1+ df 
A similar treatment of Eqs. (21) and (22), V being eliminated by Eqs. (17) and (27), 


shows that 








dU 2¢ dW ’ 
—_ = — oe a (29) 

d¢ (1+ ¢%)VM?-1 dt 
These two relations, Eqs. (28) and (29), are 

Wis a Sie : 
- the fundamental relations for the present 
ie lites nee: ee : 
Pal P an treatment of conical supersonic flow prob- 
o <az0 i" lems. 

/ \ 3. Example. Lift of a delta wing. The 
/ \ general techniques developed in the previous 
ime i sidan Ff =i} section will now be used to compute the lift 
\ is, -is, of a delta wing at a small angle of attack 
\ / for the case in which the leading edges are 
\ / inside of the Mach cone (see Fig. 1). The z 
* rw axis is taken normal to the airfoil. The con- 
“oe ditions in the ¢ plane are shown in Fig. 2. 
Note that the airfoil cuts the ¢ plane on the 
Fic. 2. Boundary conditions in the ¢ plane. imaginary axis. The boundary conditions for 


determining the vertical velocity w are then 


w=0 on s=1, 
(30) 


w= wo = — Ua on the airfoil, 


where U is the velocity of the mean flow and a is the angle of attack of the airfoil. 
This boundary value problem can be solved by conformal transformation. First, 


apply the transformation 
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This maps the interior of the unit circle in the ¢ plane into the entire §; plane with 
the region Ref >0 corresponding to the region Jmf,>0. This transformed plane is 


| wmf, 


wew, > wo + hadith Ril, 








Fic. 3. Boundary conditions in the ¢ plane. 


shown in Fig. 3. The points at the wing tips, ¢ = +iso, are transformed into the points 


(:= +1/k where 





250 en 
k=- — ee * M? — 1 tan Wo. (32) 
1+ s 
Second, apply the transformation 
f d 
be -f ) = ——, (33) 
0 VA oC — Re) 
i.e. 
f1 = sn(f2), (34) 


where the elliptic function has the modulus k. Then the region Jmf{,>0 is mapped 
into the rectangle having its corners at {2= +K, i1K’+K where K and K’ are the 
complete elliptic integrals of the first kind having a modulus of k and k’ where 


1:— si 


1+s2 


k’ k? 


=Vi- (35) 


By integrating around the slit from —1 to 1 in the ¢; plane, it is seen that the region 
Imé, <0 maps into the rectangle having corners at {:=2K+K, 2K+K+iK’. 

Now, the transformation given by Eq. (31) is double valued, i.e., two points in 
the ¢ plane correspond to each point in the ¢; plane. The ¢; plane must thus be con- 
sidered as a two sheeted Riemann surface with one sheet corresponding to the interior 
of the unit circle in the ¢ plane and the other sheet corresponding to the exterior of 
the unit circle in the ¢ plane. Furthermore, the value of the downwash velocity w 
must be equal and opposite at inverse points in the ¢ plane. This permits the analytic 
continuation of w throughout the entire ¢ plane; in particular it is seen that w= —wo 
on the exterior points corresponding to the airfoil. The two sheets in the ¢, plane 
are connected through the slit from —1 to +1. A contour cutting this line passes 
from the upper to the lower sheet or vice-versa. The sheet which corresponds to the 
exterior region of the ¢ plane is thus seen to be mapped into the rectangle having 
corners at (=K+2K, K+2K-—-iK’. The entire plane is mapped into a basic rec- 
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tangle in the ¢2 plane as shown in Fig. 4. As ¢; has periods of 4K, 21K’ [see Eq. (34) ] 
in ¢2, this pattern is repeated throughout the ¢ plane. 


SP 
T 


Pld — 
9 mR 2K 


{fe | 


ix’ 




















+—¥ 





Fic. 4. Boundary conditions in the ¢ plane. 


The function dW/df2 (but not W itself) must be doubly periodic in the ¢2 plane 
with periods 4K and 27K’, the first corresponding to a loop around the points (:= +1 
and the second corresponding to a loop around the points {;=1, 1/k or —1, —1/k. 
The only singularities of W or dW/d¢é2 must be at the points corresponding to the air- 
foil leading edges, i.e., at the points conjugate to iK’+K. Finally dW/df2 must be 
pure imaginary on the lines Imf.=nK’' and Ret.=K+2nK (n being any integer). 
All of these conditions are satisfied by the Jacobian elliptic function 


dW 


= 1Dcd?"(¢o), (36) 
dfe 





where v is any positive integer and D isa real constant. If this is integrated it is seen 
that for n>0, W has a pole of order 2n—1 at the wing tips. The cases for n>1 can 
then be discarded as the singularity at the wing tips is seen to correspond to a source- 
sink complex which has an infinite total lift. Furthermore, the case for »=0 may be 
discarded as [see Eq. (29) | it requires that U(¢) have a logarithmic singularity on the 
Mach cone. The appropriate solution is thus 


dw 





— = iDcd*(¢2). (37) 


df» 


The constant D may be evaluated from the fact that 


iK’ qW 
w = RI} f —— dts (38) 
0 dfs 


If this integration is carried out, it is seen that 


9 
“909 


“'0 


Dae ees (39) 
E(k’) 
where E(k’) is the complete elliptic integral of the second kind having a modulus k’ 
as given by Eq. (35). 
If the variable ¢ is eliminated from Eq. (37) by means of Eq. (31) and (33), it is 
seen that 
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dW 2wy 1 Pe) 
ARNE SORDID oA} aC am (40) 


dg kE(R) Ty 1\ 92 
G + 59) (*: 7 =)| 
so 


Thus, from Eq. (29), 
dU 4+wo cl + ¢) 


—= : ——— — ——___— ——— - (41) 
d¢ RE(R')\/M?—-1 F | f 1\7 
(¢? + 59) (« v :) 
s* 
Since 'U =0 at ¢=1, the integral of Eq. (41) is 
Rw; C2 — 1 
= Te. Seed... er (42) 
E(k’)\/M? — 1 Ati ee. Se 
[ « + $9) (+ a alee 
So 
On the top side of the airfoil =m where —so <9 <5o, So 
kw 1 ? 
= ° ee 5. (43) 





E(R)/M?-1 7... (1 \7e 
E — 7°) (- - ") 
s5 


This result may be considerably simplified if we introduce [from Eq. (8) and (30) | 


k = \/M? — 1 tan wo 





(44) 
w= - Ua 
and 
tan w 
ee ; 
tan wo 
Equation (43) then becomes 
4/u 4 tan w 
(=) ee (43) 
a\U E(k’)V//1—-? 


The slope of the lift curve dC,/da is given by the mean value of 4/a(u/U) over 


the surface of the wing; thus 


Oe ites a cad oven (46) 
da E(R’) Jo V1I-?f 
and, by Eq. (44), 
best _ 2m ° (47) 
da E(k’) 


In the limit for which wy or so0, k’1; so E(k’)—>1. For this case which was given 
by Jones® 


*R.T. Jones, N.A.C.A., Technical Note 1032 (1945). 
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| \ TWO- DIMENSIONAL wine | | | 
ei ee ee ee ee. 
ONES’ THEORY—, | et. | [__-<— PUCKETT'S_THEORY 
” (eq. 48) | ae er (Eq. 49) 
= Se. ia | | a 
; a — DELTA WING | 
m2-1 46, 4, | (Eq. 47) -}-; 
4 da [ i 
ae > és «3 | lial 
} | 
° | | 
0 0.2 a4 06 as 1.0 1.2 is 
ne (m2-1 tan we 
Fic. 5. Lift of a delta wing. 
2. 
TI 
\—--+— —— 
\ | 
1.6 }—_+—_ +__+\++—_++- a oe - 
: pe. 
Fa | -- \ /_TWO - DIMENSIONAL WING 
+ | 
JONES’ THEORY 
Ler\ @ 92 10° — —_—_}—_—__4 
~~ 
4G, ! } Ns 
da | | ; | ‘ee 
= | | ELTA winc—’, —— —% 
| | tS {4 
| r ; | PUCKETT'S THEORY 
ee Se See Eas Cee Oe er. 
| | 
2 a Se ee 
eS |_| 
0 
' 2 3 ~ 5 & 7 
i 
Fic. 6. dC, /da vs. M for a delta wing with w =10°. 
dC, 
—— = 2r tan w. (48) 
da 
On the other limit for which sop—1, k’-0; so E:(k’)—>7/2. For this case 
dC, 4 
= = 4 teh 6 = - (49) 
da V/M?*—-1 


This limit, the same as the two dimensional solution, had previously been obtained 
by Puckett. 

It may further be noted that the quantity }./M?—1 dCz/da is a function only 
of the parameter k = \/ M?—1 tan wo. This result is shown graphically in Fig. 5, and 
the slope of the lift curve for a particular case, w)=10°, is shown as a function of 
Mach number in Fig. 6. 
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LINEARIZED SUPERSONIC FLOWS WITH AXIAL SYMMETRY* 


BY 


WALLACE D. HAYES** 
California Institute of Technology 


1. Introduction. The study of spatial linearized supersonic flow may be aided by 
the study of some simple fundamental flows with axial symmetry. Through the prin- 
ciple of superposition, these flows may be combined to give more general flows about 
various objects and about lifting systems. It is the purpose of this paper to express the 
equations of linearized supersonic flow in a system of conical coordinates, to develop 
a theory for fundamental flows with axial symmetry, and to describe examples of 
such flows and of their combination by superposition. 

Various examples of the fundamental equations and solutions here described will 
be given in later papers, together with the development of some concepts useful in 
this field. 

2. The velocity potential. Steady-state compressible irrotational flow can be de- 
scribed by a velocity potential ¢ whose gradient is the velocity vector. Under the as- 
sumption that the velocity deviations from a uniform supersonic flow of the Mach 
number M are small, the differential equation for this potential takes the linear 
form':? 


1 1 
re + — or + — 00 — (M? — 1)¢., = 0 (1) 
r rT 


in cylindrical coordinates. 

The fundamental uniform flow is given by the potential ¢)= Vz where V is the 
velocity corresponding to the Mach number M. Equation (1) will be considered as 
yielding velocity deviations which must be added to the velocity of the fundamental 
flow to describe the net flow. 

A new coordinate is introduced to replace the coordinate r: 


t = (r/2)\/M? — 1. (2) 


This quantity is the ratio of the tangents of the polar angle and of the Mach angle. 
Equation (1) with r eliminated and ¢ introduced becomes 


1 1 
(1 — Pde + a (1 — 2#)@.+ pe + 2tzbiz — 2°ber = 0. (3) 


By separation of variables a solution of the form 


* Received March 25, 1946. 

** This paper was prepared while the writer was employed by the Lockheed Aircraft Corporation. 

1 R. Sauer, Theoretische Einfiihrung in die Gasdynamik, Springer, Berlin, 1943. Reprinted by Ed- 
wards Bros., Ann Arbor, 1945. 

2G. I. Taylor, and J. W. Maccoll, The mechanics of compressible fluids, in Durand, Aerodynamic 
theory, vol. 3, Berlin, 1935. 
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N 


o = s"@(t, 8) (4a) 


or 


@ = 3" sin (m6 + B)T(2) (4b) 
is found. The function ® satisfies the equation 


1 1 
(1 — #)@i, + 7 (1 + 2(n — 1)f)%, — n(n — 1) + — Dog = 0 (5) 
t 


and may be called the velocity potential for generalized conical flow. If m=1, the 
function ® describes conical flow. The function 7 satisfies the equation 


1 1 
1Q—-@)Taet ri (1 + 2(n — 1)#)T, — rs (m? + n(n — 1)#)T = 0. (6) 


Superposition of solutions of the type of (4a), (4b) will give a general solution. 


The velocity components are 





| : 
eee (7a) 
in the radial direction, 
Ma : 
v= - —— go (7b) 
tz 
in the azimuthal direction, and 
Y — 
w= o: — — > (7c) 


in the axial direction. The pressure in linearized supersonic flow is given in terms of 


the velocity components by 








: u? + v* 
p - —— p J Ww + a , (8a) 
and the pressure coefficient by 
: w u+o° 
C,=-—2 + — : (8b) 
V 2" 


The part of Eqs. (8a), (8b) in « and v is not necessary if w is of the same magnitude 
as u and v. In many important cases, however, u*+v? is of the same magnitude as Vw 
and Eqs. (8a), (8b) must be used in its complete form. In these cases the validity of 
the solution should be checked. 

The singularity of (5) or (6) at ¢= +1 corresponds to the two Mach cones ex- 
tending from the origin in the three dimensional flow. Various ranges of ¢ correspond 
to various regions of flow, as shown in the following table. 
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Range of t 


O<t<! 
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Region of Flow 


inside downstream cone 
—1</s0 inside upstream cone 
1<i<o 
> outside both cones 
—x<i<-—1} 


3. Solutions of the differential equation (6). The parameter n is restricted to in- 
tegers and the parameter m to non-negative integers. The differential equation for 7, 
Eq. (6), has regular singularities of exponents (+m, —m) at t=0, (0, n+}4) at 
t= +1, and (—n, 1—n) at t= ~. The solutions about the origin are 


—n+m —n+mri 
T = ("F{- : y ene ¢ le; r) (9a) 
‘n+m+1 nu+m+2 
= #m(1 — aa | Las eet 1+ m; r), (9b) 
—-nu—m —unu—-m+i1 
T=t “P( — +. ——_______- 1—m; r) (10a) 
2 2 
we~-aeti e«e-ase 
=f~(1— aa — 7 1—m; r), (10b) 


where F denotes the hypergeometric function. The solution of negative exponent, 
Eqs. (10a), (10b) is not well defined. 

It is of considerably more value to express the solutions about /*=1, since then 
both solutions are well defined and two distinct types of solution may be distin- 
guished. One type of solution, designated as type I, is the solution of zero exponent 
at #?=1 and is real throughout the range of t. The resulting solution for ¢ has no singu- 
larity on the Mach cones. The other type of solution, designated as type II, is the 
solution of exponent n+} at ¢?=1 and is real only for ¢?<1 or only for t?>1. The re- 
sulting solution for @ is defined only within the Mach cones or only outside the Mach 
cones. These solutions may be expressed as follows: 





—-n+m —un+m++i1 
I) T = (°F ( ——— » ——————— ; —-na+};1- r) (11a) 
2 2 
—n—m —n—-m+1 
= f “F (- : > ———————— ; — nt}; 1-F), (11b) 
2 2 
n+m+t+iutm+2 
IT) T=t"(1 — #)"*3F (= "iain weal sat#;1-— r) (12a) 
n—-m+1Hn—-m+2 
=o "(1 — @*) "TIF ( < of ns +2;1- r). (12b) 
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Three special cases are distinguished according to the relative values of m and m: 


case A: —2x <nS—m-—l, 
case B: —msnsm-—1, 
case C: msn<o 


The distribution of these cases for small values of m and n is shown in the table: 


bars —3 —2 —1 0 1 2 

m tac ag : Cn ea — 
0 1 1 A Ce C s 

1 1 1 B B C t 

2 1 B B B B ¢ 

3 B B B B B B 


From a consideration of Eqs. (9) to (12) the forms of the two types of ‘solutions 
in the various cases may be found. For all solutions except solutions I-A (i.e., solutions 
of type I in case A) and solutions II-C, the form is explicit in terms of a polynomial 
in ¢? or in (1—/?). Solutions I-A and II-C have logarithmic singularities at t=0 and 
are discussed later. The polynomial forms are expressed as follows: 





Order of P(t?) Equation for 


Solutions Form whichever is integral of Calculation 
I-A logarithmic ——_—_——— a 
II-A tm (1 —¢2)"+4P (12) 4(—n—m-—1) or }(—n—m—2) (9b) or (12a) 
1-B t-"P(¢*) 3(n+m) or 3(n-+-m—1) (10a) or (11b) 
II-B t-™(1 —12)n+ 1 P (72) }(—n+m-—1) or }(—n+m-—2) (10b) or (12b) 
I-C im P(t?) 3(n—m) or 3(n—m-—1) (9a) or (11a) 


II-C logarithmic —— 





There is a connection between the solutions of Eq. (6) and the Legendre functions 
with the same values of and m, except that —n—1 is used when 7 is negative. How- 
ever, since Legendre functions are customarily defined only for m<n or mSn—1, 
respectively, they are of assistance here only in cases A and C. The polynomial solu- 


tions are 


ll 


(1 — )"? Pra — #7], (13a) 
a — #)*? P*[(1 — #7). (13b) 


I-C) r 
IT-A) 4 


II 


These solutions are most easily obtained by transforming Eq. (1) into Laplace’s equa- 
tion by introducing the variable iz/./M?—1. 

The polynomial solutions may be obtained in another form from an expression 
given by A. R. Forsyth,’ and the law for differentiating the hypergeometric functions. 
When 1 is not positive, these solutions are 


? A. R. Forsyth, A treatise on differential equations, 6th ed., Macmillan, London, 1933, p. 235. 
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on. 





1) T=(t" menos EEE of. Y1—-# “ore t (§ wm > jee iE “nt m |. (14) 
1" K ) Apne 
U 1 EN, pny {—P)-*F "| (15) 
) T -_ it ew = V eo ’ —ntm —-\ is Pe —n+m x 5 
d(t*)-* 


and when +1 is not negative, 


n+l 
1) T = #™(1 — #)"+4 —_—__ [(14 1 — P)etit™ —(1i-— [<= {2) n+ 1F m ; 16) 
) d(t?)"*) | (16) 
d*ti eee ome 
Dm Teri — H—4——— (i + Vi = ee + (1 = 1 = Fr}. 17 
d(t?)"*} K | ( ) 


4. Logarithmic solutions. The logarithmic solutions I-A and II-C are most easily 
expressed in terms of the Legendre functions, as in Eqs. (13a), (13b). They are 


I-A) T=(1-#)"%0" J[a—+#)*], (18a) 


II-C) T=(1-?)""or (a —#)*). (18b) 
Since the validity of Eqs. (14) to (17) does not depend upon m being an integer, 
and since an appropriate solution of the form of these equations vanishes as a log- 


arithmic solution is approached, the logarithmic solutions may also be obtained by 
differentiating such solutions with respect to m. The logarithmic solutions in this form 


are 
q* pees a uemmiaiae 
I-A) 7 = te —— |[(114+ V1 — #)-**™ log (1+ V1 —-P) 
d(#)-* 
+ (1— V1 — #)-"*™ log (1 — V1I—P)], (19) 
qvti =: et ae 
I-C) T= #1 — #)**! — [1+ V/1—#)"*¥"log 1+ V1-f 
d(é?)**! [( \ 8 ) 
+(1— JM1— #)**F™ log (1 — V1 — #)]. (20) 


5. Generating equations. The fundamental equation, Eq. (1), expressed in Car- 
tesian coordinates is invariant under differentiation with respect to any of these co- 
ordinates. Solutions of the type of Eq. (4b) expressed in Cartesian coordinates and 
differentiated with respect to these coordinates are still solutions of Eqs. (1) and (3). 
This fact permits a given solution of parameters m and m to yield solutions of parame- 
ters m—1 and m, m+1, or m—1: 


d 
T(n — 1, m) = nT — tT, = — * (¢*T), (21a) 
( 
m d - 
T(n —1,m+1) =—T-T,= —t#"—(t"7), (21b) 
t dt 
m d 7 
T(n — 1,m— 1) = = T+ T,.= +t sh ill (21c) 
¢ 


The procedures yielding these new solutions can be considered procedures of super- 
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position; for example, one solution superposed on its negative an infinitesimal dis- 
tance downstream yields the new solution given by Eq. (21a) 

These equations are not to be considered recurrence relations, as no system has 
been established for specifying solutions with respect to the multiplicative arbitrary 
constant. 

In a similar manner solutions with the parameter m increased by 1 may be ob- 
tained by reversing Eqs. (21a), (21b), and (21c) with suitable integrations. 

6. Integral relation. An integral relation connecting two solutions whose parame- 
ters differ in value may be obtained either from the corresponding relation for the 
Legendre functions or directly from Eq. (6). If 7; denotes a solution corresponding 


to m, and m,; and 7, a solution corresponding to m2 and mp, the relation is 


(m, — m’) J C1 — #)—2ataet 7 Todt 


b 
+ (1 — me)(my + no + 1) [ HA — £2) atnet)) 7 Todt 


a 


: dT, dT. 
dt dt 


b 
+ (ny = ne) t?(1 —_ f?)—h Gut OTT ° (22) 


a 


Setting m= or m, =m, we obtain simpler equations as special cases which may be 
used to obtain orthogonality relations between solutions. 

7. Two-dimensional cross-flow. The solutions of type I for which m=|n]| are 
given by 7 =". The corresponding solutions for ¢ in cylindrical coordinates are 


@ =r" sin (| #| 4+ 8). 


These solutions give two-dimensional cross-flow because of the fact that 2 does not 


appear. This cross-flow, as a result of the linearizing assumptions, appears as an in- 
compressible flow. 

8. Conical flow. The solutions of either type for which n=1 give solutions of 
conical flow, of which only those of type II are here treated. Since the exponent of 
these solutions at /?=1 is 3/2, both the potential and the velocity vanish on the Mach 
cone. The first few solutions are: 
m=()) /1—#—tanh" 1/1—P, 
m=1) t—4/1—-2—t tanh"! /1—2, 

m = 2) t-2(1 —1?)3/?, 
m = 3) t-#(1 —#?)*/2, 


The flow about an infinitesimal circular cone at zero incidence is given by the first 
solution (II, 1, 0), the solution of type II with n=1 and m=0. The flow about the 
same cone at a small angle of attack is obtained by superposing solutions (II, 1, 1) 
and (I, 1, 1) with appropriate constants on the (II, 1, 0) solution. A standard treat- 
ment of this case will be found on pp. 46 to 49 of Sauer’s book." 

9. Infinitesimal horseshoe vortices. An infinitesimal horseshoe vortex can be 
represented by a semi-infinite line dipole in the same manner as a planar vortex 


1946] LINEARIZED SUPERSONIC FLOWS WITH AXIAL SYMMETRY 261 


system can be represented by a planar dipole system. Thus a lifting element and other 
lifting systems can be represented by solutions of type Il with m=1, as shown in the 


following table. 








} 











Solut; Semi-Infinite Line | Designation in Terms of 
cena Dipole of Strength | Lifting Properties 

(I, —1, 1): ca —#)-3 Constant “Lifting element” 

(11,0,1): t(1—e)t+4 Proportional tos | “Lifting line” 

(ii, 1, 1): ‘(1 —#2)+4—¢ tanh (1-4) +3 Proportional to 3? “Lifting infinitesimal triangle” 


The “lifting element” solution, since the potential has exponent — 4 and the veloci- 
ties —$ at /?=1, has a troublesome singularity on the Mach cones. A simpler sin- 
gularity has the “lifting line” solution, whose potential vanishes and whose velocities 
have exponent —} at t?=1. When these solutions are superposed to give lifting sys- 
tems of finite dimension, the singularity in the velocities usually disappears. The third 
solution is the same as the one which gives the lift on an inclined circular cone. Ex- 
amples of the superposition of such solutions to form a lifting system will be given in 
a later paper. 

10. The acceleration potential. Since the axial velocity component w is a deriva- 
tive of the velocity potential in Cartesian coordinates, it satisfies the same equations, 
Eqs. (1) to (6), as does the velocity potential. The acceleration potential, whose fun- 
damental theory will be found in a paper by L. Prandtl,‘ equals — p/p for linearized 
flow, and also equals Vw+43(u?+v?) from Eq. (8a). Hence the approximate accelera- 
tion potential defined by y= Vw satisfies Eqs. (1) to (6). The relation of this quantity 
to the velocity potential for a given fundamental flows is of the type which leads to 
Eq. (21a), and hence the corresponding acceleration potential is given by a solution 
with m decreased by 1. It is important to note that this does not give the true linear- 
ized acceleration potential where u?+v? is not of smaller magnitude than Vw. Thus 
the two-dimensional cross-flow described above has no approximate acceleration po- 
tential, and the acceleration potential is given incorrectly in the vicinity of the axis 
for other flows. However, the true acceleration potential may not itself be superposed, 
and often the difference between the approximate and true linearized acceleration po- 
tential disappears under superposition. 

The “lifting element,” “lifting line,” and “ lifting infinitesimal triangle” have ap- 
proximate acceleration potential solutions (II, —2, 1), (II, —1, 1), and (II, 0, 1), 
respectively. With conical flow, the approximate acceleration potential is a function 
only of ¢ and @ and can be shown to satisfy Laplace's equation in two dimensions. 


‘L. Prandtl, Theorie des Flugseugtrag fliigels in susammendriickbarem Medium, Luftfahrtforschung, 
13, 313 (1936). 
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THE OPENING OF A GRIFFITH CRACK UNDER 
INTERNAL PRESSURE* 


BY 


I. N. SNEDDON (University of Glasgow) AnD H. A. ELLIOTT (University of Bristol) 


1. The determination of the distribution of stress in the neighbourhood of a 
crack in an elastic body is of importance in the discussion of certain properties of the 
solid state. The theory of cracks in a two-dimensional elastic medium was first de- 
veloped by Griffith who succeeded in solving the equations of elastic equilibrium in 
two dimensions for a space bounded by two concentric coaxial ellipses; by considering 
the inner ellipse to be of zero eccentricity and by assuming that the major axis of the 
outer ellipse was very large Griffith then derived the solution corresponding to a very 
thin crack in the interor of an infinite elastic solid. Because of the nature of the co- 
ordinate system employed by Griffith the expressions he derives for the components 
of stress in the vicinity of the crack do not lend themselves easily to computation. 
An alternative method of determining the distribution of stress in the neighbourhood 
of a Griffith crack was given recently by one of us? making use of a complex stress- 
function stated by Westergaard.* This method suffers from the disadvantage that the 
Westergaard stress-function refers only to the case in which the Griffith crack is 
opened under the action of a uniform internal pressure; the stress-function correspond- 
ing to a variable internal pressure does not appear to be known. 

In the present note we discuss the distribution of stress in the neighbourhood of a 
Griffith crack which is subject to an internal pressure, which may vary along the 
length of the crack, by considering the corresponding boundary value problem for a 
semi-infinite two-dimensional medium. The analysis is the exact analogue of that for 
the three-dimensional “circular” cracks developed in the previous paper? except that 
now we employ a Fourier cosine transform method in place of the Hankel transform 
method used there. A method is given for determining the shape of the crack resulting 
from the application of a variable internal pressure to a very thin crevice in the in- 
terior of an elastic solid, and for determining the distribution of stress throughout the 
solid. The converse problem of determining the distribution of pressure necessary to 
open a crevice to a crack of prescribed shape is also considered. As an example of the 
use of the method the expressions for the components of stress, due to the opening 
of a crack under a uniform pressure, are derived and are found to be in agreement 
with those found in the earlier paper.’ 

2. We consider the distribution of stress in the interior of an infinite two-dimen- 
sional elastic medium when a very thin internal crack —cSySc, x=0 is opened un- 
der the action of a pressure which may be considered to vary in magnitude along the 
length of the crack. For simplicity we shall consider the symmetrical case in which the 
applied pressure is a function of |y| but the analysis may easily be extended to the 


* Received March 12, 1946. 

1 A. A. Griffith, Phil. Trans. (A) 221, 163 (1921). 

2 I. N. Sneddon, Proc. Roy. Soc. (A) (in the press). 
7H. M. Westergaard, J. Appl. Mech. 6, A49 (1939). 
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more general case in which there is no such symmetry. The stress in such a medium 
may be described by three components of stress oz, ¢, and Tzy; the corresponding 
components of the displacement vector will be denoted by u, and u,. The differential 
equations determining the stress-components are‘ 

OC, We OTsy Oly 


= 0, (1) ——* + 


+ — = 0), (2) 
Ox oy Ox oy 


The boundary conditions to be satisfied are that all the components of stress and of 
the displacement vector must tend to zero as x?+ yy? tends to infinity, and that 


Try = 0, oz = — p(y), (3) 
when x«=0 and —cSySc. 

It is obvious from the symmetry about the axis x =0 that the problem of determin- 
ing the distribution of stress in the neighbourhood of the crevice is equivalent to that 
of determining the stress in the semi-infinite elastic medium x 20 when the boundary 
x=0 is subjected to the following conditions: 

(i) t2y=0, for all values of y, 

(il) ¢,-=— p(y), | y| Sc, 

u.=0 \y| —e. 
From the symmetry about the second axis y=0 we may take as solutions of the 


elastic equations (1) and (2) the expressions? 


9) oo 
%?~"- f $(p)(1 + px)e** cos pydp, (4) 
us 0 
2 i 2) 
ry = — J Blo\(1 = pade-r* cos pyd, (5) 
TT 0 
(ae aie 
Try = f po(p)e~** sin pydp. (6) 
T 0 


These expressions satisfy the equations of equilibrium and the boundary condition (i) 
above; the function ¢(p) is determined from the set of conditions (ii). The compo- 
nents of the displacement vector are similarly found to be 


21+ 0) ¢” cos py 7 
‘= — —— f b(p)e-**{ 2(1 —«) + px} ——dp, (7) 
rE 0 p 
21+ 0) ¢” sin py 
ly = ——— f b(p)e-"*{ (1 — 20) — px} ——— dp. (8) 
rE 0 p 
When x =0, equations (4) and (7) reduce to 
2” | 
¢;= f d(p) cos pydp, (9) 
Jo 
4(1 — 0”) ¢” cos py 
“.=— ( — f o(p) ———— dp. (10) 
rE 0 p 


‘A. E. H. Love, The mathematical theory of elasticity, 4th ed., Cambridge, 1934, p. 208. 
5 I. N. Sneddon, Proc. Cambridge Phil. Soc. 40, 229 (1944). 
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If we insert the boundary conditions (ii) into Eqs. (9) and (10) and make the substitu- 


tions 


“er -{& a 
p = &/c, y = ne, g(a) = — e( — P(ne), @ =) = £'/"F(£), (11) 
2n . Cc 


~ 


we obtain a pair of “dual” integral equations 
x 


f EF(E)J_1/2(én)dé = g(n), O<n<l1| 
0 


x 


f F(&)J—1/2(&n)dé 


for the determination of the function F(£). Once F(£) has been found, ¢(p) can be 
written down and the components of stress calculated by means of Eqs. (4), (5) 


Q, n> 1 


and (6). 
3. The dual integral equations (12) are a special case of a pair of equations con- 
sidered by Busbridge;* the solution may be obtained by substituting a=1, y= —1/2 


in the general solution given in the paper.® In this we obtain 


“5 
F(é) = 4/ e[ sae) f yl2(1 — y?)!/2o(y)dy 
T 0 
] 1 
+ ef u'/2(1 — u?)' "du f g( vu) vy? “ndeyidy | (13) 
0 “” 0 


Thus if the pressure p(y) is given by a Taylor series of the form 


p(y) = po dia, (<). (14) 


n=O 


convergent for ~cSy Sc, then the corresponding expression for ¢(p) is readily found 
to be 
= T(gn + 3) 


; , 
- ° ° j ! +2? 4 { 5 
o(p) = — dpoc*x! >>. a, \Jo(ep) + cp yt J (cpy)ave . (15) 
n=0 I'( In + 2) \ 0 f 


Substituting for ¢(p) from Eq. (15) into Eq. (10) and making use of the results’ 


x 
1 
f Jo(cp) cos pydp = — <a ¥ Si 
/ g \ c* cased y? 
e c 
pJ (cp) cos pydp = rel, os 6 4 
0 te = un)? . 
we find that the normal component of the displacement along the crack is given by w, 
where 
2(1 — o?)poc & V(gn + 3) f ¢ ee ae 
PO es KB. Ae SO! te 
Vr E  j-0 T($n + 2) W/e? — y? c 1 (w? — 1)8”? 


6 |. W. Busbridge, Proc. London Math. Soc. (2) 44, 115 (1938). 
7G. N. Watson, The theory of Bessel functions, 2nd ed., Cambridge, 1944, p. 405. 
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For the case of a uniform pressure po we take d9=1, a,=0, n21 and find 


w 


ai — «* eee 
= > = c= y¥*. (17) 


If we write 


b = 2(1 — oa”) poc/E, 


Eq. (17) reduces to the form 


which shows that the effect of the uniform pressure is to widen the crevice into an 
elliptic crack. 

4. It is also of interest to determine what distribution of pressure will produce a 
crack of prescribed shape. In this case we assume that the value of the normal dis- 
placement u, is known all along the y-axis; we have 


w(y), ys | cl, x = 0, 
wml | 
0 


IV 


v2lel, x=0. 


’ 


Inverting Eq. (10) by the Fourier cosine rule and substituting this value for u, we 
have 
_ E c 
o(p) = — — —_— w(y) cos pydy. (18) 
2(1 = a~) 0 





With this value of 6(p) in Eqs. (4), (5) and (6) we obtain expressions for the compo- 
nents of stress in the interior of the elastic solid. 


; - 
wy) = ef 1 — —], 
c 


- ve sin cp 
$() = — ——_—— (= — cos or). (19) 
(1 — o*)cp\ op 


For example if we take 


then, from Eq. (18) 


Substituting from (19) into Eq. (9) we obtain for the normal component of the stress 


along x =0, 





yu 
sin # sin 
2Ee i i c 
Po ewe lee - du }. (20) 
m(1 — o°)e cJo u 
Now 
* cos gx — cos px Pr 
f —————amevome G4 = § log — 
0 x i 
so that Eq. (20) reduces to 
2Ee y CTF 
o, = — ————_| 1 — — log ——-|,  0<y<e, (21) 
(1 — o*)c 2c c-y 
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giving the normal component of stress along the crack. This stress is negative when 
y=0 but becomes positive for a value of y between 0 and c, so that if a crack of this 
shape is to be maintained the applied stress must be tensile (and very large) near the 
edges y= +c of the crack. 

5. Expressions for the potential functions w(z), Q(z) of Stevenson corresponding 
to this problem can easily be deduced from the analysis of Section 3. It was shown by 
Stevenson,’ that if we write 

0 =o.+ ay; © = g, — Oy + 2itrey, D=uzt Uy, 


then the components of the stress and the displacement can be expressed in terms of 
two “potential” functions w(z), Q(z) by means of the equations 
l+o ze ) 
D = —— E{(3 — 40)Qz) — 20/(2) — a’(2)} | 
4 
20 = Q(z) + 2'(2) | 
e : | 
— 26 = 20’"(z) + @’’(2) 
in the absence of body forces. 
It follows from Eqs. (4) to (8) that the stresses and the components of the dis- 
placement vector may be derived from the potential functions 


4 6° Go) 4 (pe) 
Q(z) = -— —— ¢ dp, w'(z) = —— (1 + pz)e~?*dp, (23) 
us 0 p Tv 0 p 


where ¢(p) is given by Eq. (15) in the case where the applied internal pressure is 
given by Eq. (14). 

6. We now consider the distribution of stress in the solid when the crevice 
—cSy Sc, x=0 is opened up by the action of a uniform pressure po. Taking ao=1, 
a,=0, n>0, in Eq. (15) we obtain for ¢(p) the expression 


1 ¢p 
G(p) = — Erpoc*p otc) if f sTi(s)dsb . 
Cp’ 6 ) 
Now, 
| 27J\(s)dz = c*p*J2(cp) 


and, by a well-known recurrence relation, 


? 
J o(cp) ++. J (cp) = : J (cp) 
cp 
so that 
d(p) = — brpocJ (cp). (24) 


Substituting from Eq. (24) into (4), (5) and (6) we obtain the equations 
2) 


A(o. + oy) = — poc f e~°* cos pyJi(cp)dp, (25) 


* A. C. Stevenson, Phil. May., (7) 34, 766 (1943). 
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oo 


(Oy — oz) = pocx f pe’? cos pyJ(cp)dp, (26) 
0 


- poce f pe’? sin pyJi(cp)dp (27) 
0 


II 


for the determination of the components of stress. 


Now, 


oo 
J pe~"*I (cp)dp = c(c? + 2*)*”, 
0 
so that writing 
s= x+ iy = re, 3s — ic = rye", ot+ ic = reei® (28) 


we obtain the formula 


oo 


C 
f Ji(cp)pe-°?(cos py — i sin py)\dp = ——— ¢~ 8/2 (+62) | (29) 
/ 0 (rir2)3/? 
In a similar way we can establish that 
7 sats 1 eee 
J (cp)e-°*(cos py — isin py)dp = — 41 — ——— ei(@-ii-it2) | (30) 
0 c \ (r\%2)!/? f 


Equating real and imaginary parts in Eqs. (29) and (30) and substituting into (25), 
(26), and (27) we obtain the expressions 


1 = j . cate — i] 
(02 + oy) Po ) wOS (0 29 242) lf 4 


\ (rire)! 
rcon 8 / oc \e" 3 
(oy —¢o,)= Po ( ) COS - (A; + M2), i (31) 
c rife 2 
r cos 6 Cc 3/2 ; 3 | 
Tzy = — po- : sin — (0; + @2) 
G rile 2 


for the components of stress. Equations (31) are agreement with those derived in the 
previous paper;’ in making the comparison it should be noted that the angles 6, 61, 42 
of this note are the complements of the angles denoted by these symbols in the paper 
quoted. 

It follows from Eqs. (23) that these equations are a consequence of the Stevenson 
equations (22) if we write 


Q(s) = 2po[Vc? + 2 — s], w'(s) = — 2poc?(e? + 2?)-4. 





ANALYSIS OF SHEAR LAG IN BOX BEAMS BY THE 
PRINCIPLE OF MINIMUM POTENTIAL ENERGY* 


BY 
ERIC REISSNER 


Massachusetts Institute of Technology 


1. Introduction. Let us consider a thin-walled box beam of web height 2h and 
cover sheet width 2w which is bent in such a way that one of the cover sheets is in 
tension while the opposite cover sheet is in compression (Fig. 1). In elementary beam 
theory the assumption is made that the normal stress in the cover sheets does not 
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Fic. 1. Sketch of spanwise element of box beam with doubly symmetric cross section. 


vary in the direction across the sheet. Because of the shear deformability of the cover 
sheets this assumption of elementary beam theory is often seriously in error for wide 
beams. In aeronautical engineering this effect is known under the name of shear lag. 

In recent papers,!:? shear lag in box beams has been analyzed by an application 


* Received Feb. 22, 1946. 
1 E. Reissner, Least work solutions of shear lag problems, Journal of the Aeronautical Sciences, 8, 284- 


291 (1941). 
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of the theorem of least work which is the basic minimum principle for the stresses. 
The present paper contains an application to the problem of shear lag of the theorem 
of minimum potential energy, which is the basic minimum principle for the strains.* 
It is shown that application of the theorem of minimum potential energy to the pres- 
ent problem leads to simpler and more general results than the application of the 
theorem of least work. While the least-work method furnishes the stresses in box 
beams with no cut-outs, application of the minimum-potential-energy method fur- 
nishes, in a simpler manner, the stresses in beams without or with cut-outs. It also 
furnishes beam deflections, and is equally convenient for beams supported in stati- 
cally determinate or in statically indeterminate manner. 

Application, in the manner described below, of the minimum-potential-energy 
principle to the problem of bending of thin-walled box beams leads to a differential 
equation for the beam deflection which is a generalization of the relation 2” = — M/EI; 
this differential equation contains an additional term proportional to the fourth de- 
rivative of z which takes into account the shear deformability of the cover sheets. 
As the order of the differential equation in this theory is higher than the order of the 
differential equation of elementary beam theory, boundary conditions appear in addi- 
tion to those of elementary beam theory. These additional boundary conditions are 
different for beAms with cut outs and for beams without cut outs. 

The manner of application of the results obtained in the present paper is shown by 
solving explicitly the following four examples. 

1. Simply supported beam. Load distributed according to a cosine law. 

2. Cantilever beam with uniform load distribution. Cover sheets fixed at the sup- 
port. 

3. Cantilever beam with uniform load distribution. Cover sheets not fixed at the 
support. 

4. Beam with both ends built in. Uniform load distribution. 

For the sake of simplicity, it is assumed in what follows that the cross sections 
of the beams are rectangular and doubly symmetrical. It also is assumed that there 
is no continuous variation of cross-sectional properties. 

The author believes that the way in which the principle of minimum potential 
energy is here applied to the problem of shear lag will prove useful in other problems 
of structural mechanics. As an example of such future application, the theory for 
combined torsion and bending of beams with open or closed cross sections is men- 
tioned. 

2. Formulation and solution of problem. In the following, we analyze a box beam 
of doubly symmetrical rectangular cross section, composed of cover sheets, sidewebs 
and flanges. A given distribution of loads is applied to the sidewebs, acting normal 
to the plane of the cover sheets (Fig. 1). To this load distribution there corresponds a 
distribution of bending moments M(x). The spanwise coordinate being x, let y be the 
coordinate in the plane of the cover sheets perpendicular to the x direction, and 2(x) 
the deflection of the neutral axis of the beam. 





2 F. B. Hildebrand and E. Reissner, Least work analysis of the problem of shear lag in box beams, 
N.A.C.A. Technical Note No. 893 (1943). 

3 For a formulation of these theorems see for instance |. S. Sokolnikoff and R. D. Specht, Mathe- 
matical theory of elasticity, McGraw-Hill Book Co., Inc., New York, 1946, pp. 275-287. 
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The potential energy of the bent beam may be considered as composed of three 
parts. The first part 1s the potential energy of the load system. This may be written 
in the form 


d*z 
II; = f mx) — dx, (1) 
dx? 


the integral being extended over the entire length of the beam.‘ The second part is 
the strain energy of sidewebs and flanges. This may be written in the form 


> 


1 d*z \? 
ll, = - f EI, ( “) dx, (2) 
2 dx? 


the quantity J, denoting the principal moment of inertia of the two sidewebs and 
flanges. 

The third part is the strain energy of the two cover sheets. If it is assumed that 
the normal strains in the chordwise direction in the sheets may be neglected, as dis- 
cussed in the reference given in Footnote 1, then the strain energy of the two sheets 


is given by the integral 
1 2 2 
ll, = 3 ff 2t| Ee. + Gy |dxdy, (3) 


where the quantity ¢ denotes the cover sheet thickness, and where E and G are the 
effective moduli of elasticity and rigidity. Spanwise normal strain ¢€, and shear strain 
y are then expressed in terms of the spanwise sheet displacement u as follows 


Q 


Ou u 
= ’ y,=—: (4) 


€o = 
Ox y 





| 


Q 


The theorem of minimum potential energy states that the total potential energy 


il = Il, + Wl, + Wy (5) 


becomes a minimum for the correct displacement functions u and 2z, if only such dis- 
placement functions are compared which satisfy all conditions of support and con- 
tinuity imposed on the displacements. 

Direct application of this condition by means of the calculus of variations leads 
to a partial differential equation for u and to a complete system of boundary condi- 
tions. In what follows, an ordinary differential equation for the beam deflection z 
and boundary conditions for it are obtained instead. This is done by making a suitable 
approximation for the sheet displacements u and by applying the rules of the calcula- 
of variations to the resultant approximate expression for the potential energy func- 
tion. 

A reasonable assumption for the spanwise sheet displacements is 


dz y” 
u(x,y)= th E “- (1 = *) via) | (6) 
dx w? 


4 Eq. (1) implies that the beam is supported in such a manner that the end forces and moments can 
do no work. This restriction shortens the developments slightly. 
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The second term on the right of Eq. (6) represents the correction due to shear lag. 
Instead of the vanishing chordwise variation of the sheet displacements of elementary 
beam theory, we now assume a parabolic variation. The relative magnitude of the 
function U is a measure for the magnitude of the shear lag effect. The form of the 
correction is such that continuity of the displacements along the flanges, that is along 
y= +w, is preserved. 

Denoting differentiation with respect to x by primes, we obtain the following ex- 
pressions for the strains in the sheets from Eqs. (6) and (4): 


y? 
ec= th & (1 a ~)u’|, (7) 
w? 


_ wae, 

eo eT (8) 
WwW WwW 

On the basis of Eqs. (7) and (8) the following expression for the strain energy of the 

sheets is obtained: 


f y? 2 2 y . 
I, = ff all B| 2" + (1 -- "| +G | “- v] \ ava (9) 


In Eq. (9) the integration with respect to y is carried out. Setting 


I, = 4wth’, J =I, +T., (10) 
we have 
1 j :. 4 ; G4. 
i, = fe. (2’’)? + —(U')? + — 2”"U’ + — — U*} dx. (11) 
2 l 15 3 E 3w’ 


Substituting Eqs. (11), (12) and (1) into Eq. (5), we obtain the following expression 


for the potential energy of the system 


¢ 
I = f rs EI(s"’)? + Mzbax 
(2 J 
1 Li > FOO ne 
+f. Ely (U)? +4 —2"U' + ude. (12) 
2 15 3 E 3w? § 
Differential equations and boundary conditions for z and U are obtained by 
making 


él = 0. (13) 


Thus, with x; and x. denoting the ends of the interval of integration, 


2 
sl = f ‘ EZ +845 BLU | bs!” 
2 - 4 


I 
f 8 2 i] 2 
+ <FEI, E U'+- "| ou = (). (14) 
» ) 
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As 62” and 6U are arbitrary in the interior of the interval (x, x2) the terms multiply- 
ing them must vanish. This gives the following two differential equations 


‘ae © M 
e’/+— — U’'+— = 9, (15) 
ce EI 
5 GU 5 s 
EI, Ge _ —-+- "| az (), (16) 
2 E w? 4 


The integrated portion of Eq. (14) defines the boundary and transition conditions 
for the function U’. At a section where the sheet is fixed, 6U=0 and 


U =0. (17) 
At a section where the sheet is not fixed and consequently 6U is arbitrary, 
EI,[U’ + $2’’] = 0. (18) 
Transitions conditions for adjacent bays with different stiffness are: 


U and EJ,[U’ + 42’’| continuous. (19) 


co 


The above boundary and transition conditions are in addition to those imposed on z 
and M in elementary beam theory, as may be verified by repeated integration by 
parts of the term containing 62’’ in the integral of Eq. (14). 

3. The modified beam equation and its boundary conditions. By climinating the 
quantity U from Eqs. (15) to (19), we obtain a system of relations containing the 


beam deflection z only. 
The differential equation for z is derived by differentiating Eq. (16) and substitut- 


ing LU’ from Eq. (15). There follows 


M te ae i eee 
o’ + — — w’- | (<" + ) aa | = 0. (20) 
El GLS EI 3] 


When the shear deformability of the sheets is neglected, that is when it is assumed 
that G= «, Eq. (20) reduces to the well known result of elementary beam theory. 


Equation (20) may be written in the alternate form 


2 E 5 I.\ 2!V M 2 E M” 
“=a ee : (21) 
5 G w? EI 3 G wl 


With the help of Eqs. (15) and (16), the boundary condition (17), which holds 
when the sheet is attached to the support, is transformed into 


D- dy M’ 
(1 —— =)" a ae, (22) 
6 71 EI 


Similarly, the boundary condition (18), which holds when the sheet is not attached 


5 i, M 
(1 -— “) s’ +— = 0. (23) 
ae EI 


to the support, becomes 
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The continuity conditions (19) may be transformed in an analogous manner. 
The values of the sheet stresses may be obtained from Eqs. (9) and (10). From 
Eq. (9) it follows that the flange stress is given by 


o, = + Ehs". (24) 


For the application of the results it may be noted that the differential equation 


(21) can first be solved for the value of 














s’’ which, according to (24), gives di- y T 
rectly the approximate value of the = bs 
flange stress o,. The magnitude of the 

deflection z can then be found from the ~ oc aemesqnet 


value of 2”’ as in elementary beam the- aT 





| 
Abbi bbb hibit tdi ts 
te 



























































ory. 
For the evaluation of the solution f 
we define the following two parameters (1) _ 
2 
7S ee eae Fs (2 5) 
1 — 5/,/61 
——-—-- tT 
F 1 Y foes G > 
7 = - — (26) 
w 4 2 E = x = (2) 
“ 
With (25) and (26) the differential equa- ‘ 
tion (21) becomes 
ae. eS ode ITT it iii iid | 
Ce ae a ree po —— 5 (TF) ae 
k? EI k? El . 
the boundary condition at an end sec- —_ 1 T 
tion where the sheet is attached to the sy 72> 
support becomes b= x bs (3) 
M’ | 
J’ = —n—; (28) i 
EI 
— 1 — 
and the boundary condition at an end 
section where the sheet is not attached 
to the support becomes 
M y 
3’ = — n— > (29) t=5 f . = 
EI 2w — or 
. " 
4. Examples of applications (Fig. _ Dic tiniertinniniannnaaiiolill 
2). 1. Simply supported beam. Load dis- 
: ; : : ooo oot y. 
tributed according to a cosine law. Desige Petiiiiiiiiiiiiiiiiiis 
nating the span length of the beam by / 7 4 








and assuming the origin of the coordi- (4) 
nate system at the center of the beam, Fy. 2. Diagrammatic sketches of beams analyzed 


we consider the moment distribution as examples of application of the theory. 
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tN 


x 
M = Mocosr r . (30) 
A particular solution of Eq. (27) is 
L\? Mo 1+ n(x/ kl)? x 
EI 1-+ (x/kl)? 


T 


As Eq. (31) satisfies the boundary condition (29) and the condition of vanishing de- 
flection at the ends of the beam, it is the complete expression for the deflection func- 
tion. When 1/k=0, Eq. (31) reduces to the expression for z in the case where shear 
lag is not taken into account. The factor 

1 + n(x/ kl)? 1 + (29°E/5G)(w/l)? 

1 + (x/kl)? 1 + (29°E/5G)(w/l)?(1 — 5/,/67) 
expresses the effect of shear lag on deflection and flange stresses. 

2. Cantilever beam with uniform load distribution. Cover sheets fixed at support. As- 

suming that, contrary to what is indicated in Fig. 2, the free end of the beam has 
the coordinate x =0 and the fixed end of the beam the coordinate x =/, we may write 


the moment distribution in the form 


x 2 
M = Mo ( ). & 
l 


The differential equation (27) then becomes 


1 M, x\? 2n 
\ = a | (- ) aulencieean | (34) 
k? EI l (ki)? 


we 
w 
~~ 


Solving for 2’’, we find 


" Mf. bie a > hk (*)+ 2(n — 1)) (35 
5 = — sinh kx 2cosh Rx — [- ———a , 35) 
EI \"' i (ki)? J 
Satisfying the boundary condition (29) when x =0 and (28) when x=/, we obtain 
Mof{/x\*? 2(n — 1) sinh kl — kl . \ 
f= — — ( ) + —--- - (cosh kx — 1) — ——————— sinh kx |? . (36) 
EI \\1 (Rl)? cosh kl § 
According to Eq. (24), the flange stress at the fixed end of the beam becomes 
M oh 2(n — 1) 1 1 q 
o/(l) = = ——<1+ ———— | tanh &/ — — + ————— ; (37) 
i | kl kl kl cosh kl 
We take for a numerical example 
ie 1 G 3 l 5 
ka ea bac ed cael Wi tee (38) 
I 2 E 8 2w 2 


so that according to Eqs, (25) and (26) 
n = 1.714, ki = 6.34, (39) 
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and we find 


Moh 
ofl) = F — {1 + .190}. (40) 


By application of the least work method! a factor 1.186 is obtained instead of the 
factor 1.190 in Eq. (40). 

The deflection of the beam is obtained from Eq. (36) by integrating twice and 
making 2(/) =z'(1) =0. In the present case, the correction due to shear lag for the 
maximum deflection is about ten percent. 

3. Cantilever beam with uniform load distribution. Cover sheets not fixed at support. 
Moment distribution and differential equation are given by Eqs. (33) and (34). The 
constants of integration in (36) are determined by satisfying Eq. (29) for x =0 and 


for x =l. There follows 


M> f z\? 2(n — 1) 
= — — : + -——-.--- (cosh kx — 1) 
EI \\1 (Rl)? 
cosh kl — 1 — 3(kl)? 


sin <cneecinnnshewacuiiememes gills Ext. (41) 


sinh kl 


Taking again I,/I =.5, we should have, for the flange stress at the supported end, a 
value twice as large as the stress according to elementary beam theory for a beam 
with sheet attached to the support. In the present solution the factor 2 is replaced by 
n =1.714. This indicates that with the assumed parabolic chordwise variation of sheet 
displacement the condition that at the support of the beam the sheet is free of stress 
is only approximately satisfied. The same difficulty arises in methods which incorpor- 
rate the ability of the sheet to carry normal stresses as effective width contributions 
to the strength of stiffners.® This difficulty is not serious when the main purpose of 
such “cut-out” calculations is the determination of the distance over which the cut- 
out is effective and its effect on the over all beam stiffness.® 

The localization of the effect of the cut-out may be seen by writing (41) in the form 


Mo» ie oh" \ 
o" x - 4(=) + @- pene os (42) 


This equation indicates that the influence of the cut-out is small as soon as the dis- 


tance /— x satisfies the inequality 


— log (n — 1) 2E 5 1, 6 1 
f= # > ————_- = w (1 _— “) log —— 1). (43) 
k 5G o # a P 


Thus, the wider the sheet and the smaller the value of the shear modulus G, the far- 
ther away does the effect of the cut-out extend in the spanwise direction. 
The magnitude of the beam deflection is obtained from (41) in the form 





5 P, Kuhn and P. Chiarito, Shear lag in box beams—methods of analysis and experimental investigations, 
N.A.C.A. Technical Report No. 739 (1943). 

6 Exact solutions of problems of this kind have been obtained by F. B. Hildebrand, The exact solu- 
tion of shear-lag problems in flat panels and box beams assumed rigid in the transverse direction, N.A.C.A., 
Technical Note No. 894 (1943). 
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Ax) = | f 2!"(x2)dxed x, (44) 
l l 


which determines the constants of integration such that 2(/) =2’(/) =0. For the de- 
flection at the free end of the beam, we have 


Mol? (1 An —1)F/1 1 1 1 y 
EI \4 kl 2° (ki)? kl (kl)? sinh RIJS 


For a beam with dimensions as in (38) and (39), Eq. (45) becomes 


Ml? | 
2(0) = —— (.25 + .083). (46) 
EI 
This indicates that for a beam with dimensions as given shear lag due to lack of sheet 
restraint at the supported end of the beam 1s responsible for a thirty percent increase 
of the maximum beam deflection as compared with the result of elementary beam 
theory for a beam fully restrained at the supported end. This increase of deflection 
of thirty percent compares with one of hundred per cent which is obtained if the con- 


tribution of the cover sheets is neglected. 
4. Beam with both ends built-in. Uniform load distribution. The distribution of 


bending moments may be written as 
x 2 . 
M = Mo 73 + M;,. (47) 


The value of Mo is determined by the load intensity, the value of MJ, in this statically 
indeterminate problem has to be determined from the displacement boundary condi- 
tions. The boundary conditions are 


l l 
:(+ =) = 0, (48) (+ 5) = 0, (49) 


( -) M’(+ 1/2) 
Fe (50) 
EI 





For these boundary conditions the moment distribution is not affected by shear lag, 
provided the moment distribution is symmetrical about the mid-span section of the 
beam. Indeed, the differential equation (27) may be integrated to give 


eC f M ln + n M’ 
f-—a — ae tae (51 
k > 2 ”_ 


9 


the limits of integration being so chosen that Eq. (51) satisfies the conditions of zero 
slope and zero vertical shear at the mid-span section. In view of (49) and (50), 


a ile 
— dx = 0, 52 
, EI (52) 


Eq. (51) implies 


1946] ANALYSIS OF SHEAR LAG IN BOX BEAMS 277 


regardless of whether or not shear lag is taken into account. A considerably less simple 
proof of the same fact by means of the least work method has been given in the refer- 
ence quoted in Footnote 2. For the moment distribution of Eq. (47) there follows, 
from (52), 
Mo M, 
24 2 


zy" if 
Uy 12 


With this value of M and the requirement that s’”’ be an even function of x, Eq. (27) 


= 0, (53) 


and hence 


is solved in the form 


Mo | su* 1 2(n — 1) 
i= — ‘ ( ) _ — ———| +C, cosh bx ; (55) 
EI \\1 12 (kl)? 


The constant Cz is determined from Eq. (50). There follows, 


Mo({(/x\? 1 n—1f[ cosh kx 1 , 
—G-Eidd « 
E l 12 kl Lsinh ki/2 = ki/2 


Taking a beam five times as long as wide, that is //2w=5, and assuming the re- 
maining parameters as in (38) and (39), we obtain the following expressions for the 
flange stresses at the built-in section and at the center section of the beam 


l Mok 1 | 
od — = + — —(1 + .283), (57) 
2 IT 6 
- Moh 1 F - 
oO) = + —_ Dp" + .106). (58) 


These results agree to within a fraction of a percent with the corresponding results 
obtained by the least work method.? It is worthy of note that, for this beam with both 
ends built-in, shear lag is considerably larger than for a cantilever beam with the 
same load, same width and half the span of the beam with both ends built-in. If both 
beams had the same span, the discrepancy would be even larger. 

The deflection z of the beam is obtained from (56) and (48) in the form 


Mf? jl x\* i/jz\* 1 
EI \12\1 24\ 1 192 
n—1f1 x\* cosh (1/2) cosh kx ) P 
+ ra techies “) ne IP (59) 
(kl)? L4 l kl sinh kl/2 s 
Corresponding to the stresses of Eqs. (57) and (58) we find for the deflection at mid- 
span 
1 MJ? 


s(0) = — — —— (1 + .145). 60 
' 192 El , (00) 
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Shear lag in this beam is thus responsible for an almost fifteen percent increase in 
deflection. This percentage increase of deflection, while appreciable, is considerably 
smaller than the percentage increase of maximum flange stress. 
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THE ANALOGY BETWEEN MULTIPLY-CONNECTED 
SLICES AND SLABS* 


BY 


RAYMOND D. MINDLIN 


Department of Civil Engineering, Columbia University 


1. Introduction. The analogy between the two-dimensional field of stress and the 
transverse flexure of a thin plate was first applied by K. Wieghardt! to the solution 
of a problem involving boundary loading of a simply-connected body. As is well 
known, the analogy establishes the proportionality of the curvatures of the surface 
of the plate to the components of stress in the two-dimensional field of stress. H. M. 
Westergaard? introduced the useful terminology of slab and slice, free slice and con- 
strained slice, and gave the boundary conditions for the slab when the slice is multi- 
ply-connected and is stressed by boundary loads having no resultant force on an 
internal boundary. Westergaard also proposed the use of the analogy in the investiga- 
tion of the stresses in the Boulder Canyon Dam,* a problem involving gravity and 
boundary loading of a simply connected body. An improvement in experimental 
technique was contributed by H. Cranz‘ in introducing an optical spherometer® for 
measuring the components of surface curvature. Cranz's application was to boundary 
load problems in simply connected bodies. 

[t is the purpose of this paper to give the general boundary conditions for the 
slab when the slice is multiply-connected and is stressed by any combination of 
boundary loading, body forces, dislocations and thermal dilatations. The analogy has, 
in fact, its most useful applications in the last three cases as they are either difficult 
to reproduce, or the resulting stresses are difficult to measure, in an experimental 
model of the slice itself, while the analogous conditions for the slab, developed below, 
are easy to handle. 

In order to proceed, it is necessary, first, to set down the general boundary value 
problem for the slice. It is convenient to do this along the lines established by 
Michell,® with the additional consideration of dislocations and thermal dilatations. 

2. Airy’s stress function and its differential equations. In a state of plane strain 
defined by setting 

* Received April 9, 1946. 

1K. Wieghardt, Uber ein neues Verfahren, verwickelte Spannungsverteilungen in elastischen Kérpern 
auf experimentellem Wege zu finden, Mitteilungen iiber Forschungsarbeiten a. d. Gebiete d. Ingenieur- 
wesens, 49, 15-30 (1908). 

2H. M. Westergaard, Graphostatics of stress functions, Transactions, Amer. Soc. Mech. Eng., 56, 
141-150 (1934). 

3 United States Department of the Interior, Bureau of Reclamation, Boulder Canyon Project, Final 
Reports (1938), Part V, Technical Investigations: Bulletin 2, Slab analogy experiments; Bulletin 4, Stress 
studies for Boulder Dam. 

*H. Cranz, Die experimentelle Bestimmung der Airyschen Spannungsfunktion mit Hilfe des Platten- 
gleichnisses, Ingenieur-Archiv, 10, 159-166 (1939). 

° E. Einsporn, Ebenhett, Zeitschrift fiir Instrumentenkunde, 57, 265-285 (1937). 

6 J. H. Michell, On the direct determination of stress in an elastic solid, with application to the theory of 
plates, Proc. London Math. Soc., 31, 100-124 (1899). 
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tu ™@ Vas PP te = 0 


and restricting the displacements u and v to be functions of x and y only, the relations 
between strain, displacement, stress and temperature in an isotropic elastic body are 


Ou 1 2 = 

e. = — = — [(1 — no. — (1 + i) oy] + (1 + nidarT, (2. 1a) 
Ox Ey 
ov 1 2 , 

€y = — = — [(1 — n)oy — (1 + ry)o2] + (1 + ridaiT, (2. 1b) 
Ov om 
ov Ou 2(1 + 7) 

Voy = — + — = ——— Ty. (2. 1c) 

Ox dy E, 


These are the relations for a constrained slice. The notations for stress, strain and dis- 
placement are the usual ones and &,, »; are Young’s Modulus and Poisson’s ratio 
for the material of the slice, a; is the coefficient of linear thermal expansion, and T is 
the temperature in excess of a uniform initial temperature. 

When the stresses are expressed in terms of Airy’s stress function (@) and a body 


force potential (V) by 


0° " 0° ? d*p 
op=—+V, o=—+V, ty=—-——: (2.2) 
dy? Ox? Oxdy 
the equations of equilibrium are satisfied and the strain relation 
0*€, Oey 07 zy 
—— + — = —— (2.3) 
Oy Ox? Oxdy 
yields the differential equation governing ¢: 
1 = 2; rows 1 + Vi aoe 
Vg = — —— VV — ——a,V°T. (2.4) 
— yy 1— vr 
In a state of plane stress, defined by 
wg ~*~ Vee ~ Tee 0, 
the strain-displacement-stress-temperature relations become 
Ou To = : 
ee a (¢, — vidy) + al, (2.5a) 
Ox E, 
Ov 1 = > 3 
€y = — = — (¢, — 10s) + al (2.35b) 
Oy FE, ; 
Ow Vy : ees o 
é->-— = oe ae (or + oy) + aT, (2.5c) 
Os E 
Ov Ou 2(1 + v1) pare 
Vey = —— + — = Try. (2. 5d) 
Ox ov EF, 


These are the relations for a free slice. If the components of stress are again ex- 
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pressed in terms of an Airy function and a body force potential by (2.2), the equilib- 
rium equations are identically satisfied and the strain relations reduce to 


V'g@ = — (1 — »)V?V — (1+ n)aV*T (2.6) 


if terms associated with the coordinate z are neglected. 
In what follows, the case of plane strain (constrained slice) will be treated, but 
the results are directly applicable to plane stress (neglecting z-dependent terms) if 
Young’s modulus £, Poisson’s ratio », and the linear thermal expansion coefficient a, 
are replaced by Ey, vi and aj where 
E\(1 + 2x) , V1 : ay(1 + ») 
1 + 21 


rea {2 , ay 

(1 + 1)? i+y 
3. Conditions on ¢ at a point on a boundary of the slice. Michell® gave the condi- 
tions to be satisfied, at each point of each boundary, by ¢ and its derivative normal 


(2.7) 


"7 


1 = 


to the boundary: 


¢ =f (Bl — Am)ds + ax + By + +, (3.1) 
do 
; = Al+ Bm + al + Bm, (3.2) 
an 


where a, 8, y are constants, in general different for each boundary, ds is an element 
of arc of a boundary, dn an element of normal to that boundary, and 


. dy ~- dx (3 3) 
Pe my m=—), J. 
ds ds 
[ = -f Yds + q Vmds, B -f Xds -f Vids, (3.4) 
0 e 0 0 0 
X =o, + rrym, Y = 72,1 + oym. (3.5) 


In a simply connected body, a, 8, y may be assigned arbitrary (including zero) 
values as the addition of a linear function of x and y to ¢ does not affect the stresses. 
In a multiply-connected body, three additional conditions on ¢ are required for de- 
termining a, B, y, on each additional boundary. Equations (3.1) to (3.5) are not al- 
tered by introducing thermal dilatations and dislocations of the type considered here. 

4. Conditions on ¢ for each boundary of the slice. The additional conditions on @ 
are obtained by assuming the strains (and hence the stresses) to be continuous and 
requiring the rotations and displacements (a) to be single-valued or (b) to have pre- 
scribed discontinuities (dislocations). Michell® gave the conditions for case (a). The 
conditions for case (b), including, also, thermal dilatations, are derived by following 
Michell’s procedure with modifications along the lines indicated by Volterra.’ 

(i) Rotation condition. Considering the rotation 


1 /dv Ou 
0, = = ~ =), (4.1) 
2\e0x dy 


7 Love, Theory of elasticity, 4th ed., Cambridge Univ. Press, Cambridge, 1927, pp. 221-228. 
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we require that the line integral of its differential have a value, say c, after one com- 


plete circuit around (and along) a boundary. Thus, 


c =f te, (4.2) 


Ow, Ow, 
f te = ¢ —dx + —dy 
Ox Oy 
1 oy, O€, de 1 Oy. 
aD een ee 
2 Ox oy Ox 2 oy 


Replacing the strain components by their expressions in terms of ¢, V and 7, we find 


E 16 0 0 OV OV 
= (2 = ng ( -(V*o)dy — - (v'6)ds) + (1 — 2y;) g ( -—dy- — i) 
1 + py; Ox oy Ox Oy 
oT oT 
+ Ea ¢ sae dy = 8 iz), 
Ox oy 
d(V*o) Fic 1— 2n, dV Ey\ay dT 
¢ ds = - - ¢ ds — ——— - ds, (4.4) 
dn 1 —p 1-7 dn 1-—y dn 


1 


Now, 


Then 


This is the first of Michell’s three conditions on ¢ for each boundary of the slice. It 
may be observed that, if the circuit of the line integral in (4.3) were reducible, the 


integral would vanish because, by Green’s theorem, 
1 OVcy de, dey 1 OY zy 
$1 22-29) co (*- yO 
2 ox Oy Ox 2 oy 
0*e, 0? ex O7Y ry 
mie — —— } dxdy; (4.5) 
Ox? a y? Oxdy 


and the surface integral vanishes by virtue of (2.3). 
(41) Displacement conditions. We admit a translational dislocation a parallel to x 


Ou Ou 1 
= g aw = ¢ dx + y= § (cas + veil) - § way. 
Ox Oy y 
g wdy = mh dw, — g¢ ydw, = yor — $> dw, 


where yo is the y-coordinate of the starting point of integration. Also 


OW, 
g ydw, -$ (= dx “0 “ ay) 
, Oy 
1 Oy. Oe, de 1 oy: 
-$ (> id te ) as + hy (- ee *) dy. 
"2 Oe dy Ox 2 ody : 


and set 


Now 
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Hence 


7 1 Oy. de, 
a+ Voc = g le. + (> —— =) Jas 
2 ox Oy 
1 de 1 Oy: 
+$| Yay + (° -——- 2) Nay. (4.6) 
2 Ox 2 oy 


We now note that 
0€; 1 OY cy 


F 1 0 1 0 
dx + —yeydy = [xez]o + — [yr] — @ (x— dx + —y— dy 
¢ € dn ; VY ry [ver Jo + r [vy cy le g (: — dx + r J ay iy) 
0€, 1 OY zy 
= - $ (« — dx + — y—_ dy}, 
Ox 2 dy 


the terms outside the integrals vanishing because of the assumption of continuous 


strains. Equation (4.6) then becomes 


1 Oy: de, Oe, de OY: 
a+ yoc = ri | »( -— ) —x— | dx +g v( _— =) dy. (4.7) 
2 ox Oy Ox Ox Oy 


When the strain components in (4.7) are replaced by their expressions in terms 


of ¢, V and 7, we find 


Ey(a + voc) 0 0 
=(1- ng »| (V*o)dyv — (v6)dx | 
1+» ax ay 
ov aV 
+ (1 — 2) b 9(— dy —- — ix) 
Ox Oy 
oT oT 
t Pash v( dy - — ix) 
Ox Oy 
aT OV oT 
-$ | a — vy) — (V*b) + (1 — 2x1) — + Eya, — |dx 
Ox Ox Ox 
d%q dp 
+ ¢ (« dx + y—— iy). (4.8) 
0x8 xd y* 


a 0% 0*p 0° 0° 0 0° 0" 
$ (: dx+y iy) = E + y | - $ ( - dx + — iy) 
dx? ~ Oxdv* Ox? OxO Vo Ox? Oxdy 
d {dd 
a -$ ( ) as, 
ds\0x 


the term outside the integral vanishing because the stresses are continuous. But, 


l (do —_ 
; ( )- Vm — Y. 
ds\Ox 


from (2.2), (3.3) and (3.5), 
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Hence (4.8) may be written: 
d(V*o) d(V") Ex(a + yoc) 1— 2m dV dV 
g¢ y— —- x - ) a = — — y — — x — fds 
dn ds i « 1— dn ds 
Eyay dT dT 
= g (» -— 2—lHMs 
1- » dn ds 
1 “ 
—_- g¢ (Y — Vm)ds. (4.9) 
1 = 


This is Michell’s second condition on @ for each boundary of the slice. 
Similarly, admitting a translational dislocation b in the y-component of displace- 


ment, we set 


and we find 


d(V°o) d(V*) E\(b— xoc) 1 — 2% dV dV 
$(: — + «— *) ds = — =a Decent ad _ —$ ¥ + x )as 
ds dn 1 — p;? 1-— » ds dn 
Eye, dT dT 
a $ (9 _ + x —) ds 
1- » ; a8 dn 
1 = 
_-—— gi X — Vijds, (4.10) 
1 = 


which is the last of Michell’s three conditions. 

Corresponding to (4.5), a similar application of Green’s theorem to (4.6) reveals 
that the right hand side of the latter vanishes for reducible circuits and the same 
result is found for the corresponding step in the development of Michell’s third con- 
dition. 

The differential equation (2.6), the boundary conditions (3.1) and (3.2), and the 
three conditions (4.4), (4.9) and (4.10) constitute a statement of the boundary value 
problem of plane elasticity for stresses induced by boundary loading, body forces, 
dislocations, and thermal dilatations. The general formulation of the problem reveals 
the analogies, discovered by M. A. Biot,® between gravity loading and boundary 
pressures, and between thermal loading and boundary pressures and dislocations. 

5. The slab equations. In the approximate theory of the bending of thin plates® 
(slabs), the deflection (w) is governed by the differential equation 


DV‘w = Z, (5.1) 


where D is the flexural rigidity of the plate and Z is the surface load, normal to the 


middle plane. 
The components of curvature in the y, z and x, z planes are given by 


5M. A. Biot, Distributed gravity and temperature loading in two-dimensional elasticity replaced by 
boundary pressures and dislocations, ]. Appl. Mech., 2, A 41—-A 95 (1935). 
* Love, loc. cit., p. 487. 
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Ow Ow 
Kk; = ae ’ Ky = - 
dy? Ox? 





(5.2) 


On a boundary of the slab, the shearing force (NV) normal to the middle plane, the 
flexural couple (G), and the torsional couple (77) (all per unit of arc length s) are 


0 
N = — D — (V*w), (5.3a) 
on 
0° w aw 1 dw = 
G = - o| = + 7 ("+——)|. (5.3b) 
on* as? p’ on 
0 fdw 4 
H = (1 — ».)D = (=), (5.3c) 
on \ Os 


where p’ is the radius of curvature of the boundary of the unflexed slab and ¥ is 
Poisson’s ratio for the slab material. 


The resultant force and the components, parallel to the x and y axes, of the re- 


sultant couple on a complete boundary are’ 


OH » 
F,= g(a -- —) ds, (5.4a) 
Os 
0H dx s 
M a 1 EG —_ -) + G =| ds, (5.4b) 
Os ds 
dy OH a 
M,= gle —-—-2 (v _ —)| ds. (5.4c) 
ds Os 


Substituting (5.3) in (5.4) we find 


0 8 8 Sow 
ay beens ees 
on Os On \ Os 
0 0 0 fdw 
M,=- vg $9]. (V2w) 4+: (1 — ve) — — (=) 
on ds On \ OS 
dx [0°w Ow 1 ow \ kone 
+ | + ve (— + a ds, (5.5b) 
ds Lon? Os* p’ on f 
jay 0*w Ow 1 dw 
M,= -p g- 3 +n (—+— — 
lds Lan? ds? p’ On 


0 0 0 fdw 1 a 
— 2x | (V2w) + (1 — ve) — —(— ¢ ds. (5. 5c) 
on Os On \ Os 


6. The analogy for singly-connected bodies. Noting the similarity between the 
differential equations (2.6) and (5.1) for @ and w, we set 


w= K®@, (6.1) 


5.5a) 


"ry 

ees 
I 
si 


10 Love, loc. cit., p. 460. 
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where K is a conversion constant having the dimensions of length/force. 
Then, from (6.1) and (2.6), 


1 + Vi e 
Vw = — ———_ KV*V — —— Ka,V°T (6.2) 
1 Pe 1 -— 


becomes the differential equation for the deflection of the analogous slab. Hence 


1 =— 24 1 + Vi _ 
La = ——— ZpvV — ———- EDaiV*1 (6.3) 


—= Fi 1-—v 


is the normal surface loading to be applied to the face of the slab. In the case of a 
steady state temperature distribution, 


wT = 0. (6.4) 


If, in addition, the body force potential is harmonic, the slab is subjected to edge load- 
ing only. If either V or T is not harmonic, transverse loading is required on the sur- 
face of the slab, and the load may vary slowly with time. 

The edge conditions (i.e., the elevation and slope at each point of a boundary) of 
the slab are specified by substituting w= K@ in (3.1) and (3.2). Thus 


. = f (BL — Am)ds + ax + By + ¥, (6.5) 
0 
1 dw 
-—— = Al+ Bm + al+ Bm. (6.6) 
K dn 


The normal components of stress in the slice are obtained by combining (2.2), 
(5.2) and (6.1), with the result 


0; = -+ /), t., = —-+ V. (6.7) 


The principal stresses and their directions may be calculated from two sets of curva- 
ture measurements at each point.‘ If the boundary of the slab is a scale model of the 
boundary of the slice, e.g., if the ratio of a linear dimension of the slab to the corre- 
sponding linear dimension of the slice is k, the stress components in the slice are given 
by 
ks k°Ky 
6, = + V, o, = + V. (6.8) 
K K 

For a singly-connected body, (6.1) to (6.8) completely specify the analogy, since 
the unknown constants a, 8, y may be given arbitrary values. 

7. Additional conditions on the slab for multiply-connected bodies. For a multi- 
ply-connected body, a, 8, y must be prescribed for each boundary. Now, it will be 
observed, from (6.5) and (6.6), that a, B, y specify a rigid body translation and rota- 
tion of each complete boundary of the slab. Such rigid body movements may be 
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effected by applying, on each boundary, a resultant force, normal to the middle plane 
of the slab, and a couple about an axis properly oriented in the plane of the slab. The 
magnitudes of the force and the x and y components of the couple on each boundary 
are determined by expressing F,, M,, M, (see (5.5)) in terms of the specified boundary 
loadings, body forces, dislocations, and temperature distribution of the slice. 

i. Resultant force on a boundary of the slab. Replacing w by K@ in (5.5a), we have 


0 0 0 fd¢ 
F,=- KD g | (V°d) + (1 — ve) —— (“)] ds. (7.1) 
on Os On \Os 


Now 
0 ¢0 fd¢ 
ft hws 
Os On \Os 
because of the assumption of continuity of the components of stress in the slice. Hence 
a = 
F, = — KD — (V*o)ds. (7.2) 
on 
Therefore, from (4.4), 
(1 — »,)F, Exc dV ; dT 3 
= — ——— + (1 — 2) —ds + Eja — ds, (7.3) 
KD 1+ 7 dn dn 


whereby F, is expressed in terms of known quantities. 
ii. x-component of couple on a boundary of the slab. Substituting K@ for w in (5.5b): 


M KD { E (V°6) + (1 — 2) ae (~)| 
a ¢ 4 Fhe "as an \as 
ix Ta? a? 1a 
i= Ee in (= ims 2) | as. (7.4) 
ds Ldn? Os? p’ dn/ J) 


0 
y — (V°o)ds 
on 
between (7.4) and (4.9), we find 


M, d(V*o) 0 0 do 
= -${ x + (1 — m)y— — - 
KD ds Os On \Os 


dx [a* 0° 1 d¢ 
+ [4 (S4— S) [has 
ds | an? ds? p’ on/ }) 


E\(a+ yor) 1-2 ( dV “~) 
net aan ae y — — x—]) ds 
~ dn 


1 — vi 1-— » ds 


Ejay; dT dT 1 - , a 
eis — gf (» — — ¢—} ds + ——_ (Y —Vm)ds. (7.5) 
1- » “dn ds 1-— » 


Eliminating 
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O(V"o) 0 0 dd 0 0d £0bO\ Wo 
g¢ E rt = wy ( )| ds = |x¥°o|o + (1 — v2) |» ( )| 
Os Os On \Os On \Os/ Jo 
dx dy 0 /d¢ 
- ¢ | V*bd + (1 — we) — — ( )| ds. (7.6) 
ds ds On \Os 


The terms outside the integral vanish on account of the assumption of continuity of 
the stress components. Therefore the first integral on the right hand side of (7.5) be- 


comes 
‘ dy 0 fd dx 0° 0* 1 0¢ 
Pee). ee oe) 
\ ds on \Os ds On- Os” p’ On f 


On a boundary 


~I 


: 0° 1 0p 0° 
7 oe tone 
on- p on Os” 


so that (7.7) becomes 
dy 0 {0d dx {0° 1 do 

(1 — po) - —j}-+- + — ds. (7.9) 
ds dn.\Os ds \ds* p’ on 


However, along a boundary, 


~ 


0 /d¢ 0* 1 d¢ 
(*) - — Tass (7.10) =f a ee, — F, (7.11) 
dn \ds as? p’ On 
dy dx 
Fag mm ge (7.12) 
ds ds 


Hence, (7.9) becomes 


ce dx 
—(1- ng (7 +V *) ds. 
ds 


Substituting back in (7.5), we have, finally, 


(1—»)M, E\(a + yoc) dV dV 
1+ y dn ds 


KD 
dT dT 
ian g (y= _ -) ds 
dn ds 
— [(1 — ve)(1 — v1) + lg @ — Vm)ds. (7.13) 


This gives M, in terms of known quantities. 
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iii. y-component of couple on a boundary of the slab. Substituting K@ for win (5.5c), 


iy [a a 14 
M, = — KD g¢ {2 = + (= +— “| 
ds Ldn? ds? p’ an 


0(V"°o) 0 0 fA 
~S eu- at Se 7 
on Os On \As/ J) 


$= (veyd 
re ds 


between (7.14) and (4.10), we have 


M, d a a (a 
iif innit oO 
KD ds Os On \Os 
ly [oO 07 1 0a 
a2 [ ae ( ae ey [has 
ds Ldn? Os? p’ on 
E\(b — xc) 1— 2 1V 1V 
= Pan hes. sein “ei te ¢ (> er + x ~) ds 
l-—», 1— ds dn 


Eya dT dT 1 
des —$ (» -+ 4 <- ds — —— gx — Vids. (7.15) 
1-—», ds dn 1-y 
Now, 


d(V°¢) a 9 (a¢ 0 [ a (~)] 
— (1 — wm)x — —[— is = |yV2 — (1 — ») |} « 4— 
¢ > ds \ ne Os On (~)] ” Ly ]o ( tot On\Os/ Jo 
ly ix 0 fd 
— g¢ E V°o — (1 — ve) aes (“)] ds. (7.16) 
ds ds on \Os 


The terms outside the integrals in (7.16) vanish on account of the assumption of con- 
tinuity of stresses. Therefore the first integral on the right hand side of (7.15) becomes 


‘dy 07 0" 1 da ix 0 (d¢ re 
f} ) [vs a (= **)) ee. < (=) rm (7.17) 
lds on? ds? p’ on ds dn\ds 


Then, using (7.8), (7.10) and (7.11) and noting that 


Eliminating 








dx — 
ese + os al XxX, (7 : 18) 
ds ds 


(7.17) may be written in the form 


(1 — rg (X¥ — Vids. (7.19) 


Substituting back in (7.15), we have 








290 RAYMOND D. MINDLIN 


(1—»)M, E,(b — xoc) dV dV 
—— = — —————_ - (1 - 20) § (» —+2 ~) ds 
KD 1+ » ds dn 
dT dT 
— Ejay $(> —+ 2x -) ds 
ds : dn 
+ [(1 — »2)(1 — 1) — ga — V1)ds. (7.20) 


8. Recapitulation. The stresses, in a multiply-connected slice, resulting from 


boundary loadings X, Y, a body force potential V, dislocations a, 6, c and tempera- 
tures 7, are related to the curvatures of a slab according to 


if the following conditions are satisfied on the slab: 
(i) The surface loading on the slab is 
1 = 2r1 F _ 1 +- Vi 2 ai 
Z = — —— KDV’V — —— KDa,V’T; (6.3) 


— 1-— py 


(ii) The boundaries of the slab are geometrically identical with those of the slice, 
with elevations and normal slopes given by 


w . 1 dw 
3 = f (Pl — Am)ds + ax+ By+y, 3 — = Al+ Bm+al+ Pm, (6.5) 


dn 


at each point of each boundary; 
(iii) There are a resultant force (F.) and resultant couples (M,) and (M,), on each 
boundary, with magnitudes given by 


(1 — »)F, Exc dV ; dT 7 
—_ = — — - (1 — 2n) $ -ds + Fea 6 = ds, (7.2) 
KD i+” dn dn 


(1 — vyi)yM. 


E,(a + Voc) : dV dV 
_- — —" +(1—-— 20) f y — — x — J ds 
KD 1 + vr dn ds 
dT rss ie 
-+ Brea 6 (» —-—vxX -) ds 
dn ds 
— [(1 — »)(1 — ») — 1] g (Y — Vm)ds, (7.13) 
(1—»)M, E\(b — xoc) dV dV 
— = —- - — (1 — 2») y —+ x—]ds 
KD i+ y ds dn 
dT dT 
= Bra 6 (» - 2 —)ds 
ds dn 


+ [(1 — )(1 — 2) - lg (X — Vids. (7.20) 
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ON AN EXTENSION OF THE VON KARMAN-TSIEN METHOD 
TO TWO-DIMENSIONAL SUBSONIC FLOWS WITH 
CIRCULATION AROUND CLOSED PROFILES* 


BY 
C, <.. 25 


Brown University 


1. Introduction. The method for treating compressible flows, as developed by 
Chaplygin,' von Karman’ and Tsien? leads to a successful solution of the flow pattern 
past solid bodies when the flow has no circulation. When the flow has a finite circula- 
tion, as in the case of airfoils, the profile shapes furnished by this theory are not closed. 
It is doubtless desirable to develop the theory so as to remove this difficulty. 

Recently, Bers‘ succeeded in obtaining flows with circulation around closed pro- 
files. As is usual in the case of a first success, the new method has a few disadvantages. 
In the first place, the mapping between the actual compressible flow and the asso- 
ciated incompressible flow is not regular at the stagnation points. Thus, if the profile 
in the associated incompressible flow is regular everywhere, angular points would 
appear in the profile in the compressible case; and vice versa. The application of the 
method is further complicated by the fact that the angle thus generated depends on 
the free-stream Mach number. For the engineer, the treatment has the additional 
inconvenience of involving the concepts of Riemannian geometry (which are avoided 
in the present treatment). 

In the present article we shall describe a method which is free from the disad- 
vantages mentioned above. The derivation is very simple, and no reference is made to 
Riemannian geometry. Yet the result includes all the previous ones as special cases. 
Indeed, the treatment seems to be now in the most natural and the most general 
form which is obtainable from the line of study of Chaplygin, von Karman and Tsien. 
It also has great flexibility. Given one incompressible flow, there is still an analytic 
function at our disposal for constructing compressible flows. This freedom of choice 
enables us to avoid much unnecessary numerical labor in constructing flows of certain 
general types. A large number of compressible flows can be derived from a given in- 
compressible flow by the present method without numerical integration. 

Apart from giving a useful method for constructing compressible flow patterns, 
the present development has the following significance. First, the freedom of dispos- 
ing of one analytic function leads to the solution of the direct problem,—namely to 

* Received May 18, 1946. 

1S. A. Chaplygin, On gas jets. Scientific Memoirs, Moscow Univ., Math. Phys. Sec. 21, 1-121 (1902). 
(English translation published by Brown University, 1944. Also NACA TM No. 1063, 1944.) 

2 Th. von Kérman, Compressibility effects in aerodynamics. Jour. Aero. Sci. 8, 337-356 (1941). 

3 Hsue-Shen Tsien, Two-dimensional subsonic flow of compressible fluids. Jour. Aero. Sci. 6, 339-407 

(1939). ; 
‘L. Bers, On a method of consiructing two-dimensional subsonic compressible flows around closed pro- 
files, NACA TN No. 969 (1945); On the circulatory subsonic flow of a compressible fluid past a circular 
cylinder, NACA TN No. 970 (1945). See also S. Bergman, On two-dimensional flows of compressible fluids, 
NACA TN No. 972 (1945). 
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calculate the compressible flow past a given profile. Indeed, the solution of problems of 
compressible flow, either direct or inverse,—namely, the construction of flows around 
profiles either given beforehand or not—1s now an a parallel footing with the incompres- 
sible case. In either case, the direct problem requires a method of successive approxi- 
mations.* Secondly, the application of the pressure coefficient formula of von Karman 
and Tsien to flows with circulation is justified on the same basis as in the circulation- 
free case. Experimentally, the formula has been found to be successful even when 
there is circulation, although the theory has so far been incomplete. The original de- 
velopment of von Karman and Tsien leaves the body not closed, while the profiles 
given by the method of Bers do not have the same analytic nature in the incom- 
pressible and the compressible cases. The present method removes these difficulties. 

Further investigations of the significance of the present method are being carried 
out. The present article contains only the essential result and its proof. It is hoped that 
a complete discussion of further developments may be published very soon. 

2. Method of constructing two-dimensional subsonic flows with circulation 
around profiles. Let p, p, u, v be the pressure, density and components of velocity of 
a steady two-dimensional irrotational flow in the x, y plane. Let p be a function of the 
density p only (given either by the isentropic relation or any other approximate rela- 
tion). Then there exist the velocity potential @ and the stream function y defined by 
the following differential relations: 


do 
dy 


The velocity components u, v and the density p are further connected by Bernoulli’s 


udx + vdy, (2.1) 
— pdx + pudy. (2.2) 


equation 
2 d 
- +f = const., (q? = u? + 0°). (2.3) 
2 p 


It is convenient to refer the density of the gas to that at the stagnation point and to 
refer all the velocities to the velocity of sound at stagnation. The coordinates x, y 
may be regarded as referred to the size of the body, and the pressure as referred to 
the product of stagnation density and the square of the velocity of sound at stagna- 
tion. Throughout this article, this process shall be implied, and all the quantities 
under discussion are dimensionless. 

As is well-known, the problem simplifies greatly if the pressure-density relation is 


approximated by 


incl —- (2.4) 


This is the basis of the method of Chaplygin, von Ka4rm4n and Tsien. A discussion 
of its physical interpretation may be found in the papers of these authors. Equation 
(2.3) leads to 

* For the incompressible case, see T. Theodorsen, Theory of wing sections of arbitrary shapes, NACA 


Rep. No. 411 (1931). 
T. Theodorsen and I. E. Garrick, General potential theory of arbitrary wing sections, NACA Rep. No. 


452 (1933). 
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cm — a Vi t¥ (2.5) 
p 
where c is the local velocity of sound. Indeed, B must be equal to unity if the reference 
stagnation quantities are calculated from (2.3) by the use of (2.4) 

Under the approximation (2.4), the following method may be used for construct- 
ing two-dimensional subsonic flows with circulation around closed profiles. 

Given an incompressible flow past a profile Py in the {-plane ({ =&+in) described 
by the complex stream function F({) and the complex velocity wo({), choose a function 
k(¢), regular in the region exterior to Py and including the point at infinity, having no 
root in Ro, and such that Ro 


| Swf) | < | k(f)| << © on Po, (2.6) 


a ae 
g k(¢)dt — cf = dt = 0, (2.7) 


where the integration is performed along any contour enclosing Po. Then 


and that 














1 fu) 
x+iy= f k(¢)dt — —f oe dg, (2.8) 
2 w 
— eee A (2.9) 
1+V/1+¢ k(¢) 
+ w = F(e) (2.10) 


gives the parametric representation of a compressible flow past a profile P in the x, y plane 
with ¢ as parameter, where P has the same general analytic nature (e.g. same number 
of corners, etc.) as the original profile Po. In these formulae, ,,y are the velocity potential 
and stream function defined by (2.1) and q, 6 denote the magnitude and direction of the 
velocity of the compressible flow. 

3. Proof. The proof consists of two parts. First, it is necessary to show that, after 
the auxiliary variable ¢ is eliminated from (2.8)—(2.10), we obtain proper functional 
relations between ¢, y, q, 8 and x, y. Secondly, we must show that the profile P is a 
closed curve and is mapped into P» by a regular mapping such that the regions R, ex- 
terior to P, and Ro, exterior to Po, are mapped into each other in a one-to-one manner. 

(a) The first part of the proof is simple. It is well known* that under the approxi- 
mation (2.4), the relations 


dF _ ’ 
isa —-— SP, Fo=PK er) (3.1) 
w* 4 
with z, F, w* defined by 
z= x-+ iy, ) 
F = | 
¢ + wy, (3.2) 
e=- — . e~ 8 | 


® Cf., for example, Eqs. (23), (26) of Tsien's paper quoted in Footnote 3. 
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would give a solution of the differential equations of compressible flow. Instead of 
trying to establish a relation between F and w* directly, we introduce an auxiliary 
variable ¢ such that F(¢) and w*(¢) are analytic functoins. It is well known that great 
simplification is obtained by taking F(¢) to be the complex potential of an incom- 
pressible flow similar to the compressible flow we desire. However, the extent of 
arbitariness in the choice of w*(f) has not been carefully examined. It is clear that 
any choice of w*({) will be sufficient so far as satisfying the differential equations is 
concerned. 
In the present case we dispose of the arbitrary function by writing 


w*(C) = wo(t)/R(¢). (3.3) 


The only requirements on k(¢) are the general conditions of regularity and the rela- 
tions (2.6) and (2.7), which will be discussed immediately. 

(b) The second part of the proof is also very simple. In the first place, the profile 
P is closed by virtue of (2.7).6 The regularity of the mapping is established if the 
Jacobian of the transformation maintains the same sign and does not vanish or be- 
come infinite in the region Ro, including the boundary Pp» and the point at infinity. 
It can be easily verified that the Jacobian is 


) x, y) 1 | wo '4 
J= nh a = Im\22%,; =| k 2) ct ae | ‘ (3.4) 
a(é, n) (16k! S 


From this expression, it is clear that the requirement is satisfied when & satisfies the 
restrictions specified in the last section.* 

4. Discussion. (a) The function k(¢). To make use of the freedom in choosing the 
function k(¢) is the essential improvement made in the present paper. Von Karman 
and Tsien gave an interpretation of w*(¢) by identifying it with the complex velocity 
wo() in the associated incompressible flow. This means that they put 


k(¢) = 1. (4.1) 


They were therefore unable to meet the requirement (2.7) for closing the profile, for 
flows with circulation. Bers overcame this difficulty by virtually taking 

k(¢) = const. } wo(¢) { '-Y'%, (4.2) 
where M,, is the free-stream Mach number. However, at the stagnation points of the 
incompressible flow 


| k(¢) | = 0, (4. 3) 


and the mapping between the profiles P and Pp is not regular there. In applying his 
method to the calculation of compressible flow past a circle, Bers had to start with 
the rather complicated problem of finding the incompressible flow past two intersect- 
ing circular arcs with the stagnation points at the points of intersection and with the 

* Note added in proof: Professor K. O. Friedricks pointed out to the author that, in the compressi- 
ble as well as the incompressible case, a mapping with non-vanishing Jacobian does not always lead to 
a useful result: the region obtained may be simply-connected but multiply covered. This difficulty has 
not been experienced so far in some numerical examples which have been worked out. 
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angles properly adjusted in relation to the free-stream Mach number. Indeed, it seems 
that after reaching the relation 





w*dF, 


the most natural development of the Kérman-Tsien method is to leave w*({) quite 
free, as we have done here, instead of connecting it definitely with wo, as was done by 
previous authors. The present method of leaving k(¢) free seems to be the most gen- 
eral scheme. 

If we deliberately want to introduce some singular points in P by starting from a 
profile Py without a singular point, | &(¢)| should be allowed to take the limiting val- 
ues specified in (2.6); e.g., k(€) =O where wo(f) =0. 

Although k(¢) cannot be taken to be unity when the flow has a circulation, it 
should not depart from unity very much if the profiles P and Po are not to differ 
much from each other. This is easily seen from a comparison of terms in (cf. (2.8)) 


1 wW 


dz = kdt — — — df. (4.4) 
4 k 


They have the ratio (cf. (2.9) and (2.5)) 


qg ; M? ” 
1¢4fit@/) (1+V/1— M?? 


where M is the local Mach number q/c. This value is much smaller than unity, except 
for values of M close to unity.’ Hence, (4.4) is approximately the identity transforma- 
tion if k is very close to unity. This approximately preserves the shape of the profile 
during the transformation. 

(b) Conformal mapping of compressible flows. lf we make a conformal transforma- 
tion of the ¢-plane into the ¢-plane by the analytic relation 


c= h(%), (4.6) 


we are merely making a change of the auxiliary variable in (2.8)—(2.10). Indeed, we 


have 
_. 1 pat) ; 
x+iy= f kdé — —f =~ af, (4.7) 
4 k(¢) 
2 Wo( ¢. 
Rete: GW oh tees (4.8) 
i+vVJi+¢ k(¢) 
+ ip = F(E); (4.9) 
where 
F(t) = F(o), (4.10) 
~ ar . d¢ 
wo(C) = F(t) = F’(¢) —, (4.11) 
dt 


7 Cf., Fig. 3 of Tsien’s paper quoted in Footnote 3. 
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and 


me d¢ 
¢) = k(t) — (4.12) 


(; 
dt 

We note that the equations (4.7)—(4.9) are of the same nature as (2.8)—(2.10). How- 
ever, the profile Po, into which the profiles P and Py are mapped, may bear no re- 
semblance at all to the original profiles. Indeed, there is no loss in generality in taking 
Py to be a circle. The relations (4:7)—(4.9) thus serve to transform the incompressible 
flow past a circle into a compressible flow past a profile of a quite arbitrary shape. Re- 
ferring to (4.6) and (4.12) and to the fact that k(¢) should be taken not far from 
unity, we see that &(¢) should be chosen so that it is not very much different from the 
derivative of the function mapping P into a circle. 

(c) Formulation of the direct problem. If we disregard the intermediate step of the 
¢-plane and drop the tilde, we have a mapping of the nature described in Section 2, 
but with k(¢) so chosen that a circle will be mapped into some profile P. The function 
k(€) must satisfy the requirements established there, but it should not be very much 
different from the function ko(¢) which maps the circle into Py by the relation 


c= f balsas. 


Thus for each profile P, the determination of the compressible flow past it is 
equivalent to the determination of a proper k({) mapping it into a circle by (2.8), 
where wo({) is the flow past the circle. There is no question about the existence of 
such a mapping function. It is clear that to each purely subsonic flow, where the ap- 
proximation (2.4) is accurate enough and therefore F = F(w*), we may find a certain 
k,(€) mapping the compressible flow into some incompressible flow. By considering 
successive conformal transformations, we can therefore always map the flow into a 
circle. The actual method of finding k(¢) is all that remains to be done in the direct 
problem. 

The theory of subsonic compressible flows (so long as the approximation (2.4) is 
valid) is now put on an equal footing with the incompressible flows. The inverse 
problem is complete, the direct problem of finding a mapping function for a given 
profile can only be solved (practically) by successive approximations, even in the in- 
compressible case. To develop a method of successive approximations for the direct 
problem in the compressible case seems to be a natural next step.* 

5. Application of von Mises’ method of generating airfoils. In the incompressible 
case, von Mises transforms a circle in the ¢-plane into an airfoil of a quite general 
shape by the transformation 


dz Ao Ai An - 
dé ft ¢ f 


8 Cf., e.g. the papers quoted in Footnote 5. 

* Note added in proof: The essential difference of the two cases lies in the existence problem. While 
the existence of the incompressible flow follows from well-known results concerning the Laplace equation, 
yery little seems to be known about the existence problem for compressible flows, 
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where Ay, - - - , A, are points inside the circle, and Xo is a point on the boundary (which 
transforms into the trailing edge of the airfoil). A similar method can be used here. 
In (2.8)—(2.10), we put 


Xo Ai An 
f f f 


with the same general restrictions on the points A. The condition that the points \ 
are inside the circle is exactly the condition required of k(¢) as stated in Section 2. 

The condition (2.7) for the closure of the body has also its equivalence in the case 
considered by von Mises. There it is® 


f Hous = 0. 


Here, it differs by the presence of another term. The condition that 
| R(S) | > | 2wo() | 


on the circle is the only additional restriction in the present case. As it is a mere 
inequality, there is no great difficulty in ensuring it to be satisfied. 

The integration required in establishing the relation between z and ¢ can be readily 
performed, as it involves only rational functions. This ease of calculation, together 
with the flexibility of the choice of the points \ in controlling the shape of the airfoil, 
are the advantages of the method of von Mises which are still preserved in the present 
application. 

The author is greatly indebted to Dr. J. B. Diaz for very helpful discussions in 
the course of the investigation and to Professors W. Prager and K. O. Friedrichs for 


their interest and discussions.* 


° Cf. W. F, Durand, Aerodynamic theory, vol. 2, J. Springer, Berlin, 1933, p. 78, Eq. (20.4). 

* Note added in proof: After the paper was completed and presented in a colloquium at Brown 
University, Professor E. Reissner informed the author that Professor A. Gelbart had recently presented 
a somewhat similar approach in a lecture at M.I.T. 








—NOTES— 


ON THE ELASTIC DISTORTION OF A CYLINDRICAL HOLE 
BY A LOCALISED HYDROSTATIC PRESSURE* 


By C. J. TRANTER (Military College of Science, Shrivenham, England) 


When a hydrostatic pressure is applied over only a small part of the length of a 
cylindrical hole extending through an infinite elastic solid, the stresses and displace- 
ments differ considerably from those caused by the application of this pressure over 
the entire length of the hole. This problem has been discussed by H. M. Westergaard! 
using an approximate method but it is not easy to assess the accuracy of his numerical 
results. It is the purpose of the present note to give an exact solution and to compare 
numerical results with those given by Westergaard. 

The analysis used here is a simple adaptation of that given by A. W. Rankin? 
for the similar problem of a band of uniform pressure applied to a long cylindrical 
shaft. The numerical calculations are not so formidable as would appear at first sight 
and a method given by L. N. G. Filon*® for evaluating trigonometric integrals has 
proved very valuable in this connection. The results for the maximum radial dis- 
placement show that the approximation used by Westergaard is rather crude. 

1. The analytical solution. We use cylindrical coordinates and consider the pres- 
|z| >c on the surface of the 


. . * | 
sure loading as being given by o,= —p, |2| <c, o,=0, 
cylindrical hole r=a. With the usual notation‘ we therefore require to find a stress 


function ¢ satisfying 


Vig = 0, r>a, — oOo < gs < oO, (1) 
where V? denotes 0?/dr?+ (1/r)(0/9r) +0?/dz? and the boundary conditions 
te) f 0" | 
o, = —4rV? — —7o = — , fi=—<c } 
Oz ] ar2t > y= a, (2) 
= Q, ee a 
re) J 0" | 
Tre = —¥(1 —v)V? - -->@=0, —x2m <3s< ©, r=—a4, (3) 
ar | dz? 


v being Poisson’s ratio for the material of the elastic solid. 

Once @¢ has been found, the stresses ¢,, T,, are given by the expressions shown in 

(2) and (3) and the remaining stresses and the radial displacement are given by 
af 1 8) af ay ity da 

; ?, = > (2-9) — Ta ee. : ’ 


r or f Oz 02?) E Oras 


(4) 


E being the modulus of elasticity. 
* Received May 8, 1946. 
1H. M. Westergaard, Kaérmd4n Anniversary Volume, 1941, p. 154. 
2A. W. Rankin, Shrink-fit stresses and deformations, Journ. Appl. Mech. 11, A77 (1944). 
3L. N. G. Filon, On a quadrature formula for trigonometric integrals, Proc. Roy. Soc. Edin. 49, 38 
(1928-29). 
4S. Timoshenko, Theory of elasticity, McGraw-Hill Book Co., New York, 1934, p. 309. 
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Following Rankin, we take 
wo 


o = f Rb sin cd sin zAdd (5) 


where R is a function of 7 only and 6 is a function of X. To satisfy (1) we must have 


and the solution, finite as r—> «, of this equation is 
R = AKo(p) + BpKi(p), (6) 
where Ko(p), Ai(p) are Bessel functions of imaginary argument, A and B are constants 


to be found and 


p = dr. (7) 
Using the well known relations 
Ké (p) = — Kilp), 
Daas (p - | (8) 
pKi(p) + Kil(p) = — pKol(p), 
we find 
V%> = — 2f BK,(p)bd? sin cd sin 2ddX. (9) 
With 
a = da, (10) 


substitution from (9) into (3) yields 


(tiie -{ |— BaKo(a) + }2B(1 — v) — A} Ki(a) Jd? sin cd sin sddd, 


0 


so that, to satisfy (3) 


A B = 2(1 — 7) — akKo(a)/Ky(a@). (11) 


We also find 


‘ A 
(0;)rxa = — f E + (2y — 1)B} Ko(a) + ( + Ba) Ka) | bd* sin cA cos sddX, 
J a 


and since the boundary condition (2) can be represented by 
— 
2p sin cA cos 2A 


(6,)rna = —— | —— ad, 
Tv 0 A 


we have 


2 A 
b = p E + (2v — 1)B} Ko(a) + ( + Ba) Kia) | ; (12) 
rr! a 
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Equations (11) and (12) yield 


. 4 Ve “a ) 
m\4D(a)bA 2p12(1 — v)aKy(a) — a?Ko(a) |, 
(a) p([2( aK 1(a) o(a) | (13) 


t\MD(a)bB = 2paK;,(a), 
where 
D(a) = fa + 21 -— v) | Ki(a) - a Kea). (14) 


bA, bB having now been found, the expressions for the stresses and radial displace- 
ment are found to be given by 


2pa = > ae + : : 2 
+ = =f [apKo(a) Ko(p) + aKo(a)Ki(p) — pKo(p) Ki(a) 
us 0 


c Zz } 








sin — a cos a 
a a 
— jp? + Ai —-— v)| Ki(a)K “1(p) J = da 
aD(a) 
Cc z | 
sin — a sin —a 
2p 7 . : : a a 
Tre = — " [aK (a) Ki() — pKo(p)Ki(a) | —— nnn See | 
Ww/o Dla ) i 
2pa " ere : 
Renee f [aK (a) Ki() + (2 — 1)pKo(p)Ki(a) 
Tr 
0 Cc « 
sin — a cos —a 
1 Kil Ki a a j ra (15) 
eo wy )K1(p) | —— ae | 
( 1(a) Ki(p) aD(a) oe | 
2p f° | 
c= ~~ f [aK ola) Kulp) + 2Ku(0) Kila) 
T 0 
c Z 
sin — @ Cos a 
K ( \K a a j 
~ Gimiig(————————— da 
1 D(a) | 
Eu 2pa 
... ao f [aK o(a)Ki(p) — pK o()Ki(a) | 
i+p us 0 } 
. c 2 
sin — a cos —a 
a 
— 2(1 — v)Ki(a)Ki(p) | — accede | 
it i(o) | aD(a) F ) 


2. Numerical results for the maximum radial displacement at r=a. When r=a, 
p=a and the greatest displacement occurs when z=0, so that we have 





ps | a 7 4pa(1 — —? se c 
1+ y - 
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If the pressure p acts over the entire length of the hole, the displacement (u’),. is 


given by 








E(t’) rama \rma 
1+, +y — 
so that 
max/r=a 4 1 -_ p 
(ss ) on ( —"f Ki(@) sin whe (16) 
(wens aD(a) a 


The numerical work was performed with v=0.3 and, above a=12, it was found 
that the first three terms of the asymptotic expansion of Kj(a)/aD(a), viz., 


Kia) 1 04 0.965 


—_—  /~ 
aD(a) a? aé at 





gave an adequate representation. Integration by parts then leads to 


” Ki(a) _ ¢ oy . id c 12c 
{ sin ada = | .08176 + .01340 — | sin —— — .01778 — cos —— 
Ji2 aD(a) a a? a a a 


cig , (12 c 12c¢ 
+ .2—|—- Si —) - E + .16083 — | — Ci ' 
a*>L2 a a’ a 
es * sin x sn * cos x 
Si (2) -f _ : dx, Ci(x) = -f dx. 
0 - x d 


The evaluation of the integral in equation (16) from a=0 to a=12 was performed as 
follows. The function Kj(a)/aD(a) was computed at intervals of a=0.2 from a=0 
to a=2 and at intervals of a=0.5 from a=2 to a=12. The integral was then evalu- 
ated by a method due to Filon® in which Simpson’s rule is replaced by the formula 


where 








B 
f F(x) sin kxdx = h[a{F(A) cos kA — F(B) cos kB} + BSxe + yS2e-1], 
{ 


where the range of integration is divided into intervals of length h, S2, is the sum of 
all the even ordinates of the curve y= F(x) sin kx between A and B inclusive less half 
the first and last ordinates, S2,_; is the sum of all the odd ordinates, and a, 8, y are 
given in terms of Y=hk by 
sin y cos p 2 sin? y 7 1+ cos*y 2sinycosy 
- ae j B= a ae Y 
y ¥° 








ce —— ; “ 


v v v' 


siny cosy 
. i i| roe | 
y? y? 





5 Filon, loc. cit. 
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This formula holds even when & is large, provided that the function F(x) can be 
fitted with reasonable accuracy over the range 2h by parabolic arcs. 


To avoid an infinity at the origin, the integral actually evaluated was 


ss l K (a) . Cc 
J = f | — sin ada, 
0 l.4da aD(a) a 


and when this had been found, the required integral was given by 


1 12c 
si(=*) = 1 
1.4 a 

Asa check that the substitution of the asymptotic series did not lead to unaccepta- 
ble errors, the range of integration was also divided into 0 to 10, 10 to infinity and the 
infinite integral was similarly computed on this basis. Little extra work was involved 
and excellent agreement was obtained. 

The results are shown below, together with those given by the approximate analy- 
sis by Westergaard. It is seen that even his second approximation is quite crude. 


Values of (tmax)roc/(t#’) rae 


Westergaard 


. j Present Method 
First Approximation | Second Approximation 
0.25 0.557 0.537 0.450 
0.50 


0.806 0.770 ' 0.633 


ON THE REPEATED INTEGRALS OF BESSEL FUNCTIONS* 
By J. C. JAEGER (University of Tasman‘a) 
It is well known that 


iJ ab)? 
LY : f = [(p2 + 1)'/? — p|", n>, (1) 


and 
_ ert - pb 


L}J,(t) | — 
(p? + 1)'? 


* Received Jan. 25, 1946. 
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where Lj f(t); is written for the Laplace transform of f(t), that is, 


© 


Li f(t); -f e*'f(t)dt. (3) 


These results have important applications in the theory of the semi-infinite dis- 
sipationless artificial transmission line with simple terminations, and thus in the ex- 
pression of the solutions of corresponding problems on finite lines in terms of multiply 
reflected waves. 

In an important class of similar problems in which the line is terminated by a 
matching resistance, the Laplace transforms of the solutions contain powers of 
[1+(p2+1)"?] or [p+1+(p?+1)"/?] in the denominator, and the functions which 
have such Laplace transforms do not seem to have been given. The object of this note 
is to show that they can be expressed in terms of repeated integrals of Bessel functions 
and that numerical values of these can readily be obtained. 


We use the notation 


(r), : ‘J, (td 
It @& -{ dt--- f : , n > 0, 
0 0 l 
t t 
Jt, Ai) = f dt-+-- f J,(t)dt, n 


for the r-ple integrals of J,(t)/t and J,(t) respectively. 
It is convenient to use both these types of integral though there are many relations 


(4) 


IV 


between them, the simplest being 

Fina (t) + Tinsr. A(t) = WnJin’ (2) (5) 
and 

Jin-1(t) — Jingir(t) = 2in s(t), (6) 


which follow immediately from the recurrence formulae for J,(t). Jto,(t) is tabulated! 
and Ji,” (t)=(1/n)+Ji,(t) where Ji,(t) is the ordinary Bessel integral function. For 
all values of » and r repeated application of the result 


f J, (t)dt = 2 Tnsomr(t) 


m=0 


gives the formulae 


pw. (m+r—1 
Jtg.AD = 2 ef ) Insane ( 


r=0 m 


~ 
— 


Pat yr! “ an ae 
nJi, (t) = 2 fe ee 3 


)Jnstmieald- (8) 
m=0 (m + Y= 1) 


m 


For integral values of f, which are in fact close enough for many practical purposes, 
(7) and (8) may be evaluated rapidly from the Tables in Gray and Mathews.’ 


1 Lowan and Abramowitz, J. Math. and Phys., 22, 2 (1943). 
2 Gray and Mathews, Treatise on Bessel functions, 2nd ed., 1922, Table II. 
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The Laplace transforms referred to above may now be written down. Firstly we 
have immediately from (1) and (2) 
—T 2 1/2 n 
p le +1) - 2] (9) 
[(o? + 1)? — p]” ‘aes 
rr + 1” 


L{nJix ()} 





ll 


L{ Jin (t)} 


Then, since 


1 -(1 1, “) 1 1 
1+ (p? + 1)!” p?/ (p? + 1)? p? 


it follows that 
1 


1+ (p+ 1) " 


L{Jo(t) + Jio.o(t) — t} = 


In the same way if »>0 
P 2 “p- 1 2 o. n 
L\Jn(t) + Jin 2(t) — nJix(t)} = Lib ea. - (12) 
1+ (9? + 1)" 


Similarly 
L{Tn(t) + 3Tino(t) + Win a(t) — 2nJix (1) — 2nJi, (t)} 
Fp ie es. se te alla: 
1+Q@+ynee 


if n>0, and if »=0 the term nJi\” (t) is to be replaced by ¢*-1/(r —1)!. Again with this 
I Z g 





convention we have 


e+" = 2F 





L{anJic (2) — 3(n+ 1) Fics s(t) = - ; (14) 
. Olmeit@e ym 
BHiadd — Ueadkd +I + Haudd ~ Mind + Ji) 
4(p? + 1)1/2[(p? + 1)1/2 — p|” a 
-— ————-- (15) 





[p+ 1+ (9? + 1)"?]? 


These expressions may be transformed in many ways using (5) and (6) and gen- 
eral results for higher powers in the denominators* may be obtained in the same way. 

As an example of the way in which the above functions arise, we consider a semi- 
infinite artificial transmission line with mid-series termination, in which the series 
elements are inductances Z and the shunt elements are condensers of capacity C. 
Suppose that all condensers are charged to unit potential, and that at time t=0 the 
line is discharged through the matching resistance \/(L/C). Then if Jo is the current 
in the resistance, J, that in the mth inductance L, and Cv, is the charge on the mth 
condenser, applying the Laplace transformation method in the usual way we find that 


2 Wilt... 5 2r 
aC [(1 + p?/a*) bet ll satlidies's (16) 
2p[1 + (1 + p*/a’)!/?] 


2 The extension of (15) is trivial; for that of (13) the results needed are given in Chrystal, Textbook 
of algebra, 2nd ed., 1906, vol. 2, pp. 204—205. 


j == 
j 
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1 [(1 + p?/a2)'? — p/a)? 


L\v,} =- ee r=1,2,---, (17) 
p p[1 + (1 + p?/a?)'/?] 
where a=2(LC)-"/2. 
It follows from (12) that 

v7, = 1 _ J ty 1,1(at) —_ J is, 1,3(at) + (2r mcs 1) Jiz,_(at) (18) 
I, = (C/L)"!2| Sie (at) + Jire.s(at) — 2rJiz, (at)} (19) 
To = (C/L)"!?| Jig.(at) + Jio,3(at) — 4a7t?} (20) 

= 4(C/L)"24 (1 + a°t?)Jio (at) — at?(1 + J,(at)) + atJo(at)}, (21) 


where (21) follows from (20) by integration by parts. 


If the line is discharged into inductance }Z and resistance \/(L/C) in series, the 


solution follows from (14) in place of (12). 


ON CERTAIN INTEGRALS IN THE THEORY OF 
HEAT CONDUCTION* 


By STEWART PATERSON (.C.J. (Explosives) Limited, Stevenston, Scotland) 


In a recent note' W. Horenstein evaluates the integrals 


t 1° 

o= f x-3/? exp (- on in) dx, (1) 
P x 
¢ a" 

y= f a—'/? exp (- —-— ix) dx, (2) 
0 x 


in terms of the tabulated exponential and error functions. The evaluation of the more 


| exp (— s* — n?/s*)ds 
U 


from which @ and y are easily derived, was given by Riemann.’ 

Integrals of the above type arise in the solution by classical methods of various 
heat conduction problems. It is the purpose of this note to point out that treatment 
of many such problems by the Heaviside “operational” or equivalent Laplace trans- 
form method leads directly and naturally to the required solution in tabulated func- 


general integral, viz. 


tions. 
Thus, to take a simple case, the classical solution of 
06 1 0°@ 
=— — b*0; @—0, t— 0, 0d 1, a 04,7" (3) 
at 4 da? 


* Received Nov. 24, 1945. 
1 W. Horenstein, Quart. Appl. Math. 3, 183--184 (1945). 
2 B. Riemann, Partielle Differentialgleichungen, 2nd ed., 1876, p. 173. 
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(where @ is a function of a and ?¢) will be 
0 = an!*g. 
If, however, 


A(a, p) ={ e~”6(a, t)dt, 
0 


equations (3) transform into 


076 n i 1 
— = 4(p + 6°)6; @——, a— 0+ (4) 

da? p 

which lead at once to 
§ = p- exp [— 2av/(p + 6%) ]. (5) 
The inversion theorem for the Laplace transform then gives 
1 vt io secuitiomaals 
6 = — A-' exp [At — 2av/(A + 8?) Jd, (6) 
2rid y~ix 


along the usual contour. 
By a series of obvious and natural steps,’ it is easy to show that this is equal to 
e2ab y'+i 


x e 2ab y/'’+ix 
— A-l exp [At — 2(a + bt)V/A] dd + re f A- exp [At — 2(a — dA) ddd, 


Amid + ix ct a 


e2ab a : e~20b a 
= | — erf (— + by ‘)] >> (1 — erf (< — by ‘)|. 
2 Vi 2 Vi 


and it can be verified that this satisfies (3). 


NOTE ON A FORMULA FOR THE SOLUTION OF AN 
ARBITRARY ANALYTIC EQUATION* 


By HERBERT E. SALZER (Mathematical Tables Project, New York City) 


In a recent note D. R. Blaskett and H. Schwerdtfeger! give a fairly well known 
expansion for a root a of the equation f(z) =0, as a power series in f(2o), where Zo is 


near a, namely, 
= (= 1)’ [f(20) ]’ Paf-(w) 
ey sesind Al fal Mo | (1) 
v=O vp! dw’ w= J (z) 
where w denotes f(z). 
Of use in connection with (1) is a paper by Van Orstrand, “Reversion of Power 
Series,” Phil. Mag., (6) 19, 366-376 (1910). Van Orstrand’s article deals with the re- 
3H. Jeffreys, Operational methods in mathematical physics, Cambridge, 1931, p. 70. 
* Received January 26, 1946. 
1 This Quarterly 3, 266-268 (1945) 


1946] H. B. MANN 307 


version of the special type of power series y =aox-+a,x?+a2x*+ --- to obtain x as 
an integral power series in z=y/do, whose coefficients are given as polynomials in 
b;=—a;/ao, as far as the term involving z*. Now the explicit expansion for (1) in 
terms of the derivatives of f(z) at zo can be written down immediately, as far as y= 13, 
from Van Orstrand’s expansion on pp. 369-370, merely by 

(A) replacing 0); in his formula by —f‘**” (g9)/(¢+1) !f’(z0), 

(B) replacing his 2 by —f(20)/f’(20), and 

(C) adding the constant term Zo. 

The truth of the last statement is obvious from the fact that when (1) is applied 
at the origin it yields Van Orstrand’s expansion and from the uniqueness of Van 
Orstrand’s expansion. 


A NOTE ON THE CORRECTION OF GEIGER 
MULLER COUNTER DATA* 


By H. B. MANN (Ohio State University) 


The correction of Geiger Miller Counter data has been considered in a previous 
paper by J. D. Kurbatov and the author.' According to the model described there 
the following result was proved: If the density of radiation is a constant a and if 7 
denotes the resolving time, B(7) the number of discharges during the time 7; then 


B(T) a + (1) 
a Up 

[+ 
where 7 is given by 


T 
n= of e(t)dt (2) 
0 


and ¢(¢) satishes the conditions 





t j 
et) = — of e(x)dx for t2r7, 
t—1 (3) 
ar 
e(t) = 1 — eo ** — ——— for 0 St Sr. 
1+ ar ) 
It was further shown that for ar <1. 
bi (ar)? 
| o| s ——> [It - (es), 
1 — (ar)? 


where s is the largest integer not larger than 7/r. In this paper an upper bound for 
|7| will be derived without the restriction ar<1. We shall prove the following in- 
equality: 





* Received May 29, 1946. 
1 J. D. Kurbatov and H. B. Mann, A correction for Geiger Miiller counter data, Phys. Rev. 68, 40-43 
(1945). 
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Let [x] denote the largest integer not larger than x; then 


- r j at a*r : i ee 
In(T)| Ss (2e% — 1) +- - T - ") (1 — e797) (Ter), (4) 
1+ ar 1+ ar T 


Proof of the inequality (4). From (3) we see that ¢€(¢) is a continuous function of ¢. 
Applying the mean value theorem to (3) we have 


e(t) = — are(i*), t—rsgf st. 
Hence e(t) changes its sign at least once in every open interval of length 7 and will 
therefore be 0 at least once in every such interval. Hence we have 
PROPOSITION 1. In every open interval of length r there is at least one point for which 
e(¢) =0. 


Differentiating (3) with respect to ¢ we obtain 


e'(t) = ae(t — r) — ae(t). (3) 


In the interval /<$tsi+7r Eq. (5) may be considered as a differential equation for 
e(¢t) with the initial condition that its solution be equal to e(é) at the point 7. Solving 
(5) with this initial condition we have, for /S¢<i+r, 


t 
e(t) = e(d)e~* "9 + ae «ff e"e(x — r)dx. (6) 
t 


Let M(i) be the maximum of the absolute value of ¢(¢) in the interval [7—7, 7], then 
e(t)| S M(d)e-*-© + eM (i)(e" — e) = M(t) for tstst+r. (7) 

From (7) it follows that | e(¢)| <M for t2% Hence we have 

PROPOSITION 2. If | e(t)| <M for t—r St Si, then | e(t)| <M fort=i-r. 


If e(?) =0 then we obtain from (6) 
(8) 


lA 
~ 
“~ 
“4 


et)| $< M(D(A — ec *?) S M((1—e°) for ists 


From (8) and Proposition 2 follows 

PROPOSITION 3. If e(4) =0 and | e(t) | <M for i—r sisi, then | e(t)| <= M(1 —e-*") 
for t>7. 

According to proposition (1) we have in the interval ar S¢S(a@+1)r at least one 


point ¢, for which e(¢,) =0. 
Consider the points /;, fs, - - 
we must have, according to Propositions 1, 2, and 3, 


-, tongs. If M is the maximum of | e(t)| in OStS<r 


M for 0sSt3Sh, 


IA 


| e(t) | 
e(t)| S M(1 —e-°") for ft, 


IA 
IIA 


lA 
lA 
IIA 


e(t)| S M(1 — e-*")* for tex 
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Hence 


T 


7 | | a=[7/r] ] 
| = a f e(x)dx | = >» e(ar) — af e(x)dx | 
J 9 | | (T/r]}r | 


i a=! 


7 
< M + 2M(1 — e"") + 2M(1 — 6°97)? + --- + oM( 7 | ") (1 — e927) (T/2r] 


T 
T +.” 
M(2e7 — 1) + aM (7 — | ’) (1 — e927) (T/?r1, 
T 


From (3) it can be seen that MJ =ar/(1+ar) and (4) follows. 
The inequality (4) is very satisfactory and shows that even for large values of a7 
the quantity 7 will be very small compared to a7/(1+ar) even if T is only a few 


minutes. 


WA 


CORRECTIONS TO OUR PAPER 


STABILITY OF COLUMNS AND STRINGS UNDER 
PERIODICALLY VARYING FORCES* 


QUARTERLY OF APPLIED MATHEMATICS, 3, 215-236 (1943) 
By S. LUBKIN anp J. J. STOKER (New York University) 


The following errors were found in the tables printed on pp. 232-235. 





| 

















( | _ >) l 
8 ta for read | ? iy ” for read | 
1.6 —0.77898 —0.77897 | | 3.6 — 2.32402 —2.32401 
1.8 —0.92281 —0.92282 | 
7.6 ~5.71537 —5.71538 
9.2 7.11974 —7.11975 | | 
8 ia for read | Pa for | read 
0.8 0.55906 0.55406 3.8 | —0.00468 —0.00464 
1.4 0.63015 0.63016 | 6.8 | —1.60383 | —1.60379 
4.4 -0. 29781 —0.29780 | | 8.4 | —2.58478 —2.58477 
7.6 —2.08644 | —2.08648 | | 
| 11.0 —4.29436 | —4.29437 | 
| ee = PS 1 a. 
; es for read Pama for read | 
2 Kelle tee SEM esol SES. ea eS 
0.6 1.12806 | 1.12810 | 0.6 | 2.26622 2.26621 
z 2.01478 2.01477 1.0 2.28515 | 2.28516 
20.0 —5.05198 | —5.05199 | | hee 2.31495 2.31493 | 
2.8 2.29660 | 2.29661 
| 5.6 | 1.85589 | 1.85591 





* Received Aug. 16, 1946. 
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a for read ns Neha for read 
S:,* = a 2a Se 
1.6 2.51308 2.51309 | | 1.6 | 4.07660 4.07659 | 
2.0 2.66777 2.66776 | 1.8 | 4.09433 4.09432 | 
3.8 3.46578 3.46579 | 2.4 | 4.15212 | 4.15211 | 
| 4.4 3.69216 3.69215 saan 4.17199 | 4.17200 | 
| 5.6 4.01149 4.01150 | | iZ.6 | 3.38817 | 3.38820 | 
| | 
ko | 
3 a(Co)| for read 8 SG) | for read | 
- * am a — a Ea ce 
1.6 | 4.09776 4.09777 | | 16.0 6.52721 | 6.52709 | 
2.4 | 4.24889 4.24891 | | | 
2.8 | 4.35867 4.35865 | | | 
; 48 (] 5.18127 5.18128 
| 6.0 | 5.74803 5.74303 | | | 
Column headings for a(S;) and a(Ss) are interchanged on p. 234. 
a rn a + | 
le we for read aw for read | 
halle.” EEER Roe ETS ase, & GRroetit Veale : 
| 0.4 | 6.25333 6.25334 3.8 | 9.20714 9.20713 
i ae 6.35487 6.35488 5.2 | 9.38281 9.38279 | 
es | 6.48591 6.48590 16.0 10.59848 10.59849 | 
|} 4.8 } 6.86185 6.86180 20.0 | 10.35813 10.35825 | 
: 32 6.99394 6.99396 | 
| 5.6 7.14093 7.14116 | 
20.0 10.33749 10.33744 











BOOK REVIEWS 


Theory of Structures. By S. Timoshenko and D. H. Young. McGraw-Hill Book Com- 
pany, Inc., New York and London, 1945. xiv+ 488 pp. $5.00. 


This valuable addition to text-book literature is based on the senior author’s earlier volume, pub- 





lished in Russia in 1926 (S. Timoshenko, Theory of Structures, Leningrad). The book is intended for engi- 
neering students with some background in mechanics. The keynote of this book is that familiarity with 
the general principles of mechanics is indispensable to a thorough understanding of the analysis of stresses 
in trusses and frames. For this reason two of the nine chapters are devoted to a comprehensive recapitula- 
tion of the rudiments of plane statics and of such general theorems on elastic systems as the Principle of 
Least Work, Castigliano’s Theorem, Maxwell's Reciprocal Theorem, etc. 

As one would expect from the authors, the book is very clearly written. It abounds in carefully con- 
structed figures and diagrams, and contains a wealth of well-graded problems. 

The chapter headings are as follows: Elements of Plane Statics, Statically Determinate Plane Trusses, 
Influence Lines, Statically Determinate Space Structures, General Theorems Relating to Elastic Sys- 
tems, Deflection of Pin-jointed Trusses, Statically Indeterminate Pin-jointed Trusses, Beams and Frames, 
Arches. 

This book will be of considerable interest to structural engineers and will be welcomed by the teachers 


of mechanics and theory of structures. 
I. S. SOKOLNIKOFF 
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Table of arc sin x. Prepared by the Mathematical Tables Project conducted under 
the sponsorship of the National Bureau of Standards. Official Sponsor: Lyman J. 
Briggs. Project Director: Arnold N. Lowan. Columbia University Press, New 
York, 1945. xix+121 pp. $3.50. 

The main tables give the values of arc sin x to twelve decimal places, the intervals of the argument 
being .0001 in the range betWeen 0 and .9890, and .00001 in the range between .98900 and unity. To 
facilitate interpolation the second (and, wherever necessary, the fourth) differences are tabulated, and 
auxiliary tables are given for the coefficients in the interpolation formulas of Newton-Gregory and 
Everett. For values of x exceeding 0.99950, interpolation by means of differences becomes unsatisfactory. 
For such values of x the use of the formula arc sin(1—v) =2/2—f(v)4/2v is recommended, and f(r) 
=1+v/12+3v2/160+5v3/896+ --- is tabulated (with first and second differences) to thirteen decimal 


places at intervals of .00001 in the range from 0 to 0.00050. 
W. PRaGER 


Network Analysis and Feedback Amplifier Design. By Hendrick W. Bode. D. Van 
Nostrand Company, Inc., New York, 1945, xii+551 pp. $7.50. 


This book is concerned with a complete exposition of electrical circuit theory, the properties and de- 
sign of feedback amplifiers, non-feedback amplifiers, and the discussion of certain problems of wide band 
transmission. A great deal of the material presented in this book has not appeared before in text book 
form. 

The book is divided into nineteen chapters. The first two chapters are devoted to the presentation 
of the fundamental principles of linear, passive, electrical circuits and to a formulation of the fundamental 
equations of these circuits from the mesh and nodal standpoints. The response of linear circuits to driving 
functions of the exponential type is considered and the very useful concept of the complex frequency plane 
in the study of the properties of linear circuits is introduced. 

In the next four chapters, the basic principles and theorems of feedback are considered in detail. A 
thorough discussion of stability, physical realizability, contour integration, Nyquist’s criterion for stabil- 
ity, and the physical representation of driving point impedance functions, occupies a central position in 
the book. The remaining chapters are devoted to the design of impedance functions, equalizers, inter- 
stage networks, single loop amplifiers, single loop feedback amplifiers, and a discussion of the relations 
between the real and imaginary components of network functions. 

From a mathematical standpoint, the material presented in this book is a beautiful example of the 
power and utility of the application of the fundamental theorems of the complex variable to a most im- 
portant physical problem. Since the subject of network analysis and synthesis is of such great importance 
not only in itself but also because it serves as a model for the analysis of mechanical and acoustical sys- 
tems, the excellent original analysis of the problem presented by Dr. Bode in this book is a great contribu- 


tion to the field of applied mathematics. 
Louts A. Pipes 
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